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X XIX. BAND SECHSTES (SCHLUSS-)HEFT 1960 
eee 


Kin neues Verfahren zur Drehschwingungsberechnung 
von Kolbenmaschinen bei Beriticksichtigung der Torsion zweiter Art* 


Von H. Bufler und F. KieBling 


1. Einfithrung. Zur Berechnung der Torsionseigenfrequenzen von Maschinenaggregaten sind 
zahlreiche verschiedenartige Verfahren bekannt!, denen die Annahme glatter Wellenstiicke zwi- 
schen den Drehmassen zugrunde liegt. Ein solches Ersatzbild ist bei den Kurbelwellen von Kolben- 
maschinen jedoch nicht ohne weiteres zuldssig, wie erstmals R. Grammel? erkannte. Danach hat 
man bei Kurbelwellen zwischen Torsion erster Art und Torsion zweiter Art zu unterscheiden. Die 
Torsion erster Art kommt durch allein an den Wellenenden eingeleitete Drehmomente zustande. 
Die Verformung der Kurbelwelle durch an den Kurbelzapfen angreifende Krafte bezeichnet man 
dagegen als Torsion zweiter Art. Da diese den wirklichen Beanspruchungsverhaltnissen bei Kolben- 
maschinen entspricht, hat eine genaue Schwingungsrechnung hiervon auszugehen. Im Unterschied 
zur Torsion erster Art hangt hier die Differenz der Drehwinkel zweier benachbarter Drehmassen 
nicht nur von dem zwischen ihnen wirkenden Drehmoment, sondern auch von allen iibrigen Dreh- 
momenten ab. Die Beitrige dieser iibrigen Drehmomente zur Gesamttorsion im betrachteten 
Wellenabschnitt werden als Nebentorsionen bezeichnet. 

R. Grammel® fiihrt die Schwingungsberechnung bei Zugrundelegung der Torsion zweiter Art 
mit Hilfe seiner Frequenzfunktionen durch und zeigt den Einflu8 der Nebentorsionen auf die 
Eigenfrequenzen. Insbesondere behandelt er eine homogene Maschine ohne bzw. mit einer Zusatz- 
drehmasse, wobei er die Nebentorsionen durch eine Stérungsrechnung beriicksichtigt. 

In der vorliegenden Arbeit wird unter Beschrankung auf die primaren Nebentorsionen (der 

-EinfluB der iibrigen ist vernachlassigbar klein) zunachst die Frequenzdeterminante aufgestellt 

und ein (schon bekanntes) tabellarisches Lésungsverfahren angegeben. AnschlieBend wird der 
praktisch wichtige Fall der homogenen Maschine ohne und mit Zusatzdrehmassen mit Hilfe der 
trigonometrischen Methode behandelt. Wahrend danach im rein homogenen Fall. die Kigen- 
-frequenzen in geschlossener Form streng darstellbar sind, miissen diese beim Vorhandensein von 
Zusatzdrehmassen aus den Nullstellen transzendenter Gleichungen bestimmt werden. Ferner 
werden die entsprechenden Formeln fiir die Momente und Drehwinkel angegeben. 


2. Ausgangsgleichungen und Frequenzdeterminante. Das Momentengleichgewicht an der k-ten 
Scheibe liefert die Bewegungsgleichung 
M,— My_1 = OF. - (1) 
‘Dabei bedeuten O, das polare Massentragheitsmoment, gy, den Drehwinkel der k-ten Scheibe 
‘sowie M;, und M,_, die Momente in den benachbarten Wellenabschnitten gemaf Abb. 1. Mit 

dem Ansatz fiir stationare Schwingungen 


, = %, snot, 


Abb. 1, Allgemeines Drehschwingungssystem mit NN Drehmassen. 


* Die Verfasser sind Herrn Prof. Grammel fiir die Anregung zu dieser Arbeit sehr zu Dank verbunden. 

1 K. Klotter, Ing.-Arch. 17 aed =: ere 

2R.G I, Ing.-Arch. 4 (19 . 287. se Fa 
Es € a Sena R. Eehihel.  Fechntache Dynamik Bd. 2 S. 410, 2. Aufl. Berlin/Gottingen/Heidelberg 
"1953. 
: 26 
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wobei w die Eigenkreisfrequenz und t die Zeit sind, geht (1) iiber in die spater mehrmals bendtigte 
Gleichung 


M,— M,_;1 = — O, wg - (2) 


Da die Differenz der Drehwinkel zweier aufeinanderfolgender Scheiben bei der Torsion zweiter 
Art auBer von dem zwischen ihnen wirkenden Torsionsmoment auch von allen iibrigen Momenten 
abhangt, gilt folgende lineare Formanderungsbeziehung : 


Pei —Pu = 91,1 My + 4,2 Mo + +++ + Op,n—1 Mg_a + 5p,4 Mia + 5,,041 Mig +++ 0p,n-1My_i- (3) 


Die EinfluBzahlen 6,,;, welche wegen des komplizierten Aufbaus einer Kurbelwelle durch geeignete 
Versuche bestimmt werden miissen, sind durch diese Gleichung definiert. Aus dem Satz von Betti 
folgt, daB die Matrix der KinfluBzahlen 0, ; zur Diagonale symmetrisch ist, d. h. 


Oni = Oi,k 


t) 


Zum Beweis betrachten wir zwei Wellenabschnitte mit den Momenten M; und M,, wahrend alle 
iibrigen Momente Null sein sollen. Bringt man zuerst das Moment M; und dann das Moment M, 
auf, so leistet das erstere beim Aufbringen von M,, die Erganzungsarbeit 6; M, M;; beim Auf- 
bringen dieser Momente in umgekehrter Reihenfolge leistet M;,, die Erganzungsarbeit 0, ;M; M,, 
Nach dem Satz von Betti miissen beide einander gleich sein, woraus sofort 6, ; = 0;,, folgt. 

Man bezeichnet in (3) 6,,, als EinfluBzahl der Haupttorsion, 6;,,_; und 0,,,; als HinfluB- 
zahlen der primaren Nebentorsion, 6,5 und 6,449 als EinfluBzahlen der sekundaren Neben- 
torsion usw. Wie R.Grammel, K. Klotter und K. von Sanden zeigen, liegen die EinfluBzahlen 
der primaren Nebentorsionen eine GréBenordnung, diejenigen der sekundaren Nebentorsionen 
um zwei GréBenordnungen usw. tiefer als die Einflu8zahl der Haupttorsion. Aus diesem Grunde 
ist es gerechtfertigt, sich auf die Haupt- und die primare Nebentorsion zu beschranken. Damit 
geht (3) tiber in 


Pr41— Vr = On,n—-1 My_y + Op,n My + Op,4-41 Mays - (4) 
Eliminiert man nun die Drehwinkel aus (4) mittels (1), so erhalt man 
. - 5 M,_ 1 oe Mp1 
On,k—-1 Mu_1 + 9p, Ma + On,441 Miya 6, oo (5, ai Sa) M,— Scie 0. (5) 


Mit dem fiir stationare Schwingungen giiltigen Ansatz 
M, = M, sinw t 


ergibt sich aus (5) folgende Rekursionsformel fiir die Momente: 


1 2 = 1s 1 2 zs 
G + @ b1.1-1) My1 Ae oe Ont ) Ont] M, + c= +o ds.u+1] Mizi=0. (6) 


Bemerkenswerterweise kommen in dieser Formel — wie bei der Torsion erster Art? — nur drei 
aufeinanderfolgende Momente vor. 

Schreibt man diese Formel fiir k = 1, 2... N — 1 unter Beriicksichtigung der Randbedingungen 
M, = 0 und My = 0 (N Zahl der Drehmassen) an, so erhalt man ein lineares homogenes Gleichungs- 
system fiir die Drehmomente, welches nur dann eine nichttriviale Lésung besitzt, wenn die zu- 
gehérige Determinante der Koeffizienten [siehe (7)] verschwindet. Diese Bedingung liefert die 
Gleichung zur Bestimmung der Eigenfrequenzen. 

Da die Auflésung der Frequenzdeterminante bei einer gréBeren Anzahl von Drehmassen auferst 
unbequem ist, wird man in diesem Fall auf ein tabellarisches Verfahren zuriickgreifen, das eine 
Erweiterung der Methode von Holzer-Tolle auf die Torsion zweiter Art darstellt. Der Vollstandig- 
keit halber fiihren wir dieses Verfahren, das bereits von A. Kimmel und M. Ldapple? angegeben 
wurde, kurz an. 


Ausgangspunkt ist Gleichung (2). Schreibt man statt k den laufenden Index J und summiert 
zwischen | = 1] und | = k, so folgt 


k 
Ni ae O19 - (8) 
1 Ing.-Arch. 7 (1936) S. 439. 


* H. Neuber, Ing.-Arch. 22 (1954) S. 258. 
3 A. Kimmel und M. Lapple, Ing.-Arch. 12 (1941) S. 320. 
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Fiir k = N ergibt sich hieraus wegen My = 0 (Randbedingung) 
N a 
d' 9,9 = 9. (9) 
1=1 


Man setzt M,, M,_, und M,,, gemaB (8) in (4) ein, lést die so érhaltene Gleichung nach @, , ; auf 
und bekommt dann die fiir das Verfahren grundlegende Formel 
k 
(1 + wd, 5 On) HH — ? (Og, 4-1 + Oegk + Oe,e+3) o, 2! 
(Le Gene) Cece OID SA Top DES ee ¢ gail d Sa (10) 
Fi 1+ w O,,n+1 O41 
Wegen M, = My = 0 sind hier 6, 9 und Oy_4, y formal Null zu setzen. 

Fiir eine bekannte Eigenkreisfrequenz w gestattet (10) die Berechnung samtlicher Drehwinkel- 
verhaltnisse, und die Bedingung (9) ist mit diesen Drehwinkelverhaltnissen erfillt. Dagegen 
wird bei einer mit einem beliebigen Wert w (der nicht einer Eigenfrequenz entspricht) durch- 
gefiihrten Rechnung (9) nicht erfiillt sein, sondern einen bestimmten Restwert R liefern. Berechnet 
man auf diese Weise fiir eine gréBere Zahl von w-Werten den Ausdruck R und tragt diesen iiber w 
auf, so sind die Abszissen der Nullstellen der Kurve R = R(w), da die Bedingung (9) mit diesen 
erfiillt ist, die richtigen Eigenkreisfrequenzen. 

Zur praktischen Durchfiithrung des Verfahrens benutzt man am besten Tabellen nach dem 
Muster der Tabelle 1. Man nimmt fiir w einen Wert an und fillt zunachst die mit Kreuzen (x) 
bezeichneten Felder aus. Die erste Zeile der Spalte 12 liefert sofort ,/~ 1, das dann in Spalte 7 der 
zweiten Zeile benutzt wird und in Spalte 12 y3/p, ergibt, das in der dritten Zeile benétigt wird usw. 
SchlieBlich erhalt man in der N-ten Zeile der Spalte 7 den dem angenommenen w- Werte zugehérigen 
Restwert R und damit einen Punkt im R,@-Diagramm. Zur Erlauterung wird am Schluf ein 


Beispiel durchgerechnet. 


3. Die homogene Maschine. Die rein homogene Maschine ist dadurch gekennzeichnet, daf 
samtliche Drehmassen sowie EinfluBzahlen der Haupttorsion und der Nebentorsion untereinander 
gleich sind, d. h. 

0,=0; Op,k = 03 On,e—1 = Onnti = O* . (11) 
Wirkliche Maschinenaggregate bestehen aus einem homogenen Teil und einer oder mehreren zu- 
satzlichen Drehmassen, die von denjenigen im homogenen Teil verschieden sind, ebenso wie die 
entsprechenden Einflu8zahlen. In diesem Fall spricht man von homogenen Maschinen mit Zusatz- 
drehmassen. 

Mit (11) vereinfacht sich die Rekursionsformel (6) zu 

2— 0? 00 


M._1— ToS" 6 


M,+My..=0. (12) 


Diese Differenzengleichung besitzt die Lésung 
M, = Csin (kx + y)sinwt » (13) 
(C beliebige Konstante), falls man statt w die GréBe x gemaB 


2—a?bO 


2 cos x = Tt ots*6 (14) 

einfiihrt. 
Man kann durch Anschreiben vou (6) fiir k = 1,2 bis (n —1) feststellen, daB® bei der rein 
homogenen Maschine dafiir Gleichung (12) gilt, sofern man 6,,) = Onn—1 = O* setzt, was wegen 


der Randbedingungen M, = 0 und M, =0 ohne weiteres zulassig ist. Infolgedessen gilt (13) 
_ bei der rein homogenen Maschine fiir k = 0 bis k = n. 
Die Auflésung von (14) nach @ liefert 


5 x 
ee 2 (1 — cos x) fe ‘ psig 2 
SOTSeo coon ene 5 (15) 
( hes 0 +25 cosa] 
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Der Kigenwert x und die Konstante y i u i i 

y in (13) miissen dabei aus den Randbedingungen bestimmt 
werden. Den Winkelausschlag g, der k-ten Masse erhalt man nach kurzer Umformung aus (2) 
wenn man M,, M,_, nach (13) sowie w aus (15) einsetzt: : 


o* 
a(t + 2-5 cose 
—C cos 


eae 
2 — | 
sin > 


PY, = sinwt. (16) 


1 
bs}ers 


Um die Konstante C aus (13) und (16) zu entfernen, ist es zweckmaBig, das Momenten- bzw. 
Drehwinkelverhaltnis M,/M, und ¢,/, zu bilden*: 
M, _ sin(kx +) 
MM, sin(x+y) ’ 


Thiet («-4)» A | (18) 


m cos (+ a ”) 


Bei der rein homogenen Maschine gelten die Randbedingungen M, = M,, = 0. Daraus erhilt 
man in Verbindung mit (13) 


(17) 


=i) 19 
Re V3 (19) 
fe Sy sly 2am lye (20) 
Mit (20) geht (15) tiber in 
2 sin ™ 
an 
= 1,2...1— 1): (21) 


o, =— = 
Jee( ae 25 00S) 
) n 
Aus (17) und (18) erhalt man mit (19) und (20) das Momenten- bzw. Drehwinkelverhaltnis 
| 1 kya 


M, oe ey ieee alt 
(an), = pee peed, 2s ee G4) 


in) alta) Gee ) a8) 


Pi y=1,2...n—l1 


0, = (Oe —o(1 +56] (24) 


ab. Hierbei bedeutet (cw,)s+ — 9 die Eigenkreisfrequenz bei alleiniger Beriicksichtigung der Haupt- 
torsion. Wie A. Kimmel? zeigt, gilt 


¢,—2 
(CS Tane 


wobei der Eigenwert ¢ mit unserem Eigenwert x zusammenhangt durch 
€ =2(1—cosx). 


* Da die Momente und Drehwinkel wegen der Homogenitat des dafiir zustandigen linearen Gleichungs- 
systems nur bis auf einen konstanten Faktor bestimmt werden kénnen, geben wir immer die Momenten- bzw. 
Drehwinkelverhdltnisse an. 

1 C. B. Biezeno u. R. Grammel, Technische Dynamik Bd. 2, S. 413. 

2 4A, Kimmel, Ing.-Arch. 10 (1939) S. 207, 
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On ce ER ell ne tee 
Beachtet man diesen Zusammenhang, entwickelt in (21) den Ausdruck (1 + 2 —5 Cos = in eine 
Reihe und bricht diese nach dem zweiten Glied ab, so erhalt man vollstandige Ubereinstimmung 


mit der auf ganz anderem Weg abgeleiteten Formel (24). 


4, Die homogene Maschine mit Zusatzdrehmassen. a) Die homogene Maschine mit einer 
Zusatzdrehmasse. Die Zahl der Drehmassen im homogenen Teil der Maschine sei* n, die Zu- 
satzdrehmasse sei mit 9,1) bezeichnet (Abb. 2). Die EinfluBzahl der Haupttorsion im hinzu- 
kommenden Wellenabschnitt werde mit 6,,, = 0, bezeichnet, die EinfluBzahl der entsprechenden 
Nebentorsion mit bn,,-1 = On—3,n = On 


"| M, | M, Mn-1| My fo 
0 0 
One1 


Abb. 2. Homogene Maschine mit einer Zusatzdrehmasse. 


Wie man durch Anschreiben von (6) fiir k = 1 bis n feststellen kann, gilt die Rekursionsformel 
(12) nur fiir k=1 bis k= n—2, sofern man 6,) = 6* setat, was wegen der Randbedingung 
M, = 0 erlaubt ist. Die Gleichung (13) ist daher nur fiir k = 0 bis k = n — 1 giiltig. Gleichung (6) 
lautet fiir k =n—lundk=n 


G =o wo 6") M,,_2 ar on Oe 5] J = (5 =F wo? or) M, a 0 > (25) 


1 1 


1 
(5 +0°8:|Mi—(5+% - 


ntl 


— wo? 6n} M, + (e + ot Snmet) Muss = 0%. > (26) 


n+1 
Aus der Randbedingung M, = 0 (links freies Ende) folgt mit (13) 
y = 0. 


Die zweite Randbedingung verlangt M,,,1 = 0 (rechts freies Ende). Beriicksichtigt man dieses 
in (26) und eliminiert M,, in (25) mittels (26), so erhalt man, wenn man auBerdem w nach (15) 
einfiihrt, folgende Beziehung zwischen den Momenten M,,_» und M,,_1: 


a,(*) M,_2 + (1 — 2 a,(x) cos x) M,_1 = 0. (27) 
Dabei bedeutet 
é* * 
(2 + 2 a } He Ue SRG ie | 3] cos] 
Lom n+l n+1 6 6 
a,,(x) O* o* o* 12 2 (28) 
i +2 3 | 2/ 3 3 cos Xx 


Wir fahren nun in (27) M,_; und M,,- 2 nach (13) ein und erhalten nach kurzer Umformung 
die Bestimmungsgleichung fiir die Eigenwerte x in der Form 


a,(x) sinn x — sin (n —1)x =0. (29) 


Von den Lésungen dieser Gleichung kommen nur Werte 0 < x < min Frage. Die Lésungen x = ma 
(m =0,+1,+ 2...) sind trivial, da sie keine Winkelausschlage und somit keine Schwingung 
liefern. AuBerhalb des angegebenen Bereiches liegende Kigenwerte fiihren wegen der Periodizitat 
der Lésungen der Gleichung (29) zu den schon vorher erhaltenen Eigenfrequenzen. 


* Wir bezeichnen die Zahl samtlicher Drehmassen einer Maschine i i A i 
: e immer mit NV. hr 
im homogenen Teil gelegenen Drehmassen bedeutet. SE Teak ee ag 
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Sind die Eigenwerte x reell, so bekommt man die Kigenkreisfrequenzen sowie die Momente M, 


bis M,__; und die Drehwinkel ¢, bis yn; nach (15), (17) und (18) (fiir y = 0): 


2 (1 — cos x) 
oO => or hae SE; ; 30 
61 + ge cosa] ee 


M,, sink x 
MM, sin (k=1,2...n—1), (31) 
cos (i — x x 
Vr (ha1,2%ein= 1), (32) 
Yi x 
cos — 
2 
M. 
Seo (33) 
Py 


Die letzte der angeschriebenen Formeln gibt die Beziehung zwischen dem Moment M, und dem 


Drehwinkel y;. Das Moment M,, muB aus (25) oder (26) bestimmt werden. Nach (26) ist in Ver- 
bindung mit (13) 


Me 1+ w? d*O sin (n—1)x 
Vie (2) sin x : Seed 
1+ — 06,9 
Or «3 


Die Drehwinkel ¢,, und ¢,,,; bekommt man am einfachsten aus (2) mit k =n und k=n + } 
bei Beachtung von (33) 


Pn s; M,, sax Mra 
1 a, M, M, i) (35) 

Care. Ca Ms 
Fi “3 One M, ; co) 


Man kann zeigen, daB der fiir die héchste Frequenz maBgebende Eigenwert unter Umstanden 
komplex wird!. In diesem Fall gilt 


x=na+if, (37) 
und die Bestimmungsgleichung fiir # lautet 
a,(B) Sinn 6 + Stn (n— 1) PB =0, (29’) 
wobei ; 
ee eer en eerie 
a,(B) pee n+1 n+1 (28’) 


oe /6* Of 2 
+2 $—2(5 —F) est o| 


Die Formeln (30), (31), (32) und (34) gehen dann iiber in 


ee 2 (1 + Goj B) ; (30’) 

so(1—25 oof) 
My 7 yess Sinks pee soak 31" 
ar = ( 1) Sint (peal Dan — 1) ; (31’) 

sin(k—4) 
Pk ( 1)F+1 (bs, 26 -5.0— 1), (32’) 
ve ein t 
M 1+026O Gin(n—1)B 34? 
eee tye tl ; (34’) 
M, ‘igi = 86 6 Sin B 
n+1 


wahrend die Formeln (35) und (36) erhalten bleiben. 


1 Vel. H. Bufler und F. Kiefling, Ing.-Arch. 29 (1960) S. 250. 
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Aus der Formel (30') erkennt man, daB die Eigenfrequenz imaginar wird, falls 
rs) . 


Das bedeutet, daB dann eine Frequenz weniger auftritt als bei der Torsion erster Art, bei der es 
bekanntlich fiir insgesamt N Drehmassen N — 1 Eigenfrequenzen gibt. 


b) Die homogene Maschine mit zwei Zusatzdrehmassen auf einer Seite. Wir 
legen fiir den homogenen Teil der Maschine wieder n Drehmassen zugrunde; die Zusatzdrehmassen 
seien O,,,; und @,2 (Abb. 3). Die EinfluBzahlen der Haupttorsionen in den beiden hinzukommen- 
den Wellenabschnitten sind 6,,, = 6, und On+41,n41 = 6n41, die EinfluBzahlen der entsprechenden 
Nebentorsionen 6n,n—1 = On—1,n = O* und 6n41,n = On,n41 = On41- 


MyM 


(abiees ahaa 


Ory 


G, +2 


Abb. 3. Homogene Maschine mit zwei Zusatzdrehmassen auf einer Seite. 


Schreibt man (6) fir k = 1 bis k =n + 1 an, so sieht man, da die Differenzengleichung (12) 
nur fir k —1 bis k =n—2 Giiltigkeit hat, falls man 6,,) = 6* setzt, was wegen der Rand- 
bedingung M, = 0 zulassig ist. Gleichung (13) ist daher nur fiir die Momente M, bis M,_; zu- 
standig. Fir k =n—1, k=n, k =n-+ 1 geht (6) iiber in 


1 , 
(6+ *6*)M,1—(g— 08) Mit (e+e 2)M,=0, (39) 


1 


nt+1 


= on) M,, + (3, + eo? os] Maia = 0, 0) 
nt+1 


1 = 1 
(6 anes 08) Mn ae +% 
1 : 1 
a Os ot1] M.—(o- 


nt+1 Orv. 


1 
=o. dns1}Mngs cA + oe énttinta) Mrny2=9. (41) 


Die Randbedingung am linken Ende M, = 0 ergibt mit (13) 
y=0. 
Die Randbedingung am rechten Ende lautet 
Mry2=90. 


Damit bekommt man aus (39), (40) und (41) nach Elimination von M,.; und M, sowie Ein- 
fiihrung von w nach (15) folgende Gleichung zwischen den Momenten M,_; und M,,_2: 


aa ; 
= — G,(%) On 4a(*) |+ Mn—a 
1423 +2(G—F) cose a 
1435 \? 
— 4a, + 1(x) + 2 ‘a Res — a,(x) a, 41(x) | cos x} M,_; = 0. (42) 
1+2 aa 3.) cose 


Hierbei bedeutet a,(x) die Abkiirzung (28), und es ist ferner 


6*\( O (2) Om 0 © -\ 6* 
(2 eae Ga 5 t+ See oe ee “y cosx} 


OO Okay 2 : . (43) 
(F053) om 


On +4(%) = 
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In (42) fiihren wir wieder M,_; und M,_» nach (13) ein und erhalten dann nach geeigneter 
Zusammenfassung die gesuchte Bestimmungsgleichung! fiir die Eigenwerte x: 


o* 2 
1+ 2 aS 
l O* ‘ OF O* ae a,,(x) Gy + 1(x) sin n x ee An + 1(X) sin (n = 1) eis (44) 
+ 2 +3 +) cos 
6 to) (6) 


Fir reelle Eigenwerte x sind wegen y = 0 die Kigenkreisfrequenzen sowie die Momente M, 
bis M,,_, und die Drehwinkel ¢, bis y,_; aus den Formeln (30), (31) und (32) zu entnehmen. Die 
Momente M,, und M,,,, sind aus (39) und (40) oder (41) zu berechnen. Wir erhalten in Verbindung 
mit (13) aus (39) 


M,, 2—o?d@0O sin(n—1)x 1+ w? 6* @ sin (n — 2) x 


M, 1+0760 sinx  1+o?0*6  sinx Cy, 
und aus (40) 
1 — wd, 0 
Mire Ee Ons+1 me eh 1+ @25*@ — sin(n—1)x (46) 
Mo cee M, © as ainix. =" 
Onn = OO" ao Gis + ow? 6*,.,0 


Gleichung (2), angeschrieben fiir k = n, k =n+1 und k=n-Q, liefert unter Beriicksichtigung 
von (33) die Drehwinkel ¢,, ¢,41 und p49! 


Ga Mh, 


iinet AT 
Pi M, Vie ( ) 
Pn+1 3 0 My +1 de M, (48) 
Zt as On M, M,)’ 
nse =. D-Mass (49) 
Got Ogg Mie 


Es ist nun méglich, daB fiir gewisse Daten der zur héchsten Frequenz gehérige Eigenwert % 
oder die beiden zu den héchsten Frequenzen gehérigen Eigenwerte x komplex werden”. In diesem 
Fall gilt (37), wobei 6 aus folgender Gleichung bestimmt werden muB: 

2 


6* 
jeer el 
: — a,(B) a; +1(8) | SinnB—a, + 1(8) Sin n—1)B=0. (44) 


je. /o* OR 
14 23—2(5—F) oe 


Hierin hat a,(f) die Bedeutung (28’), und es ist weiterhin 


CN (EC a8 Sut avian 5 tg | Got8| 
(+25) [5 4 ys alg Boa 5 + BE] Cof 6 


i * * * 2 
C) ree eal O60 1) cof 


(43') 


On + (8) = 


OE 


Fiir die Eigenfrequenzen sowie die Momente M, bis M,_, und die Drehwinkel Pi bis Paci gelten 
die Formeln (30), (31’) und (32’). Setzt man in (45) und (46) « =a +78, so erhalt man fiir die 
Momente M,, und M,,+, 
2—w?dO Gin (n—1)8B 1+ @?6*O Gin (n— 2) 


a : asl 45" 
a leaped, jem Capo | lt oO Ging | ( 
) 
Wiener eo ye Lteree Sings (46°) 
Mig a 7 ye. 6 Gimp 
ul oe ete OO a + 81 O 
n+1 n+1 


Di i i berechnen. 
Die Drehwinkel und %,+, sind aus (47), (48) und (49) zu ; 
3 Ist die Bedingune (38) fiir einen oder zwei Eigenwerte erfillt, so bedeutet das, daB eine bzw. 
- zwei Frequenzen weniger auftreten, als es die Torsion erster Art verlangen wiirde. 
1 Sinngem4B gilt auch hier die im AnschluB an (29) gemachte Bemerkung. 
2 Siehe FuBnote 1 von Seite 379, 
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c) Die homogene Maschine mit je einer Zusatzdrehmasse auf beiden Seiten. 
Der homogene Teil der Maschine umfasse wieder n Massen; die polaren Massentragheitsmomente 
der linken bzw. rechten Zusatzdrehmasse seien 0, baw. 0,11 (s. Abb. 4). Die EinfluBzahlen der 
Haupttorsionen in den beiden hinzukommenden Wellenabschnitten sind 09,9 = J) und 6,,, = On 
die EinfluBzahlen der entsprechenden Nebentorsionen 6, 9 = 69,,; = 0) und 0n,.—1 = Only 0 


j f M, Mp4 My, My +1 
0 0) 0 0 
Q, On 44 


Abb. 4, Homogene Maschine mit je einer Zusatzdrehmasse auf beiden Seiten. 


Aus (6), angeschrieben fiir k = 0 bis k = n, erkennt man, dafB (12) nur fiir k = 2 bis k=n—2 
gilt. Infolgedessen erhalt man die Momente M, bis M,_; aus (13). Mit k = 0, k = 1 baw. k = 
n—1, k = n folgen aus (6) die beiden Gleichungspaare: 


Linkes Ende: 


1 1 u ace ‘A * 
(3, + w? 40,3 M_, —l6 + eee 0] M, + (o + @? 66) wy 0, (50) 
(= + w? ot) My) — ce ow 6) M,+ c + w? oy M,=0, (51) 
Rechtes Ende: 
e +o 0*) M2 ale wo 6) May (e + or] M,=0, (52) 
il 
(6 aa wo” 62) ME Sots a (5 i (6) — 2. wo 6,] M, ate SQ i wo Santa} My +1 ae 0 ss (53) 


Die Randbedingungen lauten hier 
M_, — 0 und WE Say => 0 . 


Beriicksichtigt man diese in (50) und (53) und eliminiert M, aus (50) und (51) sowie M,, aus (52) 
und (53), so bekommt man, wenn man auBerdem noch w nach (15) einfihrt, 


[1 — 2 a,(x) cos x] M, + a(x) M, = 0, (54) 
a,(x) M,_. + [1 — 2 a,(x) cos x] M,_, = 0, (55) 
wobei a,(x) die Abkiirzung (28) und 


On 0% % @\d* 6, | 
Grg)it gee eali+ g) 5+ leo 


ox o* oF 2 
1 i) 0 
#23 42(5 3 oes 


ay(x) = 


(6) IF 


bedeutet. In (54) und (55) setzen wir fiir die Momente M,, M,, M,—2 und M,_, die Ausdriicke (13) 


ein; sodann folgt, wenn man nach sin y und cosy ordnet, 
{ao(x) cos 2 x + [1 — 2 ag(x) cos x] cos x} sin y 
+ {ao(x) sin 2 x + [1 — 2 a,(x) cos x] sin x} cosy = 0, (57) 
{@n(x) cos (n — 2) x + [1 — 2 a,(x) cos x] cos (n — 1) x} sin y 
+ {a,(x) sin (n — 2) x + [1 — 2 a,(x) cos x] sin (n — 1) x} cosy = 0. - (58) 


Die Determinante der Koeffizienten von sin y und cosy muf verschwinden, woraus sich nach 
geeigneter Zusammenfassung folgende Bestimmungsgleichung fiir die Eigenwerte x ergibt*: 


a,(x) a,,(x) sin n x — [ao(x) + a,(x)] sin (n — 1) x + sin (n— 2) x =0. (59) 
* Sinngem48 gilt auch hier die im Anschluf an (29) gemachte Bemerkung. | 
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Fir reelle Eigenwerte x miissen die Eigenkreisfrequenz sowie die Momente M, bis M,_ und 


die Drehwinkel ¢, bis y,_; den Formeln © 


2 (1 — cos x) 

oO= i ame (60) 
§0(1 +255 cosa} 
M, sinkxctgy + coskx : : 

MM, sinxcetgy+cosx ’ (A =1,2...n—]), (61) 

cos (k 7) vty —sin(k 5)» 
= 4 5 : (c= tae — 1). (62) 

= cos -> x ctg y —sin >-x 


die sofort aus (15), (17) und (16) folgen, entnommen werden. 

Im Gegensatz zu friiher haben wir in (62) den Bezugswinkel ~, eingefiihrt, da hier fiir :, Glei- 
chung (16) nicht mehr zustandig ist, wie man aus (2) mit k = 1 erkennen kann. Fiir die weitere 
Rechnung bendtigen wir noch das Verhaltnis 


M, w?@ 


sin x ctg y + cos x 


Pathe 2 sin Ghee ee we 
2 2 2 2 
Der Ausdruck ctg y kann entweder aus (57) oder (58) entnommen werden; nach (57) ist 
a ao(x) cos 2 x + [1 — 2 a)(x) cos x] cos x 
eh dea a)(x) sin 2 x + [1 — 2 a(x) cos x] sin x 
(eee) iL eas wie 2|(1 ee coss} 
6 aes @,) 8 ' 6 = 
a ox o* Of eevs — cig x. (64) 
es ae 
f ! 2 +2(5 3) esa sin x 
Die Momente M, und M,, berechnet man zweckmafig aus (50) und (53): 
M, 1 + w 6¢ O 
Me ae ee 
(Oy v 
M, 1 + w? 6% O sin (n — 1) x ctg y + cos (n — 1) x (66) 
Mi — 08,0 sin x ctg y + cos x 
n+1 
Die Drehwinkel @, 1, Pn» Pn 41 folgen aus (2) mit k = 0, 1, n, n — 1 bei Beriicksichtigung von (63): 
ee] @ M, 67 
Po — Wa) Fe (67 
Gi _ Ba 68 
2 = be) (1— 3p) (68) 
Tire, Men Mn 69 
. p= be) (ae me)? Oe 
0 M, 
Ae a tel 70 
Pp (*) On+1 Mi’ ( ) 
wobei wir zur Abkiirzung ip 
or% sin x ctg y + cos x 7] 
i) =< oe cy 
2(cos 7 x egy —sin 5-4] sin > 


einfihren. ae 
Man kann zeigen, daB fiir gewisse Daten der zur héchsten Frequenz gehérige Kigenwert x 


oder die beiden zu den héchsten Frequenzen gehérigen Eigenwerte x komplex werdent. Dann 


gilt (37), und f ist aus folgender Gleichung zu ermitteln: 


—— 


ao(B) a,(8) Sin n B + [ao(B) + an(B)] Sin (mn — 1) B + Sin (n — 2) eae (59") 


1 Siehe FuBnote 1 von Seite 379. 
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Aus (28’) ist a,(8) zu entnehmen, und aj() ist gegeben durch 
5 ) 
6*\{,  O bo ( es ac 
(+ 2%\h+€ 2< 2+ ola t 3 of B 
38 4 (* OF 2 
i +2 —2(5 = 3) ooi8) 


Fiir die Eigenfrequenzen sowie die Momente M, bis M,_1 und die Drehwinkel gy, bis Yn_1 gelten 
hier die Formeln 


(56’) 


ao(B) = 


a (60') 
s0(1—25 CofA) 
M, : Sink B Ctg y —Cojf kB (1 20) eee 61’) 
Mm, ~~" ~ Sing Gig — Gor ( : ers ( 
1 x 1 
Sin(k—1 ] a6ig7 —Goi(k—+)p 
Pe —_ (1) : —. (p= 3.308. . (en) 
Pa Sin-> B Gta y —Cof > f 
Der Wert Stq y ist dabei gegeben durch 
sir g-rbali eed 
P24 ie 2° — 2](1+—])— + | Gof 6 
ean ee a ey (64') 
| eet) f i 
p+ 2F—a(5 3) Soi] eine 


Das Moment M, ist durch (65) bestimmt, fiir M, ergibt sich 
M, ah (es Lye ( a w? dx) [Sin (n — 1) B Gtg y — Gof (n — 1) 8] ; (66’) 


os tere — 08 5, 6) (Einf Cig 7 —Col A) 
nt+1 7 
SchlieBlich erhalten wir fiir die Drehwinkel 
Vous b oO M 67’ 
fo — w(B) Oe, (67) 
Pi ey i 1 _ M 63’ 
B= H(B)(1—e) (68') 
Pn Sey Ali ‘M, mee | 69! 
pe M, M, |, ( ) 
Pn+1 wd Me C) M, 70' 
Po a og On+1 M,’ ( ) 
wobei 
Sin B Ctg »y —Cof B 
(8) = —— oy (71) 
2(Gin 5 Ein 7 — Co} > A) cof f 
ist. 


Falls die Bedingung (38) fiir einen oder zwei Eigenwerte zutrifft, so gibt es eine bzw. zwei 
Eigenfrequenzen weniger als bei alleiniger Beriicksichtigung der Torsion erster Art. 


5. Beispiel. Als Beispiel betrachten wir einen Vierzylindermotor mit Schwungrad (Abb. 5). - 
Die fiir die einzelnen Kolben maBgebenden Tragheitsmomente seien 0, = 0, = 0, = O, = O, 
die EinfluBzahlen der Haupttorsion zwischen den Zylindern 6,1 = 69,5 = 63,3 = 6 und die ent- 
sprechenden EinfluBzahlen der (primaren) Nebentorsion 01,2 = 02,3 = 0* = 0,15 6. Das Schwung- 
rad besitze das Massentragheitsmoment @, = 10 @, das Kurbelwellenstiick zwischen Schwungrad 


und viertem Zylinder habe die EinfluBzahl der Haupttorsion One Of = 2 6; die EinfluBzahl der 
(primaren) Nebentorsion 6, , = Of sei 6f = 0,10 6. 


et ec0T0 | a Fe ee = 000°0 — | or — = — pate Pes 
szr‘o— | zsor‘o ZS0T‘0 0000°T OLZTT PLS‘0— | zS0T— OT att cero 0 ecco | oo | + 
LESTI— | 8zr‘o— rer0— OLZTT SO6TT Le60— | pz0'0— T I 0S2°T oro oot | sto | ¢ 
ieeo— |. rest Le T— SO6T‘T SO6T‘T Z10°T £19°0 I I 008° STO oor | sto | z 
000° Lee‘o— 09%°0— SO6T‘T 0000°T 09" 000° I I OST St‘0 0 | 1 
ig 13 I H = 
— 2 & Or 2 . G*[°,0 + I “Be 
& é 


apanm 1z4asad T = ~ pun [ = Q 2agom *7ZT‘T =*@ =o pun sjaidsiag uajaiysv.yaq sap uainq uep nu | ayaqny, sep Sunuysasyoingy *¢ aT[oqey, 


(2H u1 zuenbosfuaty +f) jaidsiag araryon.yaq sop anf assruIPYysayayuimyaig pun -uaquawopy ‘uazuanhosfuasig +z e119qe yy, 


Drehschwingungsberechnung von Kolbenmaschinen 
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Mit den angegebenen Daten erhalt man aus (29) mit (28) folgende Gleichung fiir die Eigen- 


werte x: 
0,564 + 1,295 cosx . 2 * 

ue 2 ? f = = 0 : 1 
f@~= (1,2 + 0,1 cos xP sin 4% — sin 3 x (1*) 


Die graphische Lésung dieser Gleichung zeigt Abb. 6 und liefert in dem in Betracht kommenden 
Bereich 0 < x< az die drei Kigenwerte 
@,=0ATZ > 4, = 1,261 , . x, = 2,083. 


F(a) 


Abb. 5, Vierzylindermotor mit Schwungrad. Abb. 6, Graphische Lésung der Gleichung (1*). 


47°(p) 


B5- 


-95 


Abb. 7. Graphische Lésung der Gleichung (2*), 


Abb. 9. Restwertkurve fiir das betrachtete Beispiel. Abb. 8. Graphische Darstellung der Eigenschwingungsformen und 
der in den einzelnen Wellenabschnitten auftretenden Momente 
fiir » = 1 (Grundfrequenz) bis »y = 4 (hichste Frequenz). 
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Der vierte Eigenwert muB aus (29’) mit (28’) errechnet werden: 


cr; pp, —. 9,964 — 1,295 Gof B AW. a, 
f (3)== (2-01 Ea h Sm 4B + in 3p == (I, (2*) 


Die Kurve f’ = f’(8) ist in Abb. 7 dargestellt; man erhalt daraus 
B = 0,328 . 


Nach Einsetzen der so gefundenen Eigenwerte in (30) bzw. (30') ergeben sich die in Tabelle 2 
eingetragenen EKigenkreisfrequenzen. Die Momenten- und Drehwinkelverhaltnisse folgen alsdann 
aus den Gleichungen (31), (32), (34) baw. (31’), (32'), (34’) sowie (35) und (36). Sie sind in Tabelle 2 
zusammengestellt und in Abb. 8 graphisch aufgetragen. 

Zum Vergleich lésen wir die vorliegende Aufgabe noch nach dem erwahnten Tabellenverfahren. 
Zu diesem Zweck errechnen wir fiir eine Reihe von w-Werten den Restwert R (In Tabelle 3 ist die 
zahlenmaBige Rechnung fiir w = w, durchgefiihrt); die sich ergebende Kurve R = R(w) zeigt 
Abb. 9. Aus den Schnittpunkten mit der Abszisse folgen sofort die vier Eigenkreisfrequenzen 
© 1, Wg, W3, Ws, welche mit jenen in Tabelle 2 angegebenen iibereinstimmen. Mit diesen Werten 


liefert Spalte 7 der Tabelle 1 die in den einzelnen Wellenabschnitten iibertragenen Momente und 
Spalte 12 die Drehwinkelverhiltnisse. 


6. Zusammenfassung. Einer genauen Drehschwingungsberechnung von Kolbenmaschinen 
ist die Torsion zweiter Art zugrunde zu legen. Fiir je drei aufeinanderfolgende Drehmomente 
14Bt sich eine Rekursionsformel angeben, welche sich von derjenigen bei der Torsion erster Art 
dadurch unterscheidet, daB alle drei Koeffizienten der Momente ein mit der Higenkreisfrequenz 
behaftetes Glied enthalten. Im Fall einer homogenen Maschine vereinfacht sich diese Rekursions- 
formel auf eine Differenzengleichung mit konstanten Koeffizienten; ihre Lésung fihrt auf einfache 
trigonometrische Ausdriicke. Insbesondere lassen sich im rein homogenen Fall die Kigenfrequenzen 
streng und in geschlossener Form angeben. Ferner werden die praktisch wichtigen Falle einer 
homogenen Maschine mit einer oder zwei Zusatzdrehmassen eingehend behandelt; allerdings 
miissen dabei die Kigenwerte aus im allgemeinen nicht explizit lésbaren transzendenten Gleichungen 
bestimmt werden, wobei zu beachten ist, daB auch komplexe Eigenwerte méglich sind. Es zeigt 
sich weiterhin, daB die einem komplexen Eigenwert entsprechende Frequenz (theoretisch) imaginar 
werden kann; dieser Fall diirfte jedoch in der Praxis kaum eintreten, da, wie man durch eine 
genauere Untersuchung erkennt, hierzu die EinfluBzahlen der Nebentorsion die GréBenordnung 
derjenigen der Haupttorsion erreichen miBten. Fiir die behandelten Schwingungssysteme werden 
auBerdem Formeln zur Ermittlung der Momenten- und Drehwinkelverhdltnisse abgeleitet. Als 
Beispiel wird abschlieBend eine Vierzylindermaschine mit Schwungrad vollstandig durchgerechnet. 


(Eingegangen am 14. Dezember 1959.) 


Anschrift der Verfasser: Dr. H. Bufler und Dipl.-Ing. F. KieSling, Miinchen 2, Arcisstrafe 21, Institut fiir 
Technische Mechanik der Technischen Hochschule. 
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Geniherte Berechnung von Warmespannungen 
mit Hilfe der Variationsprinzipien der Elastostatik * 


Von R. Trostei 


1. Einleitung. Der iibliche Weg der Berechnung von Warmespannungen aus den thermisch- 
elastischen Grundgleichungen fiihrt schon in einfachsten Fallen auf auBerordentlich umfangreiche 
Rechnungen. Oft ist sogar eine exakte Lésung iiberhaupt nicht méglich. Es liegt daher nahe, 
z.B. auf der Basis des Ritzschen Verfahrens Naherungslésungen fiir Warmespannungsprobleme 
zu suchen. Hierfiir ist es erforderlich, das Variationsproblem der (als Eulersche Differential- 
gleichungen anzusehenden) thermisch-elastischen Grundgleichungen zu kennen. Die Formulierung 
dieses Variationsproblems ist mit Hilfe des Hamiltonschen Prinzips! bzw. im statischen Falle mit 
Hilfe des Prinzips der virtuellen Verriickungen grundsatzlich méglich. Allerdings sind in den 
hiermit aufzufindenden Variationsintegralen lediglich die Verschiebungen als zu variierende 
GréBen enthalten, so daB sie sich im Hinblick auf eine Naherungsberechnung von Warmespannungen 
in vielen Fallen als ungeeignet erweisen. Nun laBt sich jedoch dem Prinzip der virtuellen Ver- 
riickungen ein Variationsprinzip gegeniiberstellen?, das aus jenem mittels einer Variablen-Trans- 
formation hervorgeht, und in dem nunmehr die Komponenten des Spannungstensors als zu vari- 
ierende Funktionen erscheinen. Es entsteht ein Variationsprinzip, das auch unmittelbar aus dem 
Prinzip der virtuellen Krafte hergeleitet werden kann. Der maBgebende Energieausdruck dieses 
Variationsproblems findet sich bereits bei Miiller- Breslau®, der in einer im Jahre 1884 erschienenen 
Arbeit mit Hilfe der Arbeitsgleichung der Statik das Castiglianosche Prinzip auf thermoelastische 
Vorgange erweiterte. Eine allgemeine Fundierung der Arbeit Miiller-Breslaus unter Einbeziehung 
warmetheoretischer Uberlegungen wurde von Domke! im Jahre 1924 geliefert, wahrend in neuerer 
Zeit zu diesem Problemkreis von Parkus! durch die Formulierung des Hamiltonschen Prinzips 
fiir thermoelastische Vorgange ein vertiefender Beitrag geleistet wurde. 

Eine praktische Auswertung thermoelastischer Variationsaussagen im Sinne des Ritzschen Ver- 
fahrens scheint bisher noch nicht in Angriff genommen worden zu sein. Diesem Umstand soll 
vorliegender Aufsatz* Rechnung tragen, in dem an einigen Beispielen die Leistungsfahigkeit einer 
Naherungsrechnung bereits mit einfachsten Spannungsansatzen demonstriert werden soll. Sogar 
die Temperaturfelder selbst werden hierbei durch einfache Ansatze approximiert. Einem Beispiel 
fiir stationare Warmespannungen folgt die Untersuchung eines instationdren Problems (Ab- 
kihlung einer Kreisscheibe), wobei jeweils den exakten Lésungen die mittels der Naherungs- 
ansatze gefundenen Beziehungen zum Vergleich gegeniibergestellt werden. AbschlieBend wird 
die Naherungslésung eines raumlichen Problems an einem Hohlzylinder angegeben. Dabei wird 
in Erweiterung des Ritzschen Gedankens eine Methode zur genaherten Lésung von Variations- 
problemen benutzt, bei der man auf einfach lésbare Differentialgleichungen fiir die Naherungs- 
funktionen gefiihrt wird.> Den Beispielen vorangestellt ist ein kurzer Abri® der theoretischen 
Grundlagen. 


2. Theoretische Grundlagen. a) Allgemeine thermoelastische Vorgange. Bezeichnet 
man mit [ eingepragte Krafte je Volumeneinheit und mit 8, Spannungen an der Kérperoberflache, 
so wird langs einer Variation des durch den Verschiebungsvektor 1 charakterisierten Deformations- 
zustandes an einem Kérper vom Volumen V und der Oberflache O(V) die auBere Arbeit 


OuA, = ftdudV+ f3,dud0 (1) 
[Vv] [O(V)] 


* Gekiirzte Fassung einer von der Technischen Universitat Berlin, Fakultat fiir Bauingenieurwesen, ge- 
nehmigten Habilitationsarbeit. Der Verfasser spricht hiermit den Berichtern, Herrn Prof. I. Szabé und Herrn 
one ie Teichmann seinen Dank aus. : 

. Parkus, Uber eine Erweiterung des Hamiltonschen Prinzips auf thermoelastisch ange. - 
hofer-Girkmann-Festschrift, Wien 1950" , oe eee 

2 R. Courant u. D. Hilbert, Methoden der mathematischen Physik, 2. Aufl., Bd. 1. 

3H. Miiller Breslau, Z, Arch. u. Ing.-Ver., Hannover 30 (1884) S. 211. 

* 0. Domke, Eisenbau, 12, Nr. 5/6. 

5 Dies verdanke ich einer Anregung, die ich anlaBlich eines Seminars am Lehrstuhl fiir Mechanik (Prof. 
Szabé) von Herrn Prof. Morgenstern und Herrn Dr. André empfing. 
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geleistet!. Beachtet man hierin die Oberflachenbedingung 


3, = nS, 
wobei © der (mit Riicksicht auf das Boltzmannsche Axiom) symmetrische Spannungstensor und tt 


der Vektor der auSeren Oberflachennormalen ist, so erhalt man unter Heranziehung des Gauf- 
schen Satzes fiir das zweite Integral 2 


,dudO= f(n@)dudO= / (5uS)ndO = (V (GowaV— sive a 
abit iO f 2) 0AO=, | eyn40 = 0 one — [Fe bars Johar 


und damit fiir die Variation der a4uBeren Arbeit 
Oy A, a eee ee eV) Su dV. (2) 
Nun ergibt unter Heranziehung des Newtonschen Grundgesetzes fiir das Kérperelement 
VS+t=0% 
eine Umformung des ersten Integrales in (2) 


[(vS+) oar = fegerar=s| [ eG ouav 
[V] (V] [V] 


ra) ou 
eg Cz C= 
[(V] [Vv] 


Q au 0 /eu\? 8 
a| fea =f Sa) v= 
[Vv] 


[Vv] 


ou a 
=~ | 0 2 (ou) av 
[Vv] 
ou ou 
ou (Fe) aV 


I 


ey NE eK 


7] 
wobei 


ix Q /ou 2 
K=[3(3) 4. (3a) 
[Vv] ; 
die kinetische Energie des Systems ist. Auch das zweite Integral in (2) laBt sich umformen. Hierfiir 
benutzt man die Identitat 
t 
(SV) bu=G..(V edu) =G.. dy (V ot), 


wobei der Doppelpunkt die doppeltskalare Multiplikation® des Spannungstensors © mit der Dyade 
YV cou andeuten soll, und erhalt schlieBlich, da das doppeltskalare Produkt des (symmetrischen) 
Spannungstensors mit der antimetrischen Dyade 


x lV odu—(V o du) 


verschwindet 4, 


ub D a 
(OW ou = ©. 05 S(Vout (Veu)}} =O. D= 6, fS..dD= 4,4. (3b) 
Dabei ist 


D=s(Veut (Vou)! | (30) 


der Verzerrungstensor und A; eine EnergiegréBe, die bei gleichmaBiger Temperaturverteilung 
T = T, in die spezifische Formanderungsenergie tibergeht. Somit ist also 


t a 
f(GV) dudV = [S..6,DdV = 6, f ApdV = 0, 4,, (3d) 
[V] [V] (VY) 


so da8 man schlieBlich aus (2) 


: é ou 
 : <a aude ca 
[VY] 


1 Der, dem Variationssymbol 6 beigefiigte Index u soll darauf hinweisen, da8 die in (1) dargestellte Arbeits- 
variation allein durch eine Variation des Verschiebungsvektors u verursacht zu denken ist. : : 
2 Durch den Pfeil soll zum Ausdruck kommen, das der Nabla-Operator allein auf du anzuwenden ist. 

8 Siehe hierzu M. Lagally, Vektorrechnung, 3. Aufl., Leipzig 1945. 

4 (VY o du)! ist die zur Dyade V7 0 du transponierte Dyade. 


27 
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und nach Integration hinsichtlich der Zeit mit der MaBgabe du(t)) = du(y)) = 0 
t; = 
dy f [4;—4,— K] dt = 0 (4) 
ty 


erhalt. Hierbei ist fir die Arbeit der auBeren Krafte, wenn wir allgemein die Belastung als von den 
Verschiebungen abhangig annehmen, 
u ut 
=f (ita dV+ f J sau) a0 (5a) 
[Vv] \9 [O(Y)] \o 
und speziell bei konstanter Belastung 


A,=ftudV+ f8,udO (5b) 
; [OV )] 


[Vv] 
zu setzen. 


Fiir die Darstellung von A; baw. 
A, =f. .dD (6) 
0 
in Abhangigkeit von den Verformungen benétigt man ein Spannungs-Verformungsgesetz, das 


fiir elastische Kérper und kleine Deformationen geniigend genau mittels des durch die Warme- 
ausdehnungsglieder verallgemeinerten Hookeschen Gesetzes 


6 =755[9 + (ay V9 pam (P— TG (7) 


l+y 
beschrieben werden kann. (E ist der Elastizitatsmodul, » die Querkontraktionszahl, «, der lineare 
Warmeausdehnungskoeffizient, T — T, die Temperatur gegeniiber einer iiberall gleichen Anfangs- 
temperatur JT, des spannungslosen Ausgangszustandes, (© der Kinheitstensor, \/ u = Sp D 
=..D =e, +e, +e, =e die erste Invariante des Verzerrungstensors). Mit Beachtung von 


D..dD = + d(D..D), 
(Vu) G..dDd =(VudC.9 =Vuadyu =F ayuy, 
(T—T.) 6... dd = (T— T,) d(G..9) =(T—T,) (Vu) 
erhalt man damit aus (6) schlieBlich 


ae E Vv lt+y 


so daf das Variationsprinzip fiir thermoelastische Vorgange letztlich 
ty 


E 1 
ast) fran? + aT 2p TV 
t (UY) 
— 20%) (edu @ "(By )a7— J (ste) ao} at =o (9) 
d [o(v)] \o 


lautet. Der einem bestimmten Temperaturfeld ¢ = T— T, und einer bestimmten Belastung 
zugeordnete Verschiebungszustand macht das Integral J extremal. Diese Aussage ist, wie sich 
durch Nullsetzen der ersten Variationsableitung bei Beachtung der Oberflachenbedingungen 
sowie mit OU(to) = dUu(t,) = 0 zeigen 1aBt, der Verschiebungsgleichung 


1 \ 2 
Au + a5 VV) — So, V (TT) +7 Es 2, ee (10) 


aquivalent. 

b) Bemerkungen zum Temperaturfeld. Da das Variationsproblem (9) keinen Auf- 
schluf tiber das Temperaturfeld gibt, ist noch zusatzlich die Warmeleitungsgleichung heranzuziehen. 
Ihre verallgemeinerte Form!, die bei der aufzustellenden Warmemengenbilanz auch die Abkih- 
lungstendenz eines sich ausdehnenden Volumens beriicksichtigt, lautet 


AA(T—T) +0 =ey— (TT) +y%—* Ay 


3 a 


_ wobei y das spezifische Gewicht, 2 die Warmeleitzahl, c, die spezifische Warmemenge bei kon- 
stantgehaltenem Volumen, c, die spezifische Warmemenge bei konstantgehaltener Normalspan- 


* Siehe F. Neumann, Theorie der Elastizitat, Leipzig 1885 sowie FuBnote 4 von Seite 388. 
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nungssumme o = €..G und Q die je Volumeneinheit erzeugte (eingepragte) Warmemenge ist, 
oder mit Beachtung von! ; 
Ce hay FT 
eu Td =a 
= z. B. ‘os einem KreisprozeB unter Heranziehung des Nernstschen Theorems gefunden werden 
ann, auc 


0 at 
AA(T—T) +0 =47 5 (T— Th) +5q 35 o (VW). (11) 


Hierin stellt die linke Seite der Gleichung den Warmemengenzuwachs je Zeiteinheit im Einheits- 
volumen dar, der durch auBere Warmemengenzufuhr Q und den durch Warmeleitung verursachten 
Warmezustrom bedingt ist. Dieser Warmemengenzuwachs dient einerseits zur Temperaturerhéhung 


der zunichst konstantbleibend zu denkenden Volumeneinheit Cy iS (T — T,)| und andererseits 


zur eet der einmal angenommenen Temperatur T bei seiner tatsachlichen Ausdehnung 
Cp — Cy . . . 
P oo (V w)| - Mit (11) offenbart sich nun die grundsatzliche Schwierigkeit, daB im all- 


oo 3 at 
gemeinen der Verschiebungszustand u und das Temperaturfeld T miteinander verkoppelt sind, 
so daB in diesen Fallen eine Auswertung des Variationsproblems (9) im Sinne des Ritzschen Ver- 
fahrens bei vorheriger gesonderter Berechnung des Temperaturverlaufes unméglich ware. Nur 
im stationaren Fall (0/0t = 0) ist eine vorherige, gesonderte Berechnung der Temperaturverteilung 
strenggenommen moéglich. Nun zeigt sich jedoch, daB insbesondere fiir reine Warmespannungs- 
probleme bei den iiblichen in Frage kommenden festen Kérpern mit guter Naherung eine Trennung 
von Verschiebungs- und Temperaturzustand auch im instationaren Falle vorgenommen werden kann, 
indem man auf der rechten Seite der Gleichung (11) den letzten Ausdruck gegeniiber dem ersten 
Ausdruck vernachlassigt, so daB nur noch eine Gleichung fiir das Temperaturfeld allein verbleibt: 
Cyy @ 7 
AT—T)=S F(r—1) -Paesir_my—2, ea 8%. ay) 
Setzt man namlich in (11) in erster Naherung (Vu % + 3«,(T— T,), was zumindest im all- 
gemeinen der richtigen GréBenordnung der Volumendilatation entsprechen wird, so hat man 
2 
AA(P—T) +O0=4y [Lt poy] ge (T—M=ay ing (T—T)- (13) 
Der Wert x ist in Tabelle 1 fiir eine Temperatur von z. B. 250 °C = 250 + 273 = 523° (abs. 
Temp.) fiir einige Materialien berechnet. Man erkennt, daB, abgesehen von Aluminium und Zink, 
der in (12) durch die Vernachlassigung von x begangene Fehler sich auf die zeitliche Ableitung 
der Temperatur mit weniger als 5° auswirkt. Besonders unerheblich wird die Vernachlassigung 
bei den wichtigsten Materialien Stahl, GuBeisen und Beton. 


Tabelle 1. 


‘ keal 
Material E [kg em-?] ao, [°C] y [kg cem™°] Cp [| y % 
Aluminium 7,2-10° | 2,5-10-5 2,75 + 10-3 0,22 0,3 6,80 + 10-2 
Stahl 21a 08 1,2-10% Bo. LOm 0,11 0,3 3,21 - 10-7 
GuBeisen A 0° 1,0- 10-5 UO) kOEe 0,13 0,3 0,64 - 10-7 
Bronze 90/10 12 - 108 1,8 105% 8,80 - 10-8 0,09 0,3 4.51 - 10-7 
Nickel 20. - 10° 1,3 - 10-6 8,80 - 10-3 sone Os Ll; 0,3 3,20 - 10-7 
Zink 10-105 2,6- 10-5 7 Also 1s 0,09 0.3. 9,65 - 10-7 
Beton 251 10% 1,2). 10> 2520-104 0,21 0,16 0,39 - 10-2 


c) Statische Vorgange. 1) Das Prinzip der virtuellen Verriickungen. Ist die zeitliche Ande- 
‘rung der mechanischen bzw. thermischen Belastungen so geringftigig, da8 man die kinetische 
Energie als von héherer Ordnung klein vernachlassigen kann, so hat man aus (9) 


ty 


| E : 2 (1+) 20+» 5 
i fy fay (®- ®t reas TWP T9 V8 ale 


ty 
u 
= f (fsvat)ao}a1=o, 
forvy1 \b 


1 T ist die absolute Temperatur und J das mechanische Warmedaquivalent. : f 
2 Bei noch héheren Temperaturen sind gréGere Fehlerquellen schon in der im allgemeinen nicht beriick- 


sichtigten Temperaturabhangigkeit der Materialwerte enthalten. 


ma 
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Dies trifft streng nur bei stationaren Problemen zu, wo die GréBen f, 3,, T und damit u zert- 
unabhangig werden, so daB nach Ausfiihrung der Integration tiber t und Division durch (t, —t)) das 


Variationsproblem 


E y 2(1-+») 20+») ry 
iltab [racn(®2 + eal To ay ll —1y) VE E eee dV 
[VY] 


fo ea do =0 (14a) 
-f (jan) 


bzw. mit Einfiihrung der ersten und zweiten Invarianten des Verzerrungstensors’ 
é= VHC... D =e, +e ey: 
= 1 1 
5 = 16... D8 —D..D] = 24 89 +25 + &eu— | Whe + YPy + Yea) 
das Problem 
Ou EL: = 0 P) 


(4b) 


S _ 
sy ase ee 4») (T—Tp) ] dV (40) 


u u 
~ J (fea dV if [f+ a)u0 
[Vv] \o [O(V)] \o 
verbleibt. Wir wollen die Giiltigkeit von (14a) jedoch auch fiir langsam zeitlich veranderliche 
Vorginge voraussetzen, wie sie i. allg. bei instationdren Warmespannungsproblemen vorliegen. 
Die zugehérige Eulersche Differentialgleichung ist mit Gleichung (10) bis auf ihre rechte Seite, 
die im vorliegenden Falle durch Null zu ersetzen ist, identisch. Das Variationsproblem (14a) 1aBt 
sich auch unmittelbar aus dem Prinzip der virtuellen Verriickungen. 

dy (A; — A,) = 0 (15) 
unter Heranziehung des Hookeschen Gesetzes (7) herleiten und geht somit aus einer die Gleich- 
gewichtsbedingungen umfassenden Energieaussage hervor. 


2. Das Prinzip der virtuellen Kréfte. Ein dem Problem (14a) aquivalentes Variationsproblem, in 
dem die Spannungen als zu variierende GréBen auftreten, gewinnt man durch folgende Uberlegung?: 

Ersetzen wir in (14a) samtliche Verschiebungsableitungen nach MaBgabe des Hookeschen 
Gesetzes durch die Spannungen, so erhalt man schlieBlich, nachdem man noch den Ausdruck 


3 
| zza=m [Ea, (T— Ty) P dv 
ta 


gestrichen hat, da er einen, bei vorgegebener Temperaturfunktion festen Wert darstellt, das Pro- 
blem 


a u u 
6,7 =0, f= (re HEB. ftaay— ii [J ax) 0, 
[V] 0 [o(v)] \o 
wobei man sich bei Bildung der Variationsableitung hinsichtlich u den Spannungstensor © 
selbstverstandlich im Sinne von (7) durch die Verschiebungsableitungen ausgedriickt zu denken 
hat. Die hierbei auftretenden Variationen 6, © sind, wie bereits der Index 1 zum Ausdruck 
bringt, nicht willkiirlich, sondern von der Verschiebungsvariation du abhangig. Variiert man in re 
neben U auch nach den Spannungen selbst, so mus daher gleichzeitig (7) als Nebenbedingung 
genannt werden. Man gelangt so zu dem isoperimetrischen Problem 


dy,¢,. 7 =0, 
i= f fp +7) S..S—6..S)'] — f tdu 
[¥] 0 
Ade et Dee Ge Eee ie ee 
Tey 1l—2y~"" 1—2y°! 0) al [8a ujdo, 
(ory) \6 


* Durch eg, eg und e, seien drei zueinander orthogonale Richtungen beschrieben. 
? 'S. hierzu R. Courant u. D. Hilbert, Methoden der mathematischen Physik, 2. Aufl. Bd. 1, S. 201 ff. 


(16) 
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wobei u, © sowie der als Lagrangescher Parameter fungierende Tensor A als zu variierende Gré®en 
anzusehen sind. Die natiirlichen Bedingungen von (16) erhalt man nun durch Nullsetzen der 


Variationsableitungen hinsichtlich A, © und u. Die Variation nach A (34 IT = 0) verifiziert sofort 
die oe ees ), Nullsetzen der Variationsableitung hinsichtlich © liefert 


(7 
[1 +») S..6S—»(..6) G.66)] + A.. 66 lav =0 
Ari | 


bzw. mit 
auch 

He | tel fy) S—v(E..G) EG] +A].. 1S av =0, 
woraus vA 


eld +) S—r»G..S)G+A=0, also A=— "6-75, (17) 


gefolgert werden kann. Demgegeniiber erhalt man durch Nullsetzen der Variationsableitung 
hinsichtlich u, wenn wir 


ut u ¥ 
0 0 


sehreihen und durch den Doppelpfeil tiber 1 andeuten, da® sich die Variation hinsichtlich u allein 
auf den Verschiebungsvektor beschrankt, 


: 
6, I = by ae, D + 8,u dO 
1l+yp : 
[7] [o(v)] 


= OF 


bzw. wegen 


t 
f(@S)dudO= f(Sdu)ndO= fVGoudv= f[(v S)u+ (SV) sul av 
[O(V))] [o(v}] [Vv] [VY] 


= f[(VS) 6u+6..6,D] dv 
[Y] 
die Relation 


iI = 6,|— [ive+ni+ Sa. Dien. 26] +2..2}a7| 
[Y] 


— | (@,—n©) dudo, 
[o(V)] 
aus der einerseits die Oberflachenbedingung 


38, —1 © = 0 auf O(V) (19) 


und andererseits wegen 


A.. dy (E.. D) E =A.. Cb, (E.. D) = [(A.. ©) E].. 64D 
schlieBlich 
E Ev a; 20 
22 [we 4+ ou + 6474+ apmarm (A..6)6| +43 av (20) 
[V] 
gefolgert werden kann. Setzt man hierin nun die natiirliche Bedingung (17) ein, so gewinnt man, 
da dann 


E Ey oe 
es Giada in Me 


wird, als weitere natiirliche Bedingung aus (20) die Gleichgewichtsbedingung 
/ VSG+f=0 in (V). (21) 
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Durch Elimination von A gemaB (17) erhalt man nunmehr aus (16) 


Opell = 0; 
ee [ [-ggl +) SG CR eel TGS) aA See dV 
(Vises 
= 3,ud0, 


[OY] 
wobei wieder zu beachten ist, daB £(u) bzw. 3,(u) bei der Variation hinsichtlich u nicht mitvariiert 
werden, und mit 
t 
6..9=(VW1r= Ley wW Sa. i Oe eee 
schlieBlich 
by,¢= 0 > 
iis | re [1 +) 6..6— (G.. G2 +2 Eo, (T—T,) (G..)] +(VS + bulav 
[Vv] 
a [81 —n ©] udod. 
[(7)] 
Da hierin die Variationsableitung hinsichtlich u als natiirliche Bedingungen sofort wieder die 
Gleichgewichtsbedingungen (21) baw. die Oberflachenbedingung (19) liefert, konnen wir die Wir- 


kungsweise des Verschiebungsvektors 1! in (22) als diejenige eines Lagrangeschen Parameters 
auffassen, so daB dem Problem (22) das folgende Variationsproblem aquivalent ist: 


dg I* =0, m= | frpld +») 6..6—»G.. 6)? +2 Ea,(T—T,) (G..S)]} dV (23a) 
[7] 

mit den Nebenbedingungen 
VS+t=0 in[V], %—n@=0 auf [O(V)]. (23) 


Von samtlichen, im Sinne der Gleichgewichts- und Oberflachenbedingungen méglichen Spannungs- 
zustanden macht der tatsdchliche Spannungszustand das Integral /7* extremal. 

Mit (23a) und (23b) ist nunmehr die Variablen-Transformation durchgefiihrt. Bevor wir auf 
Anwendungen iibergehen, in denen ausschlieBlich die Form (23) Verwendung findet, seien kurz 
die Eulerschen Gleichungen des Problems (23) verifiziert. Es ist 


fll | i [Q +) S..d66—» G.. S) G.. 6S) + Ea, (T— T,) (€.. 6S)] + (V 66) u| dV 
[V] 
4 | (8S n) udo, 
[o(V)] 
woraus mit Beachtung von 
(V 0) u = V (6 u) — (6S V) u = V (6S u) — 68.. 5 [V ou + (Vou) 


sowie mit 


(22) 


; [V(SuydV= f((6Su)ndO= f(6Sn)udo 
ease: V] [7] [O(7)] 
die Beziehung 


— (PP e—(t.6€.9—fer—- T))6|—31V eu +(Vou}..08 av =o, 


[Vv] 5 


also 


5 (Vout (V ony] = +4 8-5. O-PS ce T))G]=D 4) 


af . ‘ . . . . . . 
Ebenso wie der Spannungstensor selbst, soll seine Variation 6G ein symmetrischer Tensor sein. Hiermit 


ee me Giiltigkeit der Momentengleichgewichtsbedingung am Element auch fiir die Spannungsvariationen 
gefordert, 
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entsteht. Man gewinnt den Zusammenhang zwischen den Verschiebungen und dem (durch den 
Spannungstensor ausgedriickten) Deformationstensor und damit wegen! 


VX DxXV=@x Vl x V=Eiveut+ (vow x Vix Vv =0 
schlieBlich die Relation 


SxvixV—=I, 7 (€..8) —A* (Tt 1,)) 6 x V} x m0. (25) 


aus der nach verschiedenen Umformungen und mit Beachtung der Gleichgewichtsbedingung (21) 
schlieBlich die durch die Temperaturglieder erweiterten Beltramischen Gleichungen 


p 1 ae Z 
ee IV Vy C+ AoE + Volvo) | 
o = &.. ©, (=e Ye 

(26) 
hervorgehen. AbschlieBend hierzu sei noch erwahnt, daB das Variationsprinzip (23) auch unmittel- 
bar aus dem Prinzip der virtuellen Krafte? gewonnen werden kann. Es geht somit aus einer, die 
Vertraglichkeitsbedingungen umfassenden Energieaussage hervor. 


3. Beispiele fiir die geniherte Auswertung des Prinzips der virtuellen Krafte. a) Allgemeines. 
Bei den folgenden Beispielen handelt es sich um Warme-Selbstspannungsprobleme, bei denen 
weder Volumen- noch Oberflachenkrafte beteiligt sind. Mit {= 0 ergibt sich somit aus (23), 
wenn man dort noch gemafB 


©..€ =o, +0, +0,=0, 


Lie Po 5 = 27 
= (S.. 6? —S.. S] <0, Og 1 626-1 Oy Og —Aing 1 Tay + ty 5) =O (27) 
die erste bzw. zweite Invariante o bzw. o des Spannungstensors einfiihrt 
él1* =0, I = | [02-21 +») F +2 Bo, (T— 7) o] dV (28a) 
[Vv] 
mit den Nebenbedingungen 
VS=0 in(V), 8,=nGS=0 auf O(V). (28b) 


Zur Auswertung von (28) im Sinne des Ritzschen Verfahrens werden nunmehr fiir die jeweils in 
Frage kommenden Spannungen gewisse, den Randbedingungen sowie den Gleichgewichtsbedin- 
gungen gentigende und mit willkiirlichen Konstanten c, versehene Ansatze formuliert. Nach Ein- 
setzen in das Variationsintegral folgen die Konstanten als Lésungen eines sich aus den Extremal- 
bedingungen g//*/éc,, = 0 ergebenden linearen Gleichungssystems. Damit ist dann ein Spannungs- 
zustand gefunden, der den Gleichgewichts- und Oberflachenbedingungen sowie im Mittel auch 
den Vertraglichkeitsbedingungen geniigt. Die in den beiden folgenden Abschnitten b) und c) 
behandelten Beispiele wurden in dieser Weise untersucht, und die erhaltenen Resultate den exakten 
_Ergebnissen gegeniibergestellt. Das in Abschnitt d) untersuchte raumliche Problem hingegen zeigt 
eine mégliche Erweiterung des Ritzschen Verfahrens in dem Sinne, daf anstelle von Gleichungs- 
systemen fir Konstanten einfach lésbare Differentialgleichungen fiir die, die Spannungen approxi- 
mierenden Naherungsfunktionen aufgestellt werden kénnen. 


1 Die vektorische Multiplikation einer Dyade mit einem Vektor ist durch (a0b) X ¢= ao(6 x c) definiert 
(S. hierzu das bereits zitierte Werk von M. Lagally). Daher ist zunachst 


1 1 
Dx V=slVeow x V +(e V) x VI=5IVeG@x Vt tlV xVI= — => VelV xu), 
und damit weiterhin 


@ x VEX V=—SUV X wo VIX V=—F(V XW e(V x V)=0. 


2 K. Marguerre, Neuere Festigkeitsprobleme des Ingenieurs, Berlin 1950, 
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b) Achsensymmetrische Erwarmung einer gelochten Kreisscheibe (Abb. 1). Ein 
achsensymmetrisches Temperaturfeld ruft einen achsensymmetrischen Spannungszustand hervor, 
der sich, wenn wir die Scheibe als diinn annehmen, auf Radialspannungen o, und Tangential- 
spannungen o, beschrankt (ebener Spannungszustand). Fiir diese beiden 
Spannungsgréfen folgt aus der radialen Gleichgewichtsbedingung am Klement 


=< (ro). (29) 


Sind f3() den Randbedingungen f(0) = f(1) = 0 geniigende Funktionen, so 
erfullt der Ansatz 


o,(8) = 3 en SSPE) (30a) 


Abb. 1. 


die Randbedingungen fiir die Radialspannungen [o,(0) = o,(1) = 0], und man erhalt mit (30a) 
aus (29) den auf Grund der Gleichgewichtsbedingung fiir o, zwangslaufigen Ansatz 


ol6) = e000) = ae [(e—z +8) 0] = Dd eege 


ka 


o—F+8)8] = > af (30) 


ae 


Der Warmespannungszustand ist somit bereits bis auf die willkiirlichen Konstanten c, bestimmt. 
Fir diese folgt schlieBlich, nachdem man die Spannungsansatze in (28a) eingesetzt hat (wobei 
fiir das Volumenelement 


dV =2ardr= 2n(R—3 +2} dx = 2a (o—z +8) ag 
gesetzt und die Dicke der Kreisscheibe mit 1 angenommen ist, da sie sich als gemeinsamer Faktor 


in der folgenden Rechnung heraushebt), aus den Extremalbedingungen ¢//*/éc, = 0 mit den 
Abkirzungen 


1 
aoe = | (USP +J2) UP + £62) — (0 + 9) USP SEP + £8 f60 [e—y +8) aE 
0 
1 


(31) 
: ; 1 
B= | Ea d6) (62 +(e +s), O= TT 
j 7 J 
das lineare Gleichungssystem 
Pet 2: C, Oi, = O (i =1,2,..n). (32) 
=1 
Begniigt man sich speziell mit einer Approximation durch nur eine Funktion, d.h. setzt man 
d 1 
= afer), GW=a ge oe +6) 2 =, fete) (33a, b) 
so folgt aus (32) 
1 
[Bx 0 Ve +17) [o— > +8] ae 
Br P i t 2 
2 eet a 4 (33) 


1 
[UP +1920 + 6272] [o—4 4 &] as 
0 


Wir untersuchen nun den konkreten Fall einer unter vorgegebener Innen- bzw. AuBentempera- 
tur 0; bzw. @, stehenden Kreisscheibe. Fiir das Temperaturfeld innerhalb des Kreisringes machen 
wir, ungeachtet des tatsaichlichen logarithmischen Verlaufs, vereinfachend den linearen Ansatz 


HE) = Ou + F028,  te=16,48),  4d=0,—9,, (34) 


wahrend fiir f7?(€) mit Riicksicht auf die beiden zu erfiillenden Randbedineungen (0) = f(1) — 0 
ebenfalls der denkbar einfachste Ansatz, nimlich eine quadratische Parabel pee 


foe =é-& (35) 


ee © 6 
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und damit 
: d il i} F 
0) = |(@—z +8) E—8)) =0—5 + B—20) 8-38 (35) 


verwendet werden. Die hiermit durchgefiihrte Auswertung der in (33c) vorkommenden Integrale 
liefert schlieBlich, wenn wir konstante Stoffwerte E und «, voraussetzen, 


a Ba, AGE Ee 20'9* 


Fur die Warmespannungen hat man damit als Naherungslésung 
Ea A’ 1+ 20¢? 
6, & ey fr(é) =— SS EE EB), 
0. 9+200 a 
re a = Ea, Ad 1+ 20 0? eae ES, ») 
O, © e7 fo, (€) = — 9 9 + 20 02 Tes 9 3 (& “) p) 
wahrend sich exakt fiir denselben Fall 
272 
Eade |" @ ee 1\2 —(4) | 
64 — 
; (37) 
p—(4) | Pe es 
Eu Abo| | 20 | ieectr| ay ies (20—1+ 2 
ES) errr We tee ee logis 2o+1 
yeas, Teel In 
20 Q | 20—1 


ergibt. Die Formeln (36) und (37) sind fiir den Fall @ = (R/h) = 5 in Abb. 2 graphisch dargestellt 
(Naherungslésung gestrichelt). Man erkennt eine besonders gute Ubereinstimmung bei den Radial- 
spannungen, wahrend die Tangentialspannungen etwas schlechter approximiert werden. Der Fehler 
bei den Radialspannungen betragt nur etwa 5%, bei den Tangentialspannungen otwar EN. Die 
gréBere Ungenauigkeit der Tangentialspannungen erklart sich aus dem Umstand, daB f5?(é) aus der 
angesetzten Naherungsfunktion f3?(€) gemaB (30b) durch einen DifferentiationsprozeB gewonnen 


lon/é)| 
Ea, Ae 


i i i i ie die extremale Ringspannung 

de. In Abb. 3 sind die extremale Radialspannung |o, rein | sowie SSP : 

: Peiiencarande |O:imin| in Abhangigkeit von 9 = R/h aufgetragen. Man erkennt, daB die Uberein- 
stimmung der Naherungswerte mit den exakten Werten bei ed 0 sehr schlecht wird. Dies 
1aBt sich jedoch durch verfeinerte Naherungsansatze fiir for baw. fey, die z. B. noch die Kompa- 

- tibilitatsbedingung 


d — 
&, — g(r &) = 0 
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an verschiedenen Stellen der Scheibe befriedigen, sowie durch eine genauere Beschreibung des 
Temperaturfeldes beheben. 

c) Abkihlung einer diinnen Kreisscheibe. 1. Exakte Lésung. Ausgehend von der 
Differentialgleichung der Warmeleitung fiir eine an ihren Oberflachen Warme in ein Medium der 


Temperatur Null abgebenden diinne Platte! vom Halbmesser R und der Dicke h 
ae fiee nae, 2 _ 2% 38 
A RI eee Cy Boe ( ) 


erhalt man fiir den rotationssymmetrischen Fall mit dem Produktansatz 
@\2 
rif Ge) e (*) (39) 
aus (38) die gewohnliche Differentialgleichung 


af AGO ee ee i Be 40 
4-240 fe0, w=e+ bP, (40) 


deren Lésung mit Riicksicht darauf, daB in Plattenmitte die Temperatur endlich bleiben muB, 


f = CIg(x 1) = C Io(E @) » Q E=xR, (41) 


r 


= Be 


lautet. Der Warmeiibergang in ein Medium der Temperatur Null am Plattenrande r = R fordert 
unter Heranziehung des Newtonschen Abkihlungstheorems die Randbedingung 


1 0 A df 


OR, 3) ob elas lee Be TSE gree (42) 
woraus mit (41) die Eigenwertgleichung 
§ J,(8) =" I) = N Jl) (43) 


hervorgeht. Ihre Lésungsschar §,(N) bestimmt nach (40) auch die Abklingungswerte w = w,, 


& 1 RaR 
OR = Hi Sd ae meee ae (44) 
und es ist somit die Lésung (41) in der Form 
f= 3G, INEso) » 
also entsprechend der urspriingliche Ansatz (39) in der Form 
eras 
Ula: s) = De Ge Jolene a (45) 


verallgemeinerungsfahig. Aus der Anfangsbedingung 


He, 0) = Cy JolEe) = % 


erhalt man nun die C, als Fourierkoeffizienten der Entwicklung der Funktion @ = 9 nach den 
Besselschen Funktionen J,(é, 0) in der Gestalt 


ek 
2 Jo(z) dz 

9-4 fe LJe(Sn) + JIE)” 
f z Jie) dz 
0 

so daB schlieBlich fiir das Temperaturfeld 
99, $1 SiG) Jono) (4) 
i) . =) Sk) YO\Sk a 
(= 2% DS STG) FRG ° (#7) 


hervorgeht. 


_* Melanu. Parkus, Warmespannungen stationarer Temperaturfelder, Wien 1953. Es bedeutenc die spez, 
Warmemenge je Gewichtseinheit und 2 bzw. « die Warmeleit- bzw. Ubergangszahl. 
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Zur Ermittlung der zugehérigen Warmespannungen geht man von den Gleichgewichtsbedin- 
gungen des ebenen achsensymmetrischen Problems 


0 
ap fOr) — % = 0 (48) 
und den thermisch-elastischen Spannungs-Dehnungsbeziehungen ! 
ov 1 1 
é =a = 7 (O,—¥ 0) + (P —B) &===7(—r6,) +o, (9 —%), (49 a) 
bzw 
E ov v E é 
a= la ty 7 —O +9)a0—%)|, =p Let? +9) 4 (0—8)] (9b) 


aus, und erhalt, indem man (49b) in (48) einsetzt, schlieBlich die Verschiebungsgleichung 


xlr ee} = 0+) a5. (50a) 
mit det Lésung : 
u(r, #) = (I +9) a= | r B(r,t) dr + Cyr + SO | (50b) 


Unter Beachtung von 
P R 
| 1 SolGie) dr = 52 Alé.o) 


gewinnt man nun aus (50b), wenn man (47) einsetzt und die Integration hinsichtlich r ausfiihrt, 
nachdem man noch die willkiirliche Zeitfunktion C,(t) wegen v(0, t) = 0 gestrichen hat, 


e a RV Seo Idee) — (2) 
u(r, t) = C(t) r + 2 (1 +7) a, 9 RD, BED ETE ° : (50c) 
so daB fiir die Radialspannungen nach (49b) schlieBlich 
O; ==; (G0 + 04, I) — (1 —») [roar 
. 3 co * (0) 2 
Pees eh os debe R ul) (Eve) — (“E)"* ea 
ae Ey C(t) as OX, ON 2 (1 v) r a ER [J2(Ex) we J2(E,)] e ( 0 ) 


hervorgeht. Aus der Randbedingung o,|r,, =0 folgt nun C,(t), und man hat endgiiltig die Warme- 
spannungen 


fe) E J1(&;, 0) — 1) Ji(&) eee os 
0,(Q, t) =— 2 Ea, 9, >, 7 ea ee é 


k=1 


(51) 


eye 

oo ne Oe (= Ji(Sx 0) + n0)| SiEn) (2H)*« 
fe &, \ 0 

(0,1) = —2 Ex Dy D! &, iG) + JiGol : 


k—1 
Die Reihenentwicklungen (47) und (51) wurden fiir eine Stahlplatte vom Halbmesser R = 50 cm 
und der Dicke h = 5 cm mit den Materialwerten c, ~ c, ~ c = 0,11 cal/g °C, A = 0,111 cal/emsek 
°C, y = 7,85 g/cm, dXstani-Lutt = 20 keal/m? h °C = 0,556 - 1073 cal/cm? sek °C zu verschiedenen 
Zeiten t ausgewertet. Es ergaben sich die in Abb. 4, 5 und 6 dargestellten Kurven. Man erkennt, 
da® die in Plattenmitte auftretenden extremalen Radialspannungen bzw. die am Plattenrande 
herrschenden gréBten Tangentialspannungen, von Null aus anwachsend, im Laufe der Zeit einen 
Extremwert erreichen, um hernach wieder langsam abzunehmen. Die Extremwerte werden fiir 
beide Spannungen nicht gleichzeitig erreicht, sondern z. B. fiir die gréBte Tangentialspannung 
~ 6,(1, toc,) = 90,0326 Ea, 9 zur Zeit too, ~ 1200 sek = 20 min und fiir die extremale Radial- 
spannung 0,(0, toc,) = — 0,0122 Ea, J aur Zeit too, © 2100 sec = 35 min. 
2) Naherungslésung. Wir gehen wieder von (28a, b) aus, also mit dV = 2ahrdr,o =o, + Op 
6 =0, 0, von 


t= Fr [(o, +0)? —2 (1 +) 0,0, +2 Ex 0 (o, +0, dr = Extr., (52a) 
r=0 


1 Der Ausgangszustand # = @ sei spannungslos. 


~ 
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wobei als Nebenbedingungen die Gleichgewichtsbedingung 
7) 


Bp (Or) —% = 0 (52) 
sowie die Oberflachenbedingung 
o,|z,2 = 0 (52 c) 
zu beachten sind, und da wir aus (52b, c) noch 
R R 
R 


r= 


f ra) 1 
[rocar= | rong (re) dr = (ro,))| = 0 
0 0 


folgern kénnen, verbleibt schlieBlich als maBgebende 
Beziehung aus (52a) 


: ) 
*= frl{(o, +o)? +2 Eo, d (o, +.0,)] dr = Extr. 
tre (52d) 


Zur Erzielung einer gré®eren Genauigkeit wollen wir 
von den Spannungsansatzen nicht nur die Erfillung 
der Bedingungen (52b, c) fordern, sondern noch zu- 


48 _ B(0t)-3(1t) 


4, % 


G05 


2400 sek 
Abb. 4, 


901 


__%(04) __ (at) 
Eu,d, Fo, 


001 


2400sek 


40 


Abb. 5. 


Abb. 6. 


satzlich die Erfillung der Vertraglichkeitsbedingung am AuBenr 


sie mit einem die Gleichgewichtsbedingung (52b) befriedigenden Naherungsansatz von selbst 


ande [in Plattenmitte ist — 
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erfiillt]. Die Vertraglichkeitsbedingung des ebenen rotationssymmetrischen Problems lautet 
bekanntlich 


o 
a, (re) —€ na 


rT 


woraus schlieSlich unter Beachtung des Hookeschen Gesetzes (49b) und mit (52 b) 
A (ro) —¢, + Bar =0 (52e) 
hervorgeht. 
Wahrend (52e) wegen (52b) fiir r = 0 stets erfiillt ist, ist dies fiir alle anderen Plattenpunkte 
zunachst nicht der Fall. Wir fordern daher von den Naherungsansatzen zusatzlich fiir den AuBen- 
rand r = R die Erfillung von (52e), also wegen o,(R, t) = 0 


Gor od 


o,R,t) +R = fis E«, ares ae 05 
bzw. mit 
ob 
Fa = — K(t) (52 f) 
die Relation 
7) 
KZ A | = fa, Ke (528) 
(Oma 
Fiir o, wird nun ein Naherungsansatz in der Form! 
0, = C(t) + ¢2(t) o + eg(t) 0° (53a) 
benutzt, womit auf Grund der Nebenbedingung (52b) auch der Ansatz fiir o, feststeht: 
é Q 
Oe op (FOr) =], (0 Or) = Hor 3 cag? +A cQ0? (53b) 
bzw. 
Fi = 6 e0 + 2, o* (53¢) 


Von den hierin enthaltenen drei willkiirlichen Koeffizienten cy, c, und cz sind nun an Hand der 
beiden Nebenbedingungen (52c) und (52g) zwei in Abhangigkeit der dritten ausdriickbar. Ein- 
setzen von (53a) bis (53c) in diese beiden Nebenbedingungen liefert das Gleichungssystem 
Cy + C3 =— Cy, 9c, +16c,=—cq,+ Ea,K, 
aus dem sich schlieBlich 
15 1 8 1 
C= —7%—z Ea, kK, C= 7% +a Eo, Kk (54) 


ergibt. Damit lauten also die Spannungsansitze 
0,(0,1) = 2 (715 2 + 8¢%) —5* Ka) —o), 
(t) Ea (3) 
C 
o(o, t) = 2 (7 — 45 ¢? + 320°) — = Kip) Bo? — 40%), 


und der hierin noch enthaltene Koeffizient c,(t) wird aus der Extremalbedingung 0L/*/dcy = 0 
bestimmt. Man erhalt schlieBlich aus 


1 1 

d r d r+ o 
[eG +o cet ee dg + Ba, | 9 00 ) do =0 
=0 e=0 


mit Verwendung von (55) 


col!) = qq Ea, K() — fq Em | (70 — 30 0? + 20 94) (0, t) do. (56) 


“a 


0 


1 Der Ansatz einer linearen Funktion von g kommt wegen der Rotationssymmetrie des Problems, die 
20,/20|o=0 = 0 bedingt, nicht in Frage. 
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Um zu einer einfachen Naherungslésung zu kommen, wird 0(Q, t) nicht in seiner exakten Form 
(47) eingesetzt, sondern auch fiir das Temperaturfeld selbst wieder ein Naherungsansatz verwendet. 
Wir setzen niherungsweise 
Ho, t) e Ap(t) + Ag(t) @ + A3(t) g° (57a) 
und driicken die Koeffizienten Aj, A, und A, durch die Temperaturen (y(t) bzw. Pp(t) in Platten- 
mitte bzw. am Rande sowie durch K(t) aus. Man erhalt aus den Bedingungen 

| oo 

Do = 0,t =e Ou(®) ? J =1,t See I(t) ? R er 
schlieBlich 

Ags= Dips A, =—3(98y— Or) + K=—3A04+K, A, =240—K, Ad = By — Pr, 


r=R,t 


und damit in erster Naherung fiir das Temperaturfeld 
H(o, t) ~ Vy(t) — Ad(Y) [3 o — 20°] + KY [e* —@'] - (57b) 
(Eine erste grobe Abschatzung, zu welchen Zeiten (57b) als Naherungslésung angesehen werden 


kann, wird durch die Aussage (0, t) < Py(t) geliefert. Nachdem man fir 0(0, t) nach (57b) eine 
Extremalrechnung hinsichtlich 9 durchgefiihrt hat, erhalt man schlieBlich mit (58e) die Aussage 


/ N 


In unserem Falle ist N = 0,25, so daB sich die Abschatzung yy = 1,08 Bp ergibt. Aus der Lésung 
fiir das Temperaturfeld entnimmt man, daB dies ungefahr fiir t = 900 sec erfiillt ist. Die Naherungs- 
lésung fiir das Temperaturfeld stimmt sogar schon von t = 600 sec ab mit den exakten Werten 
nahezu iiberein). Einsetzen von (57b) in (56) liefert nach Integration 


e(t) = — 2% [7 40) — Ko). (576) 
so daB schlieBlich fiir die Spannungen aus (55) 
(ort) » — 5% AGW) [7— 15 0 + 805] + 5K) L—5o*§ + 405], (58a) 
oo, t) ee — 2% AG) [7 — 45 o* + 32.05] +2 K(e) [1 — 15 g* + 1605] (58b) 
hervorgeht. Fiir die Extremalspannungen! folgt damit 
min o,(t) = 0,0, t) = 0,0, t) = —5 [7 A6@) — KO], (586) 
max o,(t) = o,(1, 1) = 4% [6 Ad@) +2 Kio). (58d) 


wobei K(0) = 0 sowie fiir K(t) im Sinne des Newtonschen Abkiihlungstheorems 


ob R R : 
KW) =—5 ae OL, t) = Ox) =NOR(t) far 2> 0 (58e) 


zu setzen ist. 

Die Ausdriicke (58a, b) sind zum Vergleich mit den tatsachlichen Spannungswerten (51) in 
Abb. 7 aufgetragen worden. Wahrend die Genauigkeit fiir t = 300 sec nicht mehr als befriedigend 
angesehen werden kann, was jedoch im Hinblick auf die Einfachheit der Naherungsansatze in 
Kauf genommen werden muf, erkennt man schon von t = 900 sec ab praktisch iiberhaupt keine 
Unterschiede mehr. Bis auf ein verhaltnismafig Kleines Zeitintervall nach Beginn der Abkiihlung 
gibt also die Naherungslésung den tatsachlichen Spannungszustand mit ausgezeichneter Genauig- 
keit wieder. Die (ebenfalls zur Zeit tps, ~ 1200 sec auftretenden) griéBten Tangentialspannungen 
ergeben sich mit ; 


A (toed) _ 9,954.07 K(to.,) = NO 0 
om Tew) > (toc,) =S R(too;) — ,170 416 Do 
schlieBlich zu 
max 0; it 
Ea, > 20 [6.0,054077 +- 2.0,170416] = 0,033 262 , 


-wahrend der Wert aus der genauen Lésung 0,0326 betrug. 


* Wie Abb. 7 zeigt, nimmt das Extremum in Plattenmitte erst von t ~ 900 sec ab einen gegeniiber dem im 
Plattenbereich bis dahin vorhandenen zweiten Extremum gréBeren Wert an. 
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Kine weitaus einfachére, wenn auch grébere Naherung, die aber trotzdem die richtige GréBen- 
ordnung der entstehenden Warmespannungen abzuschatzen erlaubt, gewinnt man, wenn man 
in (52d) die einfachsten Ansatze 


0, = ¢,(t) [1 — eo]. 
8 = F (0%) = a(t) [L—3 4), 
Yo, t) = Py(t) — Ad(t) 0? 


einfihrt. Man erhalt schlieBlich aus 
OlI* /écy = 0 


ot) = — 2 40(0) 
und damit 
o(Q,t) » — =" Ady) 1 —e'}, 
olo,t) » — = Ad() 1—3 eI 


mit den Extremalspannungen t=7200 sek 


min o, = 0,(0, t) = min o, 


= 9,0, t) =— 7 Ea, AH) , 


maxo, = 0,(1, t) = a Ex, A(t) . 
07 

Diese sind, im Gegensatz zu (58c, d), allein : 
von der Temperaturdifferenz abhangig, so 
daB mit dieser Naherung die extremalen 
Spannungen zu demjenigen Zeitpunkt ft, 
berechnet werden, zu dem A(X(t) ein Maxi- 
mum wird. Um ohne eine Auswertung der 
Reihenentwicklung fiir das Temperaturfeld o 
zur extremalen Temperaturdifferenz A? 
= Oy — 0p zu kommen, betrachten wir 
anstelle der alternierenden Reihe 
AKt) (0, t) — AL, t) 

ae 3, 


2 [1 — Jo(Ei] Ju(Ex) ,—(Z*) * 
§, Ja) + 2G) 


02 


t=7500 sek 


t =2400sek t=2100 sek 


of 
Ot 


or 
aa MT | 
Ea,9, Eo, 


Abb. 7, Darstellung von 
Ok \? 
(— 1}F opal A; a (=): (Naherungslésung gestrichelt), 


als Naherungslésung nur ihre beiden ersten Glieder 
@, 2 pai Oy 2 
Lo ue aah 
_wobei fiir den hiermit begangenen Fehler R(t) 


? 


RO) < 4ge7) 


: gilt. Aus einer Extremalrechnung fiir den Naherungsansatz erhalt man schlieBlich, daB zur Zeit 


a A, [@,\* 
b= a—am lala) 


fiir die extremale Temperaturdifferenz 


(@,/«,)? 


max A? WW, a A, ‘Wy 2) 1— (@,/@,)? 
= [1—(G)| fala) 


hervorgeht. 


a 


404 R. Trostel: Genaherte Berechnung von Warmespannungen Ingenieur-Archiv 


Um den hierin enthaltenen Fehler abschatzen zu kénnen, berechnen wir uns eine Zeit t)() aus 
der Relation- ‘ 


— (2) 
A, e Sey ek, i 
@, \2- @, \2— 100 
ee 
zu 
P a 100 A, , A, 
EOS: wz — w? | a Ag a 


und wegen W3 > Wz, also 


le) a ee 
d.h. 


@, 2 
IRG)) aoa _u 
@, \2- @, \2— @, \2- ©; \2>) Selo 
= Ae Ise Se EEE a Ne 
Aye a eg! () Aye a) — Age & 

folgt, daB zu diesem Zeitpunkt t, das vernachlassigte Restglied R(i)) kleiner als ~ % der Summe 
der beiden ersten Reihenglieder ist. Setzt man also i, = tg, so folgt fiir die Prozentzahl y des Fehlers, 
mit dem die naherungsweise berechnete extremale Temperaturdifferenz behaftet sein kann, 
As 1 
A, (w;/01)?—1 ° 


b= 100 
Mit den im vorliegenden Falle sich ergebenden Zahlenwerten 
A, = 0,1206, A, = 0,1143, A, = 0,02346, (w,/a)? = 0,279-1073 sec, (w/a)? = 1,029 - 1073 sec™ 


hat man damit 


0,279/1,029 
max Ad 0,279] [0,1206 0,279 1 — (0,279/1,029) 
Ses pre? > 2 = 0,05 
Do 0,1206 L 1035 bars 7,029] DoS 
mit einem Fehler von héchstens 
= 190 102346 1 — 16%, 


0,1143  (1,029/0,279) —1 
was in Anbetracht der groben Approximation fiir die Spannungen als ausreichend angesehen werden 
kann. Es folgt somit 


min 6; 


Ea On 


max Oo; 0,0553 
Ea 6 ~~ 2 


ee 2 0,0553 = — 0,01383 , = 0,02765 , 


wahrend die Werte aus der genauen Liésung 0,0122 bzw. 0,032 6 betragen. Es sei noch vermerkt, 
daB auch diese grobe Naherung fiir groBe Zeiten t die tatsichlichen Werte sehr gut wiedergibt. 


d) Naherungslésung fiir die Warmespannungen in einem Hohlzylinder mit 


freien Enden. Wir betrachten einen Hohlzylinder der Lange / mit dem Innenhalbmesser r; und | 
dem Aufenhalbmesser r,, dessen Temperaturfeld ((r, t) speziell von der Zylinderachsenrichtung z _ 


unabhangig sein mége, was z. B. bei an den Zylinderenden isolierten Rohren mit von z unabhangigen | 
thermischen Einwirkungen auf den Mantelflachen der Fall ist. In erster Naherung identifizieren — | 


wir den entstehenden Spannungszustand mit denjenigen Spannungen o, of) und of, die sich — | 


infolge des Temperaturfeldes im mittleren Bereich eines unendlich langen Rohres mit freien 
Zylinderenden ergeben, und fiir die 


or, t) = — 22 ; xa) i fr oO dr— f = (2) fr o | ; (59a) 


of)(r, t) = — — fo —- (GE 1 fr 0 dr + i + (F) fr O a ; (59b) | 


_ be 
jt 


ro | (59c) 
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folgt?, wobei insbesondere wegen der radialen Gleichgewichtsbedingung des ebenen achsensym- 
metrischen Spannungszustandes 


0 
ay (rol) = of (60a) 


und mit Heranziehung der partiellen Integrationsregel bei Beachtung der Randbedingungen 


or, t) = 0, a Mr, t) = 0 (60b) 
die Beziehungen 


fr o) dr = Be oY dr = r2 or, iy f ro” dr =0, f: ro o® dr = 0 (60c) 


r, r 


gelten. ; 

Fiir ein endlich langes Rohr werden mit (59) jedoch die Randbedingungen o, = 0 fiir z = + 1/2 
nicht erfillt, und es muB daher, um die Randspannungsfreiheit an den Stirnflachen herzustellen, 
ein zusdtzlicher achsensymmetrischer Spannungszustand 2 herangezogen werden. Die Gleich- 


gewichtsbedingungen fiir seine Spannungskomponenten o?), o(?), o® und 72) lauten 
0 é a] 0 
peo) 0) ae) = 0, a Oar) = 0, (61a, b) 


und diese haben nun folgenden Randbedingungen zu geniigen: 


OD rs = +1/2,t — = on(r, t) , (62a, b) 
Ol, seus. oo Os, i ie rh | Bee == OG (63a—d) 
OP parpnt =, G2 5 eee Ore (64, b) 


Wir verwenden fiir o’) den Ansatz 


P(r, 2, t) = —oM(r, ) fC); C= a : (65a) 
wobei f(¢) mit Riicksicht auf (62a,b) eine zunachst willkiirliche, jedoch den Randbedingungen 
fa )=1 (65b) 


genigende Funktion bedeutet. Damit hat man aus der Gleichgewichtsbedingung (61b) unter 
Beachtung von (60c) fiir die Schubspannungen 


= 2 (frowan sid +222 <2" ow pg + 462. (66a) 


Den Randbedingungen (63c, d) kann mit (66a) wegen (60b) durch A = 0 geniigt werden, und 
damit verbleibt also 


i) — 2-0 f'(C)., (66b) 
wobei mit Riicksicht auf die Randbedingungen (63a, b) fiir f weiterhin die Randbedingungen 
f(t 1) =0 (66c) 


zu fordern sind. Aus der Gleichgewichtsbedingung (61a) folgt, nachdem man auch fiir of”! mit 
- einer zunachst willkiirlichen Funktion g(¢) einen (die Randbedingungen (64a, b) befriedigenden) 


t tz in der Form 
Produktansatz in ler af a oi ie i 
benutzt, fiir die Ringspannungen 


7] 7) 
gues. (ro) et) (68a) 


mit Beachtung von (60a) und (66b) schlieBlich 
of?) = oft e(C) + 4(7) ofS"). (68b) 


1 S.z.B. R. Trostel, Ing.-Arch. 24 (1956) S. J. - 
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Die endgiiltigen Spannungen sind damit 
o, =o! + of?) = of [1 + a(C)], 
o, =o) + of =o [1+ 9(0] +4 (7) Pre 
a, = ol) + of) = o™ [1—FCMs 
r= et =2(7)o FO, 


Sie beschreiben einen, den Gleichgewichtsbedingungen und sdmtlichen Oberflachenbedingungen 
geniigenden Spannungszustand. Einsetzen in (28a) liefert mit 


o=0,+06,+0,=o0 [2+¢ Sister yea , 
T= Or yee Oso; tre =O Af) (+g) ol of (seg, 
+4(7) Mom sa ay joe fd + g) ey | 


wenn wir zunachst die Integration iiber r ausfiihren, was auf Grund der bekannten Funktionen o; 
of und of ohne weiteres moglich ist, ein nur noch die Integrationsvariable ¢ pn thalredaee 


Variationsproblem von der Form 


oe 
JTF (C), 8(2)] a Fle(2),f(0).f'(0),f'(0)] dC = Extr. (70a) 


(69) 


(1), 


Hierbei ist (mit Herauskiirzen des gemeinsamen Faktors z 1/2 E) 


F=2a,f'®—a,[f"(1 + g)—f"]—a,f" +a [2 +g—f)—2 (1 +7) (l—f) (1 + 8)] 


—2(2+8—f)4;, (70b) 
wobei 
r: Te = Tag et 
eae eh irene = Hasty [ (Pf roprar, 
Sag : 2 (70c) 
a im 8 a Ta 
a =—8 Eo, | (7) rold-dr 2-2, a= | roar, eel eae 


bedeuten und insbesondere mit (59c) 
rT 4 r, 2 
an a a 2 2 a 
froa—zha(froa) lays, (70d) 
folgt. 


Durch Nullsetzen der Variationsableitung von (70a) erhalten wir die Eulerschen Differential- 
gleichungen 


oF oF d oF d? OF 
Fee F mee ae ye (71a, b) 


Zunachst folgt aus der ersten Gleichung (71a) mit F nach (70b) | 
Sg tag t= pad ee (72) 
Weiterhin ist mit Beachtung von (72) 


Fry) a; +24 (1—99) (f—1) ty af" =2a,(L—*) ft vas", (73a) 


of aa HE of 7a” (73b) 
oF fy a3 a, a 
of’’ =f (4 1—7,,| 7 a taf), 


@ oF (73¢) 


az 
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so da8 damit aus der zweiten Eulerschen Gleichung (71b) 


f(a %—5e,\—-20—9) anf" + 2a, (1—*) f =—2[(1 +») a,—(1—2*) a] =0 (74a) 


a, 
hervorgeht. Mit den Koeffizienten 


vee 2 (1 —») a, a, 4 (1 —?*) a? 
+ x ee Nes a Bt 


Syne 
8 a, a, — a2 


, c= (74b) 


8 a, a, — a3 
laBt sich (74a) auf die einfache Form 
fe aT! + § f= 0. (74c) 


bringen. [c? bzw. cj sind stets positiv, denn es sind nach (70c) a, > 0 sowie a, > 0. Zum Nach- 
weis, da der Nenner der Koeffizienten positiv ist, benutzen wir die Schwarzsche Ungleichung. 


Man erhalt schlieBlich mit Beachtung von (59c) und (60 c) 


Ta, 4 Ta tH 
8 a, a, — a5 = 64 f (7) rot de f of!) dr —(1 +9) [rota 


i i 


Ta 4 Te Ta 
= 64 J (7) r g(t) dr (J r of)? dr — » (2 +) frog? i) > 0 5 


da o{”” bis auf einen kleinen Bereich, in dem das Vorzeichen von o) wechselt, erheblich gréBer 
als o” ist]. Der Ansatz 
f= Ge : E (75a) 
liefert nach Einsetzen in (74) eine biquadratische Gleichung fiir die charakteristischen Exponenten, 
aus der man schlieBlich 
Ai, 4 4 se Nee 9) (75b) 


erhalt, wobei 
1 


Ye 


ist. Hierbei wurde c} => ct vorausgesetzt. Im Falle ct = c§ sind samtliche Wurzeln reell: 


ti ———————— a 
ele 2 yeR + e?, Re, Vo3—e} (75) 


4 


i ct yt—ch = + m,, (75d) 
+Vei—fa— 4 = + py. 


Dementsprechend lautet die allgemeine Lésung fiir ct <c} mit den Integrationskonstanten 


A,, As, B, und B, 


Bs 4 (: Ee 
M1234 =EVit fe—d = 


f(© = e! (A, cos x6 + Ay sin xy) 4+ e— 5 (B, cos x6 + B, sin x2 C) , (76a) 
wahrend im Falle ct = c} 
f(©) — A, ela S a Ay es + Ay PEAS + A, em aS (76b) 


hervorgeht. Mit Riicksicht auf die Symmetrie des Problems zur Ebene z = 0 verbleiben allein 
die in € symmetrischen Lésungsanteile, also_ 


fare; Ss cf: f (0) = A, Gof #6 cos x6 + A, Sin x, C sin x6 , (77a) 
_ fir cf Sc: f(C) = A, oj ao + A, Gof pee - (77b) 
Nach Bestimmung der Konstanten aus den Randbedingungen (65b) und (66c) hat man dann 
schlieBlich 


VRE Cof 4 sin % ++ Y%g Sin H1 COS Xe (776) 
: "4 1 #1 SiN %y COS Hq + H%_ Sin % Cof x4 
eee 4 ‘ é 
fiir cy S ¢3: , Hy Cof #1 Sin %,— % Sin x Cos xy (774) 
2 941 SIN %y COS Hy + H%_ Sitt xy Cof x ” 
es fg Sit fy see : (77e) 
ee 1 be Sit fg Gof wy — oa Sin pr Gof My 
eect Age pa Ot (77f) 
2 


~My Sith fog Cof wr — fr Sit pra Cof ry” 
, 99% 
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Berechnen wir nun noch aus (72) 


as 


a= f" of, (778) 


so ist damit der Zusatz-Spannungszustand vollstandig ermittelt. Selbstverstandlich sind in dieser 
Lésung samtliche Grenzfalle enthalten: 


«) Die dimne Kreisringscheibe. Aus (74b) und (70c) ersehen wir, da wir die Koeffizien- 
ten c? baw. c$ auf die Form 
2 (1 —») a, a, - ey eS 
pruere Ss o eet ee eae A—ath 
Mode, ie ee °2 Biggs tt 2 


bringen kénnen. Damit werden fiir cf < ¢3 


Leaps Se “ be wees = 
ecrmns 2 Py oe , ET ene. 
Maa Vote? = 1A, Hy = —— ce? — ci = 1, 


und entsprechend fiir c} > ¢} 
2 =i; Se 


wobei die 4, und f/, nun in entsprechender Weise nur noch die quergestrichenen, von! unabhangigen 
Konstanten ¢, und ¢, enthalten. Fiir /—> 0 erhalten wir nun mit sin %2 & %9, COs % & 1 und Sin x1 
we x4, Cof x, ~ 1 aus! (77c, d) Ay=1,A,=0 und damit aus (77a) mit oj x6 ~ cos xo = 1 
schlieBlich f(¢) = 1 und aus (77g) g(¢) = —v. Die Spannungen sind damit nach (69) 


6=0, +, =0; -¢,=o0%(1—?), = og, =of (l= 9). 


Axial- und Schubspannungen verschwinden, und die Radial- und Ringspannungen nehmen die 
Werte des ebenen Spannungszustandes an. 


6) Das unendlich lange Rohr. Mit 1—o0o streben die Werte x, 
und x, ebenfalls nach Unendlich. In (77c,d) kénnen wir nun 
Sin x, ~ Cof x, ~ e%/2 setzen und erhalten 


emer e! (x1 sin %_ + %2 COS 2%) 
= 


5 x 2 
piece #1 SIN %5 COs x sega 
Ug —> co 1 2 tes 4 Y 


und entsprechend 4A,—0, so daB als Lésung f(2) = g(¢) =0 und 
somit fiir die endgiiltigen Spannungen nach? (69) 


Co, = oft) ) (Oy of) ? One of) 2 Tz = 0 


entsteht. Es ergibt sich also die Lésung des Ausgangszustandes. 


y) Diimne Rohre. Ist die Wandstarke h =r,—r,; klein ge- 
gentiber dem mittleren Rohrhalbmesser R = (r, + r;)/2 (Abb. 8), 
so liegt das Problem der Kreiszylinderschale vor. Wie man zeigen 
kann®, streben z. B. im Falle einer Temperaturdifferenz zwischen 
Innen- und AuSenmantelflache die (naherungsweise nach einer qua- 
dratischen Parabel verteilten) Radialspannungen o!!) mit h/R gegen Null, wahrend die Axialspan- 


nungen o{") den Charakter von Biegespannungen annehmen. Setzen wir also mit? ¢, = E «,/(1 — ») 
naherungsweise 


Abb. 8. 


0) AG h ae a5 
Me) =—F~[L—2z)|, ofa) = add z, 


? Dasselbe folgt (da ja c{ = c$ = 0) auch aus (77b) mit (77e, f). 
* Auch hier folgt dasselbe Resultat aus der fiir cf = c4 giiltigen Lésung. 
3S. z. B. R. Trostel, Ing.-Arch. 24, 1956 S. 1. 
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so erhalten wir aus (70c) 


lap + h/2 
én 8 | 75 go)" dr wz 8 RS o)"(x) dx = : ot(A8)? h3 R3 , 
rj «= —h/2 
+ h/2 
% =8(1+») / 3 ot) dr aw 8(1 +9) RB | olt"(x) dx = *A* (49) WR, 
x = —h/2 
eee 
iis [re of dr aw olt"(x) dx = 55 o%(A0? Rh, 
x = —h/2 
also aus (74b) mit 
- _» (6) 49h RY! h\t] 2 
SECT eel Rr oar [10 — (ev) (Fe) | RY Ge (OAV h R)* 
schlieBlich 
5) = 5 (1 — v) 101—+ /R\? 1 (l—»?)? 
enna” a= —3mR = 16 (a) >4=16 RS 
Aus (75c) hat man dann 
= Zz es, 4 /2,5 (1 =) 
ty RY hy wn ww Vesa, Se ame 


4/7 RE 
ae Veet a = 0,795 |/ 7. (78a) 


Aus einer bekannten Naherungstheorie fiir den rotationssymmetrischen Biegespannungszustand 
einer Kreiszylinderschale folgtt 


Q= jee Tay = 0,761 poe said (78b) 


was mit dem obigen Grenzwert nach (78a) fir den Abklingungskoeffizienten praktisch tberein- 
stimmt. 


4, Zusammenfassung. Mit Hilfe von einfachen Naherungsansiatzen, die in die, einem jeweiligen 
Problem zugehérigen Variationsintegrale eingefihrt werden, kommt man im allgemeinen zu aus- 
gezeichneten Approximationen. Neben der dem Ritzschen Verfahren entsprechenden Methode, 
Naherungsansatze bis auf willkirliche Konstanten vorzugeben, die dann aus Extremalbedingungen 
berechnet werden kénnen, scheint vor allem das im letzten Abschnitt geschilderte Verfahren, das 
auf einfach lésbare Differentialgleichungen fiir die Naherungsfunktionen fiihrt, eine besondere 
Bedeutung auch fiir andere Probleme der Elastizitatstheorie zu haben. 


(Eingegangen am 14, Dezember 1959.) 


Anschrift des Verfassers: Privatdozent Dr.-Ing. R. Trostel, Berlin-Charlottenburg 2, Hardenbergstr. 34, 
Lehrstuhl fiir Mechanik der Technischen Universitat, 


1S. z.B. I. Szabé, Einfiihrung in die Technische Mechanik, 4. Auf!., Berlin/Gottingen/Heidelberg 1959. 
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The Stability of Motion of a Rotor with Unsymmetrical Shaft 
on an Elastically Supported Mass Foundation 


By A. Tondl 


1. Introduction. Several authors have made study of the dynamics of rotors with unequal 
shaft stiffness in two mutually perpendicular directions. In almost all their work the solutions 
are performed under the assumption of perfectly rigid bearing supports upon an absolutely rigid 
foundation (see for instance’ » * 4 5 678), The author of this paper has also solved the 
case of a rotor mounted upon elastic, massless supports” 1°. 

With the tendency to construct light and more elastic foundations for turbogenerators (parti- 
cularly low-tuned steel foundations) the question arises in the case of two-pole generators whether 
the stability of the rotors is not adversely affected by low tuning of the foundation. 

The solution to be presented in this paper is carried out on a simplified system; several examples 
are used for elucidating the effect of the foundation tuning upon the stability of the rotor motion. 


2. Solution of the Problem. Let us assume a one-disk rotor with mass m mounted on a massless 
shaft whose section at the place of the disk attachment has stiffnesses c, and c,. The rotor rotates 
with a constant angular velocity w. For the sake of simplicity let us neglect the elasticity of the 
bearing supports. Let us further assume that the mounting of the rotor upon the foundation can 
be represented in a simplified form by the scheme of Fig. 1, in which the foundation can be replaced 
by a mass mounted elastically upon massless springs whose total stiffness in the vertical direction 


Fig. 1. Scheme of the system. 


is cy, and in the horizontal direction cy,. Let us also assume that the disk is located on the shaft 
midway between the bearings and that the foundation is completely symmetrical relative to the 
disk plane so far as both distribution of the mass of the base plate and the rigidity of the mounting 
is concerned; thus in standstill position the disk centre of gravity (assuming perfect balance) lies 


: a asl, Bane a Gasturbinen, 5. Aufl., Berlin 1922, p. 931. 
. B. Biezeno an . Grammel, Technische Dynamik 2. Aufl. Bd. 2, Berli Otting i o 
3 R. Soderberg, (1931) ASME paper APM-54-4. a eee 
* J. Dick, Phil. Mag., 39 (1948) p. 946. 
° F, M. Dimentberg, Izdat. ANSSSR, 1959. 
° 0. N. Romaniv, Lvovskij politechniéeskij institut, Lyov, 1957. 
 H. D. Taylor, J. Appl. Mech. 7 (1940) p. 71. 
° W. Kellenberger, Ing.-Arch. 26 (1958) p. 302. 
° A. Tondl, Rozpravy Cs. akademie véd, Yada TV, No 5; 1955, 
0 A, Tondl, Vydavatelstvo SAV, Bratislava, 1958, 
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on the same vertical axis as the centre of gravity of the foundation mass m,. This mass encompasses 
both the mass of the base plate and the mass of the stator, stands etc. Under these assumptions 
the gyroscopic effect of the rotor need not be taken into consideration. 

Let us denote the coordinates of the disk centre of gravity as x,y and those of the centre of 
gravity of the foundation mass as Xp, yy. Let us assume linear damping, namely damping pro- 
portional to velocities x), y) with coefficient k, (the damping of the foundation mass), and the 
damping of the rotor (induced by the motion of the rotor in the stator and in the bearings) which 
is proportional to the relative velocities « — x) and y — y,, with coefficient k. The order of magni- 


tude of the external damping proportional to x and y is lower than both previously mentioned 
dampings and can, therefore, be neglected. 


For writing the differential equations of the motion let us use the Lagrange equations. Let 


us, therefore, establish the expressions for the kinetic, potential and dissipation energies. 
We obtain for the kinetic energy: 


2 T = m (x2 + y*) + mg (%5 + 9) » (1) 
for the potential energy: 


2 U = cy x8 + Cony + 5 (ea + 9) [ee — m9)? + (yo) 
a5 oi Co) {[(@— Xo)” — (¥ —¥o)"] cos 2 wt + 2 (x — x) (y — yo) sin 2 w tl, (2) 


and for the dissipation energy: 
2W = ky [x + yo] + & [(% —%)? + (¥—Yo)"] - (3) 
Then using Lagrange’s equations: 


*-(3r) OL oU ‘oW ae 


= — = (), j= 1, 2,3,4 4 
dt \0q;) aq; | qj * 0a; (7 ) (4) 
where 
WM=%> F2=¥> 3=%> a= Yo> 


and introducing the small parameter 
Co — Cy 


a 5 
Careenen (5) 


the following differential equations of motion are obtained: 


i + by (@—&) + 9B (em) =p EX, 


2 
i+ 6— Fo) +920 — ye) =e GY, B 
Hy + 840 (Ho —%) + Ooty + AE ay + AE (ay — 2) = —p SEX, 
a + S10 Wo — 9) + S00 + 2% + SH (Yo —¥) = pty, 
where 
j=, bet, =e, 


X = (x— x) cos2mt+(y—yp)sin2ot, 
Y =— (y— yp) cos2wt + (x—%)sin2ot. 


The system of differential equations (6) is a quasiharmonic system where all members con- 
taining the periodically variable coefficients are multiplied by the small parameter. The approxi- 


mate determination of the instability intervals (relative to the angular rotational velocity «) is 
carried out by means of a method described in detail in a former paper?. 


1 4, Tondl, Anlikace matematiky. 4 (1959) p. 278. 
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In the solution by means of the above mentioned method, we assume low damping so that the 
small parameter pt can be factored out. We then transform the system to normal coordinates relative 
to the abbreviated system (for zero value of the small parameter), and introduce substitution 


@ t =T, in our case preferably 2 wt =T. 


Thus we obtain a system of differential equations in the form 


x —Atx, =m S Qaret Pram), (8 =1,2,3,---57) (2) 
k=1 


where 


4atits =H), 


Q, are the natural frequencies of the abbreviated system and Q,, = Q,(t), Psz = P, ,(7) are con- 
tinuous periodic functions of t with period 2 or constants which can be developped into Fourier’s 
series. Let us assume that /, differ from zero and from each other (A, /,, s~ f). 


Substituting 
Me NN a De X,—A, X= Neo (8) 


we obtain the system in the form 
=A tu Fes } G2 ee (9) 
is =—A.n, + BF, 


where 


P= 3 [Oude +m) +7, Pa emo]. 
k= k 


On the basis of Floquet’s theorem, we can write 
es =U, el4 + ivx(u)|r nH, =, el4 + innit , (10) 
where u,, v, are periodic functions 7 with period 2 7 and J is the root of the characteristic equation 


of the abbreviated system. Introducing substitutions (10) into systém (9) we obtain differential 
equations 


= (—2u, +u(—inu, +z Fi), : 
i (Ss 1 5 2,-5 560s y IU) (11) 
yu=—(, +A)», tu[-ine, +75) 


The instability interval can be expected in the vicinity of values! ®® 


Ge aes (2; > 2,) (12a) 
or 
Agee Ap ON ae 5 F he, Sy tek By; (12b) 


where N > 0 is an integer. Then we speak of the instability interval of the N-th order. 
When searching for the instability interval in the vicinity of value 


> Q)4 2, 
ae ce 


? 


we introduce substitution 
Aras we alae (13) 


where « is a real function of M. 


1 I. G. Malkin, Gostechizdat, Moskva 1956, p. 337. 
2 V. A, Jakubovié, Prikladnaja matematika i mechanika, 21 (1957) p. 707. 
° A. Tondl, Strojnoelektrotechnicky éasopis, 9 (1958) p. 384. 


XXIX. Band 1960 A. Tondl: The Stability of Motion of a Rotor with Unsymmetrical Shaft 413 


Substituting A, for A in equation (11) and using substitution (13), we obtain the following system 
of differential equations: 


ue (—inu, +7 F, 


= — iN tel-i@ +a)y + 2A, 
. 1 a2) 
m= Gi — 2) u, + u(— inn, +7 F) (sie ss2-7)", 


vf = — (AF LAP)», +e(-izat oF) (asl, 0 eee). 


Now, let us find a periodic solution of u,, v, period 2 7. The solution is to be in the form of a series 
relative to parameter w carried out by means of an auxiliary system according to Simanov!, 
From the conditions of the solution periodicity, coefficient x is determined with any arbitrary 
accuracy. 

The interval of unstable solutions will be associated with that coefficient x whose imaginary 
component is negative. Coefficient « and consequently the instability interval will be determined 
from this condition. 

According to the procedure referred to above and assuming small values of the damping coeffi- 
cients, the following substitutions are introduced into the system (6): 


0, =U, Oy =H, Oo =UH, 2ot=T. (15) 
For the system under investigation, the general equations (14) are obtained in the form 
Ps ae (A,—A) 9%, eles 


Y=—A,+4)y,+uF << a 


where with k = 1, 3 


mag 
A= 2@’ 
Bu Loren =F ee Cte, +2¢o1 | G1 Tale Cy + Cy + 2 Co1 ure (cy + €2) (er + 2 — 2 1) 
= eee ay ea ie a 2 IM ar 2m 2 mo 2mm { 
a8 lla+e, Cy + cy + 2 Coe Cy + €,\? Cy + eg + 2e90\? , (cr + Cg) (C1 + Cg— 2 Cp 9) 
Oe =F 2m iG 2 Mg £( pascal 2 1M BP 2m my y 


1 K, 
PF, = — 75 19c (Pe + ¥) + Or Pegs + Yesad] + 4g 85 ee 


1 K; 
Fi 41 = — 7 lOvs1 Pasa + Perr) + Pegs Pe + Yel + Tp aa Ae > 


1 
. 0, = ea, [(1 — ay) (4441 91 + Pro) — % Fol » 
O141= ae [4x41 — 1) 91 + P10) + 244190] » 
D,, = as [(l siF a), , 1) (44,44 a = Dro) Ts U1 Po] 2 


1 
D,.4= karo [(a,, — 1) (@ F1 + P10) + % Pol 


1 
Ky = a (tens om a 


iy Rave 


1 §, N. Simanov, Prikladnaja matematika i mechanika, 16 (1952) p. 129. 


414 A. Tondl: The Stability of Motion cf a Rotor with Unsymmetrical Shaft Ingenieur-Archiv 


Cy + Cg 

2 Mo 
rsa Ct Cot 2%1 op 
Tr 


2M 


(r=) 


cy +e 
2m 
9 SO 
J Cy + Co + 2 Coe Q2 
J 


2 Mo 


(j = 3,4), 


(P— wy) +2 a—0)] (C+ 0") 


a4 


—] ‘ t 
MT (G0) (0 *), 


. 3— 1 
— 1/97 @o— vs) + 
3 


A, = ; 7, (3 — Ys) + : ar = (7s— 0) (ef ==) 


a, 


C 1, — ail it —it 
=i [FF rw) +9 v9 | (0 "9. 


When the labile region is to be determined in the vicinity of wy = 2,/N (the labile region of 
the first kind of the N-th order; according to notation used in my former paper, labile region of the 
first kind is the region in the vicinity of w) = 2,/N, while labile region of the second kind is the 
region in the vicinity of @» = (Q, + 2;)/2 N, when Q, > 2;;kAj), then root A, is taken in 
lieu of A and the following substitution introduced: 


24, = iN + pialu). (17) | 


where N= 1,2,3... 

Replacing A by root A, and using substitution (17) as well as adding member W, to the first 
and member W, e—'* to the second equation of system (16), (W,, W, are constants), we obtain 
the following auxiliary system: 


9, =H (Fi, —ixgy) + W,, 
y= —i Ny +4[Fi—i eo +2) yy) + We, 


a 


a | 

w= (5—H)e, +H (Fixes | 9 9 a oF 

Ys =— (CF +At)y. +e (Fi xy,) ae i 

Let us assume that ‘ 
At Lit =Mi (Mis an integer; s=2, 3,4). (19) 

} 


Let it now be required to find periodic solutions for q,, y, (k = 1, 2,3, 4) with period of 27 in 
the form of series 


Pe = Pro +L Per + Pea tess 
Ye = Veo +L Ger + Peg ++. 
Likewise, constants W, and W are to be expressed in the form 
Wi = Wy +e Wy, Fee 
Wy = Woy + Ue Wo, +? Woe +e. 
Selecting the initial conditions 


71(0) = A, ,(0) = B, . (22) 


then the remaining six conditions will depend, according to the Simanov’s theorem, on A and B. 
In harmony with the initial conditions (22) we can select the initial conditions 


P1909) = A, 4 9(0) = B, 


9100) = 9, 0) =0 (G=1,2,3,...).f 23) 
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Introducing into the auxiliary system (18) expressions for q,, y, and W,, W, from equations 
(20) and (21), we obtain for the first approximation 


Pio = A, \ 
De erence 5 
ee tig — 0 (ss 28,4, (24) 


Wi9 =9, Wy =0. 


From the condition of the solution periodicity, the second approximation yields for N = 1 
the following expressions for W,, and W,,: 

1 

Wir aaa ke 


4 Wy 


‘: K —1 
Otindtia eB. 


25) 
Ss a= . K, a, — 1 ( 
War = |p 5,0 +4 +09) Si oy Sea ae 
If the first two approximations are satisfactory, i. e. if 
W,—W,,tuW,,=0, 
1 1o TM My (26) 
W.=W., + UW =, | 


a system of two homogeneous linear algebraic equations is obtained for A and B. From the con- 
dition of a non-trivial solution we can write the following equation for x: 
. (OF Ky (ay — 1) 2 0? Oo 
2 — 1 & — tb —_ | KH — —_ B Sey) xu == 
a (> ios) ( 16 w2 A# } Loman 4a, 7! Co) 
Using the Routh-Hourwitz criterion we get the following condition for the instability of the 
solution, i. e. for the negative imaginary part 


oe ieee 0. (28) 


4-@2 ' 16 w? A* 
Using equation (17) and inserting Aj 1 2,/2 wo, we obtain the appropriate instability interval. 
Using only the coefficients of the first powers of the small parameter jy, the conditions for the 
‘instability intervals are established as 


[Qj+ Q| _ |Q3 & Q,| 
2 


HAn< oO <9 +H Aj (9. = 1,2, 3,4) (29) 


where 


2A 1 0; se ©; Kj K, 
An 4a eS Vt zee —Y u—1)—6; 6. 
The interval limits are real as long as the discrimant under the radical sign is real. Only one 
sign can come into consideration in inequalities (29). Thus a maximum of ten instability intervals 
can be obtained. When higher orders (VN = 2,3...), are also considered, real limits of unstable 
intervals are not obtained for the number of approximations used above. Only additional approxi- 
mations would enable us to also obtain intervals of higher orders as long as the damping remains 
satisfactorily small. 

Let us now analyze the effect of the foundation tuning and the position of the foundation 
natural frequencies relative to the rotor natural frequency. Principally, three cases can be distin- 
guished. In the first case, both natural frequencies of the foundation are higher than the natural 
frequency of the rotor itself in absolutely rigid supports. In the second instance, one of the foun- 
dation natural frequencies is higher and the other lower than the rotor natural frequency. And in 

the third case, both natural frequencies of the foundation are lower than the rotor natural frequency. 

The effect of the foundation tuning is best elucidated by means of examples. Let us assume 
that the mass of the disk and that of the foundation is the same in all three cases, namely m = 0,1 
kg cm? sec?, my = 1 kg cm sec?. Furthermore, let us assume that the average shaft stiffness is 


also the same in all three cases, namely sa = L10e ke em *, 


Let it be required that the stiffnesses of the foundation mounting have the following values 
in the individual cases: 


Co1 [kg cm=}] 4,104 2.1104 6.10 
Cog [kg em=}] 2.104 4.104 2.103 


Ingenieur-Archiv 
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Fig. 2. Boundaries of the instability regions for case I with zero damping. 


Fig. 3. Boundaries of the instability regions for case II with zero damping. 
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Fig. 4, Boundaries of the instability regions for case III with zero damping. 
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The appropriate natural frequencies of the system are 


I. iene II. 
2; 97.8 95.6 | 109.4 
Qs 201.5 147.5 10.7 
Qs 95.6 107.3 105.8 
oe 147.8 58.8 42.8 


Let us select equal damping coefficients for all three cases. In order to evaluate the effect 
of the damping magnitude, let us select four different values of the damping coefficients while 
maintaining their mutual relations: 


k (kg em“ sec) 
ky (kg em=! sec) 0 


Figs. 2, 3 and 4 show the diagrams of the instability regions for all three cases of the foundation 
tuning and for k = 0, and ky = 0. Figs. 5, 6 and 7 show similar diagrams for damping coefficients 
k = 0.4 kg cm™! sec and k, = 6.0 kg cm™! sec. Figs. 8, 9 and 10 contain the same diagrams both 


ee azo \ 020 
O18 G15 015 
G10 0,10 0,10 
OQ) a G05 005 
@o 0 : @ 0 ; @ 
100 10 100 120 700 100 140 
Fig. 5. Case I. : Fig. 6. Case II. Fig. 7. Case III. 


Boundaries of the instability regions for damping coefficients k = 0,4 and k, = 6. 


Yad a 
920 920 a 
O15 015 _| 
: 0,70 010 
; 005 ae 005 
oO 0 @o 
a 700 10 Q 00 120 100. 10 ‘120 
i Fig. 10. C IIL. 
Fig. 8. Case I. Fig. 9. Case II. ig ase 


Boundaries of the instability regions for damping coefficients k = 0,6, ky=9 (light cross-hatching) and k = 0,8, ky = 12 (heavy cross-hatching). 
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for damping with coefficients k = 0.6 kg cm™ sec, ky = 9.0 kg cm™ sec (light cross-hatching) and 
for coefficients k = 0.8 kg cm™! sec, ky = 12 kg cm™! sec (heavy cross-hatching). In the last in- 
stance there exist three merging regions for the case I, two for the case II for the largest damping 
within the given range of parameter yw. For these values of damping, there exists no instability 
region for the third case within the given range of uw. When there is no damping or the damping 
is slight, the case II appears as the most unfavourable in view of the total width of labile regions 
for a given yu. With damping, the coefficient for which labile regions exist no longer attains the 
lowest value ju for the case I, and the highest value for case III which is the most favourable. 


3. Conclusion. Let us first compare the results obtained in the known case of an absolutely 
rigid foundation. For the simplest type of rotor (one disk on a massless shaft) there exists a single 
instability interval of the shaft angular rotational velocity which is limited for the undamped 
motion by the natural frequencies corresponding to the different shaft stiffnesses. This interval 
becomes narrower or disappears completely due to the effect of damping. Thus the instability inter- 
val is spread in the vicinity of the rotor angular velocity corresponding to the average natural fre- 
quency of the rotor [i.e. the frequency corresponding to the average stiffness of the shaft, + (cy + €9)]. 

The case of an elastic mass foundation is more complex. Owing to the foundation elasticity — 
which is practically always different in horizontal and vertical directions — four average natural 
frequencies of the system Q,, Q,, 23, Q, are obtained in the above simplified case. The instability 
regions of the first order which are always the most important are obtained partly in the vicinity 
of these average natural frequencies Q;, (k = 1, 2, 3, 4) — the unstable regions of the first kind — 
partly in the vicinity of combination frequencies 


(Wp => JQ; + 2,| (j= 1,2, 3,4; kA j) OO) 


which can be called the regions of second kind. Thus in our case four instability regions of the 
first kind and six of the second kind are obtained. 

For the cases discussed above all instability regions of the second kind lie in the interval limited 
by the lowest and the highest instability regions of the first kind, i. e. they lie in the vicinity of 
the combination frequencies where only the plus sign is considered (see eq. (29)). 

Thus three substantial instability regions (two of the first kind, one of the second kind) are 
obtained in the vicinity of the natural frequency of the rotor itself, which practically merge into a 
single region. The remaining instability regions are very narrow and completely disappear for not 
too large differences in the shaft stiffnesses. 

For the case I when both partial frequencies of the foundation itself are higher than the average 
frequency of the rotor itself, these three substantial instability regions (merging practically into 
one) for the damped system reach to lower values of yu than in cases II and III. Thus a smaller 
difference in the shaft stiffnesses is necessary for complete disappearance of instability regions for 
the same damping. 

In case II when one partial frequency of the foundation is lower and the other higher than 
the average natural frequency of the rotor itself, these three regions are more extended so that 
for a given which is greater than the minimum y (for which labile regions disappear owing to 
damping) the total width of regions will be greater than in cases I and III. 

In case III when both partial frequencies of the foundation are lower than the average natural 
frequency of the rotor itself, this case is more favourable for a damped system than cases I and II, 
particularly at high damping of the foundation mass. 

The following holds true for equal damping in all three cases: 

If, for a given value of the difference in the shaft stiffnesses, the rotor motion is stable in case I 
within the entire range of rotational speeds, the motion is always stable in case II and case III 
for the same rotor damping; the opposite, however, does not apply. 

This leads to a practical conclusion, namely were the motion of the two-pole turbogenerator 
stable when mounted on a high-tuned foundation, it would be much more stable on a low- 
tuned foundation, as long as the values of damping are not reduced. Since the values of damping 
of the rotor relative motion are identical for the same machine, this means that the damping of 
a low-tuned foundation must not be lower than that of a higher-tuned foundation. 


(Eingegangen am 18. Januar 1960) 
Anschrift des Verfassers: Dr. Ale Tondl Praha 1 (CSR) Zborovské 41. 
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Die Spannungen in der nur am Rande belasteten, 
um einen Durchmesser verbogenen Kreisringplatte 


Von Jiirgen Raasch 


1. Einleitung. In der Verfahrenstechnik werden auf Hohlzapfen gelagerte rotierende zylindrische 
GefaBe verhaltnismaBig haufig verwendet (Kugelmiihlen, Mischtrommeln, F iltertrommeln). Ihre 
Stirnplatten werden, wie Abb. 1 erkennen laBt, durch die Auflagerkraft A auf Zug und Druck und 
j durch ein Moment M, = A ss auf Biegung beansprucht. Im allgemeinen diirfen die von der Auf- 
) lagerkraft verursachten Zug- und Druckspannungen gegeniiber den Biege- 
spannungen vernachlassigt werden. Im folgenden soll daher — als ein Problem 
von erheblichem praktischem Interesse — die Beanspruchung einer Kreis- 
ringplatte durch Biegemomente am Innen- und AuSenrand untersucht 
werden, wobei zunachst eine vollkommen starre Einspannung an den Randern 
angenommen wird. 

Ein ahnliches Problem ist im Jahre 1929 von H. Reifner! behandelt 
worden. Er untersuchte die Spannungen in einer flexiblen Kupplung in Ab- 
hangigkeit vom Verdrehungswinkel. Die Ergebnisse lieBen sich nicht in ge- 
) schlossenen Ausdriicken darstellen und wurden deshalb in Form einer kleinen 
_ Tabelle mitgeteilt. 

Demgegeniiber ist es in der vorliegenden Arbeit gelungen, eine vollstandige 
Lésung in geschlossener Form zu erhalten. Fiir den Héchstwert der Span- 
nungen wird ein einfacher Ausdruck gefunden, der als Berechnungsgrund- 

_ lage fiir den konstruierenden Ingenieur von Interesse sein diirfte. Weiterhin 
| wird hier untersucht, wie sich der Spannungszustand in einer Kreisring- 
platte andert, wenn die Bedingungen einer starren Einspannung am Innen- 
und AuBenrand nicht erfillt sind. 

Die Grundlage der mitgeteilten vollstandigen Lésungen wie auch der Lésung von H. Reiner 
ist die Kirchhoffsche Plattentheorie mit ihren bekannten Naherungsvoraussetzungen. 


ud 


SSSSSSSSSSESSS ENNIS] 


Abb. 1. Kugelmihlen- 
Stirnwand, 


2. Ausgangsgleichungen. Fir die nur am Rande belastete Platte lautet die Grundgleichung 
der Durchbiegung wy 
AAwy = 0 2 (1) 
Bei Verwendung von Zylinderkoordinaten r, py hat das Symbol A die Bedeutung 
a Uo il ae 
=75, +-—-=s+3c 2 

eae ay ae er re (2) 

Mathematisch gesehen besteht die Aufgabe darin, eine Lisung von (1) zu finden, die den jeweiligen 

_ Randbedingungen geniigt. Ist die Durchbiegung wy als Funktion der Koordinaten r und g bekannt, 
so kénnen die Spannungen fiir jeden Punkt der Platte berechnet werden. Es gilt 


26 OW 1 dw 1 d*w 
One Fy = z|m Bee cord art ag? |’ 
2G OW 1 dw, 1 wo 3 
= pet | pas a 8 ag? }|’ (3) 
* 0/1 aw 
Trg = — 2625, (5 Gp): 


Fiir die in radialen und konzentrischen Schnittflachen je Langeneinheit ibertragenen Krafte und 
| Momente erhalt man damit 


N Ow, 1 dw 1 dw era as oe 1 &wo 
eee Lacuna a Lee or cas ay? |’ Bee on taal Orhan pare OF ees ag? ||? 
N @ /1 éu agi. Aw, Were ny ee 
M9 =—(™— Va oe ( }? <a Sete a op 


14, Reiner, Ing.-Arch, 1 (1930), S. 72. 
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Die Spannungsresultanten wurden angegeben, weil sie im folgenden zur Formulierung der Rand- 
bedingungen verwendet werden. Abb. 2 gibt itber Richtung und Vorzeichen Auskunft. 
Mit 
mG h® 
N= 6(m—1) 2 
ist dabei die sogenannte Plattensteifigkeit eingefiihrt. Die Buchstaben G und m bezeichnen wie 
iiblich, Schubmodul und Poissonsche Zahl (Stahl m ~& 10/3), und h bedeutet die Dicke der Platte. 


Abb. 2, Spannungsresultanten an einem Plattenelement. Abb. 3, Koordinatensystem. 


3. Lésung fiir die Kreisringplatte mit starrer Einspannung am Innen- und Aufenrand. Die Lage 
des Koordinatensystems ist aus Abb. 3 zu entnehmen. Mit 6 und a werden Innen- und AuSenradius 
der Kreisringplatte bezeichnet. Es wird angenommen, dai die Platte durch ein Moment vom Be- 
trag M, um die x-Achse verbogen wird. 

Wenn fiir die Lésungsfunktion der Ansatz gemacht wird 


wy = (47° + Br ++ Drinr) cos, (6) 


so entspricht dies allen Forderungen der Symmetrie. Man kann leicht zeigen, daB (6) die Grund- 


gleichung (1) erfiillt. 
Zunachst sollen an Hand der Abb. 2 die Gleichgewichtshedingungen fiir den Innen- und Aufen- 
rand formuliert werden. Man erhalt die Gleichungen 


22 
ferdg=0, 
f 
2a 
al GTR OSAP lett Nee ROe P) dp — Mg 5 (7) 


27 
Ft Teue Goi, 5608 —- mabin g) 7 dip =). 


Mit dem Ansatz (6) erhalt man die Spannungsresultanten aus (4) zu 


=—N(84A—“yJeosp, | 
(ere. BG D eae 0 
m =— (m6 Ar+ e+) t24r— “Ft Ploeg, (8) 
Meg =m —1) 5 (2Ar— Ay + P)sin g | 


Setzt man diese Ausdriicke in die Gleichungen (7) ein, so erkennt man, daB zur Befriedigung der 
Gleichgewichtsbedingungen die Festlegung einer einzigen Konstanten geniigt. Man findet 


te MG 
AGN - (9) 


| 
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Weitere drei Randbedingungen ergeben sich auf Grund der Vorstellung, daB Innen- und AuBen- 


rand der Platte eben blei i 
See in eiben, wobei der Innenrand als Ganzes gegen den AuBenrand um die x-Achse 


1) KOs Ne 0°, 
2) Pasa = = 0, | 
i l (10) 
3) i b s Ow = bw cos pene dw sin = Wo 
dy ér " rp Lanes rcosp 
Fur die Konstanten erhalt man so 
M, 1 M py e: 
iis ages 2 Me iat =e M, ab? 
82a N a2 + b?’ - 8a N (Gea 2ma), C= aaN@LE (Hy) 
Damit lautet die vollstandige Lésung 
See Ge ie) re 8") r 
t= Fh | r@+P) +2rln 7] cos. (12) 
Die Spannungen folgen durch Einsetzen in (3) zu 
Pe 3 My [{[r*b(3m-+ 1)—r*b(@ + B’) (m+ 1) — a? b? (m — 1) 
as ma h?b Sina ieape cos, 
2 mr (a? + 6?) 
ge 3 Mo [ tb (m + 3)— 7b (a? + BY) (m+ 1) + Pb (m—}) 
0 a Wb 2 mr (a? + 6%) ae eS) 
3 My[m—1 b(@—Pr)(27—b%)].. 
Tro ib | 2m r° (a? + 6?) san 


Dabei ist z = h/2 gesetzt; d.h. die Gleichungen (13) geben die Spannungen an der Plattenober- 
flache an. 


Spannungsverlauf 
bei Einspannung am Innen-und AuBenrand. 


Abb. 4, Spannungsverlauf bei Einspannung am Innen- und Aufenrand. 


Abb. 4 zeigt den Verlauf der Spannungen fiir ein Verhaltnis von AuBen- zu Innenhalbmesser 
a/b =3. Der Klammerausdruck ist dabei jeweils tiber dem Radienverhiltnis r/b aufgetragen, fiir 
die Poissonsche Zah] wurde m = 10/3 eingesetzt. Der Kurvenverlauf ist auch fiir andere Werte des 
29 


SEA 
: 


“A = 

ts 

ae tar 
n 
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Verhaltnisses a/b charakteristisch. Man erkennt, da8 der Héchstwert der Spannungen durch die 
Spannung o, am Innenrand gegeben ist: 


3 M, a? — Bb? 14 
Cmax = |Grlr=b a h? b a? + Bb” cS 


Als auffallig ist zu vermerken, da die Poissonsche Zahl in (14) nicht auftritt. Am gleichen Ort 
erhalt man fiir die Spannung o, 


eee (15) 


Aus beiden ist evtl. eine Vergleichsspannung zu bilden. Geht man von der Mohrschen Schubspan- — }j 


nungstheorie aus, so ist (14) die fiir die Berechnung einer Kreisringplatte maSgebende Beziehung. 
Abb. 5 macht deutlich, wie die Beanspruchung einer Kreisringplatte vom Verhaltnis a/b abhangt. 


a*-b* 


a’+b? 


70 


08 


06 


04 


02 


7 15 2 3 ¥ 5 
Abb. 5, Abhangigkeit der Beanspruchung vom Verhiltnis a/b. 


Der Betrag des Winkelsa, um den Innen- und Au®enrand gegeneinander verdreht werden, ergibt 
sich aus (12) aufgrund des Zusammenhanges 


= Mechel Rett i 
sa ~ ay reosy ce, 
zu 
Meee b 
x inn (eee ts =a , (17) 


4, Lésungen bei veranderten Randbedingungen. a) Beide Rander nicht eingespannt. Im 


praktisch gegebenen Fall ist die Einspannung der Kreisringplatte niemals vollkommen starr, son- 17 
dern mu als mehr oder weniger elastisch angenommen werden. In diesem Zusammenhang inter- | 
essiert daher die Frage nach dem Spannungszustand, der sich einstellt, wenn die Bedingung der | 


Einspannung ganzlich fallengelassen wird. 


Man kommt auch in diesem Fall zu einer Lésung, wenn man bei gleichem Ansatz die in den 
Gleichungen (10) formulierten Randbedingungen ersetzt durch 


1) Ps Wy =0, 
2) tee i, == 0), (18) 
3) r= om, = (., 
Die Konstanten andern jetzt ihren Wert. Man findet 
ees M, m-+1 1 
82N 3m+1 a+ Bb’ 
a MG m+] a m-+ 1 b2 ; 
ee eet porere 2a), (19) 
‘gg M, m+1 ab 


> 8a N m—1 ote’ 


sit edd pid i tle OA Re El Ln 


" iG saad 
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wahrend D erhalten bleibt. Die vollstandige Lésung fiir die Durchbiegungen lautet in diesem Fall 


ata M, jicrac cae Cateeraad Mandel eo i ae r = 
Oo Aa See ete be, m= 1 (@—B)r 2rln ako (20) 
Fir die Spannungen erhalt man damit 
Gee Sess Wisdom b (r? — a®) (72 — b2) 
r ay Oe ear oa 73 (a? + B2) cos@, 
og. = 2 Mo m-+1 f/m+3 rb a? 63 b\ 
ea Wb] 2m \3m+lae@+e (a+ 8%) _ ol (21) 
pee 3 M, [m—1 m+1 a? 68 m+ 1 rb ae ; 
Tee heb 2a \ apie P(@+e) 3mtle@+ep ! sing. 
Thr Verlauf ist aus Abb. 6 zu erkennen. 
XK 
hoe 


o8 


66 


Spannungsver/auf 
|_ohne Einspannung am Innen-und AuBenrand _| 
a/b=3 


Abb. 6. Spannungsverlauf ohne Einspannung am Innen- und AuBenrand. 


Kine gewisse Schwierigkeit ergibt sich insofern, als die Schubspannungen am Innen- und Aufen- 
rand nicht zu Null werden. Dieses Problem besteht indessen immer, solange man von den Ansatzen 
der elementaren Plattentheorie ausgeht. Es gelingt nicht, am freien Plattenrand auch das Scher- 
moment zum Verschwinden zu bringen. Man trifft daher die Annahme!, daB das Schermoment 
am Rande durch eine geeignete Verteilung der Querkrafte ersetzt werden darf. Diese von Kirchhoff 
zuerst formulierte Beziehung lautet 


(a eels (22) 


Man setzt dabei voraus, da sich die damit verbundene Stérung des Spannungszustandes entspre- 
chend dem Prinzip von de St. Venant auf eine schmale Randzone beschrankt. Weiterhin ergibt 


sich aus Abb. 6, daB als fiir die Berechnung einer nicht eingespannten Kreisringplatte maBgebend 


die Schubspannung Tg in der Nahe des Innenrandes anzusehen ist. Man erhalt die Beziehung 
3 M, a@(3m-+ 1)+ b?(m—1) ; (23) 
Tmax hb Gm+l)@+e) 


Mit (23) ist die Beanspruchung einer Kreisringplatte im ungiinstigsten Fall gegeben. 


1A, Nadai, Elastische Platten, S. 34. Berlin 1925. 
29* 
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Bei einem Vergleich mit (14) ist zu beachten, das die maximalen Spannungen nicht an den glei- 
chen Stellen des Innenrandes auftreten (Verlauf von Sinus- und Cosinusfunktion). Daher bedeutet 
eine in geringem MaBe nachgiebige Einspannung noch keine Erhéhung der durch (14) gegebenen 
Beanspruchung. 

b) Nur ein Rand eingespannt. Durch Variation der Randbedingungen in (10) und (18) 
erhalt man weitere Lisungen fiir eine Kreisringplatte, die entweder nur am Innenrand oder nur 
am Aufenrand eingespannt ist. Der Verlauf der Spannungen ist fiir diese Falle in Abb. 7 und 8 
dargestellt. 


2b a h?h 
to 3 ie te 3 My 
40 10 | 
48 08 
06 06 T =: 
OY 04 
en ee %(¢-9) 
Jno (0-2) \ 
02 
o9 75 2 3 
02 
{ 
Op (9-0) 
-0,4 
-06 | 
Or (Y=), , 
(9-9) Spannungsver/auf Spannungsverlauf 
-98 zal bei Einspannung.am Innenrand -08 | bei Einspannung am AuBenrand. 
afb=3 alb=3 
-40 -10 > 
Abb. 7. Spannungsverlauf bei Einspannung am Innenrand. Abb. 8. Spannungsverlauf bei Einspannung am Aufenrand. 


Herrn Professor Dr.-Ing. H. Rumpf michte ich an dieser Stelle fiir die Anregung und die 
Unterstiitzung der Arbeit herzlich danken, ebenso Herrn Professor Dr.-Ing. R. Sonntag fir 
wichtige Hinweise. 


Mitteilung aus dem Institut fiir Mechanische Verfahrenstechnik an der Technischen Hochschule Karlsruhe. 


(Eingegangen am 18. Januar 1960.) 


Anschrift des Verfassers: Dipl.-Ing. Jiirgen Raasch, Karlsruhe, Technische Hochschule, 
Institut fiir Mechanische Verfahrenstechnik. 
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Uber die Stabilitat der Prazessionshewegung eines symmetrischen Kreisels 


Von O. Bottema 


1. Wir betrachten einen symmetrischen Kreisel, der sich um den auf der Symmetrieachse | 
gelegenen festen Punkt O drehen kann. Das Tragheitsmoment beziiglich I sei gleich C und fiir eine 
Achse durch O senkrecht auf / gleich A. Die Lage des Kreisels werde durch die drei Eulerschen 
Winkel ?, y und y in bezug auf ein rechtwinkliges raumfestes Koordinatensystem OX YZ dar- 
gestellt; dabeiist der Winkel zwischen I und OZ und y derjenige zwischen der Ebene (I, OZ) und 
der X-Achse; ist die Eigendrehgeschwindigkeit um I. Fiir die kinetische Energie des Kreisels 
findet man 


1 2 1 . 5 52s 1 - 
i ee 0 - (A sin? # + C cos? #) y? + Cyq@eos 3 + oF Co. (1) 
Beziiglich der Krafte, welche auf den Kreisel wirken, nehmen wir an, daf dieser eine potentielle 


Energie V besitzt, welche nur von @ abhangt. Dann hat das System die zwei zyklischen Koordina- 
ten y und ¢, und es gelten die Integrale 


(A sin? @ + C cos? 3) py + Cpcos 3 =c,, (2) 
Cycos8 + Cy =e, (3) 
wo Cy und c, Integrationskonstanten sind. Fiir die dritte Bewegungsgleichung findet man 
(4) —(A—C)#sin I cos F + Cppsind + =0. ae 
oder, wenn man yp und @ vermittels (2) und (3) eliminiert, 
a : 


Der Kreisel fiihrt eine Prazessionsbewegung aus wenn @ konstant bleibt. Ist 0 =a (a + 0, 
« + 7) und dV/d? fiir 0 = « gleich p, dann hat man 


(c, — ¢, cos #) (cg — c, cos 9) + A psin? ? = 0. (6) 


Der Ausdruck (c, — cy cos #) (cy — ¢, Cos #) ist eine nicht-definite homogene quadratische Funktion 
von ¢, und c,; er kann also alle positiven und negativen Werte annehmen. Das heift: es gibt immer 
(undenklich viele) Zahlen c, und cy, welche bei vorgegebenen « und p die Gleichung(6) befriedigen. 
Oder: der Kreisel kann fiir jeden Wert ® =« und in jedem konservativen Kraftfeld V() Pra- 
zessionsbewegungen ausfiihren. 
Aus (2) und (3) geht hervor, daf fiir eine solche Bewegung und @ konstant sind. Man hat 

dabei 

° Ci — Cy COS & 

baci a sin? a (7) 


Wenn wir noch den Wert w + 0 dieser Prazessionsgeschwindigkeit vorschreiben, dann ist 


; ¢,— ¢,c0s 0 = Awsin®a, 
wahrend aus (6) folgt 


ela B De 
Cy — €, COSA = ee 
Man hat also 
are (8) 
oe igeeeearall| w? sin % — p cosa), 
ae aR (A @* sin & cos o — p) 2 2 (°) 


so daB c, und c, eindeutig bestimmt sind. Das heifit: bei vorgegebener V, « und w gibt es immer eine 
wohlbestimmte Prazessionsbewegung. 
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2. Wenn V, a und w gegeben sind, ist auch p bestimmt. Aus (3) und (9) folgt 
__ (A— C) o* sin « cos « — p 10 
? Cwsina : 0} 
Dies ist eine gebrochene quadratische Funktion von w; ihre Gestalt wird durch das Vorzeichen von 
k =(A —C) psina cosa bestimmt. Fir k > 0 nimmt ¢ bei veranderlichen w jeden Wert von 
—oo bis +00 an. Ist jedoch k < 0, dann hat g extreme Werte fiir ~? = — p*/k und es ergibt sich 
die Ungleichung 


—Ak 
; eee Th 
P > Bs (11) 


n?@” 
3. Wir wollen jetzt die Stabilitat der Prazessionsbewegung untersuchen. Die Bewegungs- 
gleichung (5) hat die Gestalt : 

A+f(9) =0, (12) 
wobei jetzt wegen (6) f(a) = Oist. Schreibt man ? =a + ¢, f(?) =f(a) +ef'(a) +:--, dann gilt 
in der Nahe der Gleichgewichtslage ¢ = 0 
Aée+fi(a)-e=0, 

d. h. man hat Stabilitat falls f(a) > 0 ist. 
Die Prazessionsbewegung ist also stabil wenn 
sin? « (c? — 2 c, c, cos a + c2) — 3 cos a (c, — cg cos X) (Cg —¢, cosa) + Aqsinta>0 (13) 


ist, wo q den Wert d?V/di* fiir 9 = « bedeuten soll, oder wenn wir c, und cy mittels (6), (8) und (9) 
eliminieren, 


A? «4 sin? « + Aw? sin «(p cosa + qsina) + p?>0. (14) 
Fiihrt man die neue Variable u = cos ? ein, so daB V eine Funktion von u wird, dann ist 
; dV = dV @V 


i a dV 
ee oink AEs Pee. ae 
7B qgz = sin oO hie cos } sie 


Bezeichnet man dV/du und d?V/du? fiir 9 = « mit P und Q dann lautet die Bedingung (14) noch 
etwas einfacher 


=— sind. 


F= 4! + Aw?(—2 Pcosa +Qsin?«) + P?>0. (15) 
4, Die Diskriminante der in w? quadratischen Funktion F ist 
D =— A’ sin? « (P + Q cos? 3) [P — Q sin? 2) 3 (16) 


F ist positiv fiir jeden Wert von «, also fiir jeden positiven Wert von w?, wenn F keine positiven 
Nullstellen hat. Dafiir gibt es zwei, sich iibrigens nicht ausschlieBende Méglichkeiten: 


2. P cos «+ Osin* « > 0 ee ee 


ad 
(P + Q cos? +) (P—Osin 3) +0. (17’) 


Im ersten Fall sind etwaige Nullstellen negativ, im zweiten Fall gibt es gar keine Nullstellen. 
Man hat immer 


—l1 2 cos & i (18) 
“i 3 
a : 
cos? —_ tes sin? * 
2; ia 
Bezeichnen wir kurz 
1 —] 2 cos « 
ee Oe a == bas sin? @ eee 
sin? oy cos? — 


so ist : 
be 0 bee. bias = be 
In der (P, Q)-Ebene ziehen wir die Geraden L; mit den Gleichungen Q = b; P (Abb. 1). In zwei 
_von den durch L, und L, gebildeten Winkeln liegt die P-Achse, und innerhalb dieser Winkel ist 
offenbar die Ungleichung (17’) erfiillt. Innerhalb dieser zwei Winkel liegt auch die Gerade Lz. Die 
Ungleichung (17) wird erfiillt von den Punkten, die oberhalb von L,liegen. Das Gebiet, wo wenig- 
stens eine der Ungleichungen (17) und (17’) gilt, ist in Abb. 1 schraffiert. 


ie Tee 


Pret plied BE 


Op es 


so muB 


sein. 
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Wir haben also folgendes Ergebnis: Wenn der Kreisel fiir 9 = a (a == 0,« =: 2) eine Prazessions- 
bewegung ausfithrt, V(u) die potentielle Energie darstellt, u = cos 0, dV(a)/du = P und d?V(a)/du? 
= Q ist, dann ist die Bewegung fiir jede Prazessionsgeschwindigkeit w stabil, falls Q > 0 und falls 
Q=0, P+ 0 und falls g<o, 2 


s iB 
zs 0 > sin? > und-falls 0 = 0, = < + cos? — ist. Dagegen ist fiir 
eer ee 
Q< 0, — cos? rs <— =< gin c die Bewegung nur dann stabil falls «2 < wij, oder w? > w? ist, wo 


2 


wo; und «3 die (positiven) Nullstellen von F sind. 


Abb. 1. 


5. Bewegt sich der Kreisel unter Einflu8 der Schwere, ist m seine Masse, S sein Schwerpunkt, 
OS =h und OZ vertikal nach oben gerichtet, dann ist V = mghu, also P=mgh, Q = 0, und 


es geht der bekannte Satz? hervor: Die Prazessionsbewegung ist immer stabil. 


Man hat weiterhin k = — mg h(A — C) sin? a cos a, und (11) sagt aus: Wenn (A — C) cosa 
< 0 ist, dann kann ¢ jeden Wert annehmen; ist dagegen (A — C) cos « > 0, dann ist bei der Pra- 
zessionsbewegung immer . 

Gee EE OES (19) 


Als zweites Beispiel betrachten wir das Feld, wo die Kraft senkrecht auf der Z-Achse steht 


2 : : aa: 
und mit dem Abstand r des Aufpunktes von der Z-Achse proportional ist. Das Potential ist = ar, 
wo a eine Konstante ist, und man hat also 


Va5a(A sin? } + C cos? #) 
und damit ‘ : 
a =—o(A—C)u, P=—al(A —C)eosa, Q=—a(A—C). 
Aus den Bedingungen geht nun hervor: Fiir a(A — C) > 0 ist die Bewegung fiir jedes w nur dann 
stabil, falls cos a >> oder cosa < 5 wahrend fiir a(4 —C)< 0 die Prazession immer 
fiir jedes q stabil ist. 

Fiir dieses Beispiel ist k = (A — C) psina cos 4 = a(A — C)? sin? 4 cos*a; das heiBt wegen 
(11): Falls a > 0 ist, unterliegt die Rotationsgeschwindigkeit y keiner Bedingung; falls a < 0 ist, 
} ; —4a(A —C)* cos? a 

gy > @ 


: (Eingegangen am 28. Januar 1960.) . 
Anschrift des Verfassers: Professor Dr. 0. Bottema, Delft (Holland), Julianalaan 132. 


. 1 R, Grammel, Der Kreisel, 2. Aufl. Bd. 1, S. 96—97 und S. 240—243. Berlin/Géttingen/Heidelberg 1950. 
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Die Inversion der Torsionsspannungen prismatischer Stabe 
in allgemeinen Koordinaten* 


Von H.-J. Franeck 


1. Einleitung. Allen elastizitatstheoretischen Untersuchungen prismatischer Kérper, die einer 
Verdrehung unterworfen sind, liegen die de St. Venantschen Torsionsgleichungen fiir ein orthogonales 
kartesisches Koordinatensystem zugrunde!. In der folgenden Arbeit werden die Differential- 
beziehungen der Torsionstheorie prismatischer Kérper fiir allgemeine Koordinaten abgeleitet, wo- 
bei von den Grundgleichungen der Elastizitatstheorie ausgegangen wird’. 

Fiir spezielle Koordinatensysteme hat die Torsionsdifferentialgleichung in der transformierten 
Form dasselbe Aussehen wie in kartesischen Koordinaten. In diesem Falle gelingt es, aus bekannten 
kartesischen Lésungen geschlossene Lésungen fiir Querschnitte zu berechnen. Die Tatsache der 
Invarianz der Differentialgleichung war fiir eine bestimmte Koordinatenbeziehung von F. Levi- 
Civita? und J. H. Michell4 fiir die Bipotentialgleichung entdeckt und spater von A. Timpe® und 
W. Olszak® fiir ebene Spannungsprobleme angewendet worden. Dabei machte es sich nachteilig be- 
merkbar, daB die dabei verwendete Inversion bekannter Lésungen in einer Abbildung innerhalb 
des kartesischen Koordinatensystems besteht und die sich ergebenden Resultate unter Umstanden 
uniibersichtlich werden. Fat man dagegen die Inversion als eine Transformation in spezielle 
krummlinige Koordinaten auf, so vereinfachen sich die Ergebnisse. 


2. Geometrische Grundlagen. Die rechtwinklig kartesischen Koordinaten eines Raumpunktes 
sind x‘, xJ, x*, x7, seine krummlinigen Koordinaten x’, x“, x”, x7, wobei die Indizes 1, j, k, q bzw. A, u, 
y, o jeweils an Stelle der Zahlen 1, 2, 3 stehen. Die kartesischen Koordinaten %' seien differentiier- 
bare Funktionen der allgemeinen Koordinaten x’ und umgekehrt. Durch partielle Differentiation 


axa 


BU) (irs 20 = () 


erhalt man die Tangenten- bzw. Gradientenvektoren 
Cy = x45 cl = xi. 
Thre nochmaligen Ableitungen werden mit 
5 Sens Yigal a 
Chu = lu > Cin = Glu 


bezeichnet. Setzt man voraus, daB iiber zwei gleiche in einem Produkt vorkommende Indizes zu 
summieren ist, dann ergeben sich die Komponenten des Maftensors in kovarianter bzw. kontra- 
varianter Darstellung zu ; 


Sore: ie 
Bin = 4,6,» ee Tet 

Die kovariante Ableitung eines Vektors wird durch einen langen senkrechten Strich gekennzeichnet 

und mit dem Christoffelschen Dreizeigersymbol 


1 
7 heen ae aie (Siol» oe SAr\o — Bav)a) 
aus 


MS At A’ ry, 
berechnet. 


* Auszug aus der gleichnamigen von der Bergakademie Freiberg genehmigten Di ion ; 
Exot. Dr.-Ing. habil. D. Riidiger; Korreferent Pool ae ele: habil. A. ‘iethte : ee pe ae 
; Vel. z. B. C. Weber und W. Giinther, Torsionstheorie, 1. Aufl., Berlin 1958. 
; vel Le o E. ro Deuker, sou ne Mathematik 5 (1940/41), S. 94. 
- Levt-Civita, Sopra una trasformazione in sé stessa della equazione = i 
‘ J. H. Michell, Proc. Lond. M. Soc. 34 (1901/02), $. 134. fae Senet 
> A. Timpe, Z. Math. Phys. 52 (1905), S. 348. 
6 W. Olszak, Ing.-Arch. 6 (1935), S. 402. 
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Die technische Anwendung der Ergebnisse elastizitatstheoretischer Untersuchungen macht die 


Kenntnis der physikalischen Koordinaten yon Tensoren erster und zweiter Stufe erforderlich. Man 
erhalt fiir einen Vektor 


A(A) = V/g,, A* (nicht summieren), 
wahrend sich die physikalischen Koordinaten eines Tensors zweiter Stufe aus der Relation 
A(A u) = ye A‘u (nicht summieren) 


ergeben. 


Fir die Darstellung eines prismatischen Kérpers mit beliebigem Querschnitt ist es zweckmabig, 
die x?-Linie der allgemeinen Koordinaten mit der x°-Achse eines kartesischen Koordinatensystems 
und der Stabachse zusammenfallen zu lassen. Dann wird der Stabquerschnitt x? = x° — konst. 
durch die ebenen krummlinigen Koordinaten «1, x? beschrieben, und man erhilt mit x! = x (21, 2°) 
x’ = «3 und a, B, y, 6 = 1,2 die Beziehungen 


? 


B= bes=9, Ss 1, Mags=ITa,=0, Ade, = Alm, 


3. Differentialgleichung der Verwélbungsfunktion. Die fiinfzehn Gleichungen zur Berechnung 
des symmetrischen Spannungstensors S’*“, des symmetrischen Deformationstensors D’“ und des 
Verschiebungsvektors V* kénnen bei nicht vorhandener Massenkraft in der Form 


Sil =0, (En) 
Div == (Vale 4 ver), (1.2) 
Sau =26(D +28 Ds) (1.3) 


angeschrieben werden. Aus ihnen erhalt man die Grundgleichungen der Elastizitatstheorie 


1 
1—2yp 


Veale ae yer = 0 (2) 


_ zur Bestimmung der drei Komponenten des Verschiebungsvektors. 


e \\ 


=a 


Fur das vorliegende Torsionsproblem erfolgt die Lésung dieses simultanen Differentialglei- 
chungssystems mit dem Ansatz 


V4 = & [64 @ (x1, x?) — 8,9 8? cf x' x3]. (3) 
Hierin sind # eine vorerst beliebige Konstante, @ (x1, x*) eine noch zu bestimmende Funktion und 
€,4g der bekannte ebene e-Tensor. ae 
Mit dem Ansatz (3) erhalt man unter Beriicksichtigung des Satzes von Ricci 
ak are Binly =0 
fiir die kubische Dehnung 


Damit vereinfacht sich (2) zu 
Vee 0, 


und es wird 


8 02 @ (x1, x?)|# =(). (4) 


Die Gleichungen (4) stellen drei Differentialgleichungen dar, die fiir 1 = a identisch erfullt sind, 
wahrend sich fiir 1 = 3 die Beziehung 


p(x, x*)|4 = Ap(x*, x7) = 0 (5) 


-ergibt. Somit ist die Lésung der elastischen Grundgleichungen (2) auf die Integration der Potential- 


leichung (5) zuriickgefiihrt worden. apts : 
: Um die Bedeutung der Konstanten # zu untersuchen, wird die Winkelanderung D(12) ermittelt. 
Man erhalt 


\/ , ie Lv) gaa ; Bee 
D(12) = St pe — 2 |/ vip yap) = 21/8 (gt 4], +2 V2|,)=0, 
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und dieses Ergebnis gestattet es, die Formanderung eines Flachenteilchens in z. B. orthogonalen | 
kartesischen Koordinaten nach Abb. 1 darzustellen. Hieraus liest man fiir kleine Verschiebungen | 


sind =a = Von» 
sind = a4 = — Vij2 
ab und kann damit die Beziehungen = < 
= Vij + Var» 
opie be ae 
a(x}, x", eh ————— oy €3 ij Vii; 


bilden. Transformiert man diese Gleichungen in allgemeine Koordinaten, so entspricht die erste der 
verschwindenden Winkelanderung D(12) = 0, wahrend sich fiir die zweite 


Stal Sik MPA Bt Mao a 


a(x, x2, x3) = — es Ep V*|6 (6) 


Abb. 1, Formanderung eines rechtwinkligen Flachenteilchens. Abb, 2. Randbedingung fiir den Spannungsvektor am Querschnittsrand, 


ergibt. Bei Beachtung der Rechenregel 
E*P ey5 = 6% 68 — 63 68 
folgt mit dem Ansatz (3) aus (6) 


a(x}, x2, x8) = a(x3) = 0x3. 


Die zunachst willkiirliche Konstante @ ist also mit der Verwindung des Stabes 


a 
0 == as 
identisch. Fiir die Komponenten des Spannungstensors S*“ erhalt man 
Sl1 — 522 — g12 — 933 — 9, (7.1) 
$34 — GP [p(x1, a) — 6B c, rAr ‘ (7.2) 


Als Randbedingung muB fiir die Lésungen von (5) gefordert werden, daf die resultierende Span- 
nung aus den physikalischen Komponenten des Spannungsvektors S** die Berandung des unter- 
suchten Stabquerschnittes tangiert (Abb. 2). 

Die Randbedingung lautet 

dx(1)_ S(31) 
dx(2) ~ 8(32) 
und kann in der Form 


a 1 PaaS 
: ap Plot, 22)|* dx? = 2d, B) | (8) 
angegeben werden. 


_ 4. Differentialgleichung der Spannungsfunktion. In kartesischen Koordinaten laGt sich jede 
Funktion 9(x',x*), die der Laplaceschen Differentialgleichung Ag(x!, x2) = 0 geniigt, durch die 
‘Cauchy-Riemannschen Differentialgleichungen 


og(z!, 32) dy, 22) op(e,3*) y(t, 32) 
Ce Raped yf Oe Oe TE Ne eee 
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“ine zweite Funktion p(%1,%?) so zuordnen, da diese ebenfalls die Laplacesche Differentialgleichung 
rfillt. Die Funktionen g(x, x2) und (x1, x2) sind Skalar- bzw. Vektorpotential eines ebenen 
Laplaceschen Feldes, so daB man fiir dessen Feldvektor mit eis = — &, = 1 

A, = 921, ¥*) i = €1j p(2", X*)); (9) 
schreiben kann. Der Vektor A; ist nicht-normierter Gradientenvektor der Kurven der Kurven- 


char ¢(x1, x") = konst. und gleichzeitig nicht-normierter Tangentenvektor der Kurven der Kurven- 
schar w(x!, x”) = konst. 


Die Transformation von (9) in allgemeine Koordinaten liefert fiir die kovarianten bzw. kontra- 
varianten Komponenten des Feldvektors 


A, = 9(x", x*)|, = €4 8°” p(x}, 2%)|,, (10.1) 

A® = o(sxt, x2)|* = e*F y (x1, lee (10.2) 
Als Randbedingung fiir die Funktion y(x1, «?) erhalt man aus (8) 
1 


Exp E*" (xt, x?) Es dx? — dy (x1, x?) = a d(x! x:) ; 


p(x, x?) — B(x, A) x(x ox°) == komt, (11) 


Zur weiteren Berechnungsvereinfachung ist es vorteilhaft, eine Funktion F(x, x2) gemaB 
F (2, 22) = 6 |y(o, 22) 22s, (12) 


sinzufithren. Dann folgt mit den Beziehungen (10.2) aus (7.2) fiir die Komponenten des Spannungs- 
vektors 


Se ae Lees (13) 
und ein Vergleich mit (10.2) zeigt, daB der Schubspannungsvektor S** in jedem Punkt des unter- 
suchten Querschnittes Tangente an die Funktionen F (x1, x?) = konst. ist, diese somit als Schub- 
spannungslinien bezeichnet werden kénnen. Gem48 (13) hei®t die Funktion F(x}, x) dann Span- 
nungsfunktion. 

Bildet man weiterhin mit (7.2) 


EBa see = &B a Gd [p(x*, 20?) | eae c, *]|? 2 
30 folgt wegen der Identitat 
Exp () |"? =0 


EB « Saale =26%. 


die Beziehung 


Daraus erhalt man mit (13) 
B(x, x7)|/% = ARG, A) =— 260 (14) 


ils Differentialgleichung des Torsionsproblems fiir die Spannungsfunktion F(x1,x?). Die Rand- 
bedingung wird nach (11) in der Form 
F (x1, x”) = konst. 
amgeschrieben. Hat der untersuchte Querschnitt nur einen Rand, dann kann die Konstante zweck- 
maBig Null gesetzt werden. 
_ 5. Schnittgréfen des Stabquerschnittes. Die in jedem Querschnitt x? = x? = konst. wirkenden 
Schubspannungen S?* kénnen sowohl zu einem resultierenden Moment M als auch zu einer Quer- 
craft Q! zusammengefaBt werden. Mit den Beziehungen 
: dM = dM(3) = dM? = dM, = dM, 
| df = df(3) = df? = df, = df 
wird fiir das differentielle Moment : 
dM = Eq c x' S36 df 


halten. Unter Beriicksichtigung des aus den beiden krummlinigen Koordinatenabschnitten dx} 
ind dx? gebildeten Flachenteilchens 
df = /g dx dx? 
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ergibt sich das tiber den gesamten Stabquerschnitt f wirkende Torsionsmoment zu 
M =f /g ex, C0 Shae ars 
if 
und mit (7.2) und (13) erhalt man 
M=GPSf Vg [eag ct x (x1, x2)| + x! x,] dx} dx*, 
M = — ff gt xt F(a}, x*)|, dx! dx?. 
yo 


Die an einem Flachenteilchen df, = df = df wirkende Schubkraft dQ (i = 1, 2) 1a8t sich in der 
Form 
dQ’ = Sidf = /g C9" dx dx 
anschreiben. Es folgt 


G=GCOSS Vg lei, ole, 22)|* — e%8 et, ch, x] dx! dx? , : = 
f 


Qi = f ci, xP F(x, x*)|, df = Sf Vgc, ex? F(x, x*)|, dx? dx? . (15) 
f if 
Bei Anwendung des Stokesschen Satzes 


—J eer f= $0 ee 


auf die Beziehung (15) wird 3 
Qi = — f F(x), x?) dx’, 


c 


und da die Spannungsfunktion F(x!, x”) gemaB der geforderten Randbedingung entlang des Quer- 
schnittsrandes C konstant ist, ergibt sich 


iO =20. 


6. Lésung der Torsionsdifferentialgleichung durch Inversion des Spannungszustandes. In kartesi- 
schen Koordinaten lauten die Potentialgleichungen der Torsionsuntersuchung fiir die Verwélbungs- 
funktion (x1, x?) und die Funktion y (z!, x?) 


Ap (#2) = 9@,®)11 + 9H) po = 0, (16.1) 
Ay (2,3) =p 2 )n1 + yp 2 )pps = 0. (16.2) 


Die Funktionen ¢(x!, x?) und y(x1, x?) stellen den Real- bzw. Imaginarteil eines komplexen Tor- 


sionspotentials ®(x! + ix?) dar und sind durch die Cawchy-Riemannschen Differentialgleichungen 


miteinander verkniipft. 


Gelingt es nun, durch die Wahl spezieller krummliniger Koordinatensysteme die Form der | 


Differentialgleichung (16) auch nach ihrer Transformation in allgemeine Koordinaten invariant zu 


gestalten, so kénnen Lésungen von (16) mit Vorteil zur Bestimmung neuer Torsionsquerschnitte — 


herangezogen werden. In diesem Falle hat man lediglich fiir die Variablen x‘ die allgemeinen Koor- 
dinaten x* zu setzen und sichert damit die Erfiillung der Differentialgleichungen 


A(x}; x*). = 0, Ay oh a?) == 05 


Dieser Vorgang wird als Inversion des Spannungszustandes bezeichnet. Ebenso 14Bt sich die In- 
version zwischen krummlinigen Koordinatensystemen untereinander oder von allgemeinen in karte- 
sische Koordinatensysteme ausfiihren. Mit den auf diese Weise erhaltenen Lésungen w(x}, x”) und 
p(x", x") kann nach (12) die Spannungsfunktion F(x!, x2) ermittelt werden. Die Kurvenschar 
F(x1, x*) = konst. liefert dann die Schubspannungslinien bzw. den Rand des Stabquerschnittes. 


Fur eine skalare Funktion f(x", x*) laBt sich die Differentialoperation Af(x1, x2) vorteilhaft in 
der Form 


Af (x22) = f(x, 22)|2 = = [Ve 8°? f(x, 2°)iphe (17) 


angeben. Soll obige Forderung nach Invarianz der Differentialgleichung erfiillt sein, dann diirfen 
nur solche krummlinigen Koordinatensysteme Verwendung finden, deren MaBtensoren den Be- 
dingungen , 


811 = 82> 8i2=—9 (18) 


5 cies LSA tis toe 


{ 


Ue eek Pig eens cee: 


tg ees 


Te ce ae a ane a oe 


) 
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gentigen. Man erhalt mit 
SP ety 2 g,4 
hieraus 


Ves'=Ves*=1, Yeg?=0, 

und damit ergibt sich fiir den Ausdruck 4 f (x1, 2?) 
I 
Af (x1, x?) = Ve. [Ff (x*, ina + f(x, x?)jo0] - 

Nun sind sowohl die Verwélbungsfunktion g(x, x2) als auch die Funktion y(x1, x”) Lésungen der 
Laplaceschen Differentialgleichung, so da® diese auch in allgemeinen Koordinaten in der Form 

p(x, Ca FT + (x1, x") i99 == (Ne (19.1) 

p(x", Coa eT + y(x!, x*)io9 = 0 (19.2) 


angegeben werden kénnen. Damit ist fiir krummlinige, die Bedingungen (18) erfiillende Koordina- 
tensysteme die Invarianz der Laplaceschen Differentialgleichung (16) bei einer Transformation 
nachgewiesen. 


Als Koordinatensysteme, die den Beziehungen (18) geniigen, werden vier Beispiele angefiihit: 


A. Parabolische Koordinaten 


a =e (x1)? panes (x?)? x2 =r x) x2 


2a ; a= 


B. Kreiskoordinaten 


Fil 2 x? =2 2 x? 
a x* = —a : 
(x1)? —- (x?)? > . (x1)2 + (x?)? 5 
C. Elliptische Koordinaten 
x! = acosh x! cos x? , x2 = asinh x! sin x? ; 
D. Bipolarkoordinaten 
1 1 7 2 
ert eee sinh x nee a sin x ; 
cosh x! — cos x2’ cosh x! — cos x? 


7. Beispiele. Fiir die durch die Koordinatenbeziehungen 


pee xt a. Face st eee a 
ee + a?’ GP GP’ 

1 2 xt yee eee ery | eer: 

Se oy Pe eee G+ 


beschriebenen Kreiskoordinaten wird die Lésung des auf der x?-Achse verschobenen Kreises in 
kartesischen Koordinaten invertiert. Sie lautet 


y(t, 2) =e-+ex, (za?) —e + em, 
und man erhalt durch Inversion , 
(xt, x?) = c + Cx, Kee) = c -+- Cx 
als Lésungsfunktionen fiir einen Torsionsquerschnitt in orthogonalen Kreiskoordinaten. Die Span- 
“nungsfunktion ergibt sich dann nach (12) zu 
F(a) = 68 e+ é¢x2 — 


a‘ 


He 
9 (x1)? ae (x2)? 
und kann nach Einfiihrung einer neuen Konstanten 7, einer neuen Funktion K(x1, x?), dimensions- 


5 ; oe i 
loser Koordinaten und eines Festwertes fiir c gema8B 


f a Ca eg : tee 2 F(x!, x?) 1 wet 2 
Cin, 5 eh K(«1, x*) =l——paa > ; gl a’ < a 
in der Form i 
K(é, €7) Peg esta (é?)? =U ) 
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angeschrieben werden. Diese Gleichung liefert fiir verschiedene Parameterwerte 7 eine Schar von = 
Querschnittsrandkurven, fiir die F(x1, x?) = 0 und damit K = 1 gesetzt wird. In Abb. 3 sind-die }} 
Randkurven fiir einige Parameterwerte dargestellt. Fiir7 = 0 erhalt man den Kreis. 


-45 
Abb. 3. Querschnittsrandkurven fir K = 1. Abb. 4. Torsionsquerschnitt fiir 7 = 0,5. 


Hat man aus dieser Kurvenschar eine spezielle Querschnittsform ausgewahlt, dann ergibt (20) 
fiir 7 = konst. die Schubspannungslinien fiir verschiedene Werte von K. Fir 7 = 0,5 ist in Abb. 4 
ein solcher Torsionsquerschnitt mit seinen Schubspannungslinien gezeichnet. Die physikalischen 


Koordinaten des Verschiebungs- und Schubspannungsvektors folgen mit ¢=02 


ws &2 &8 = & & es 7 
V(1) a 3 a» + (&p? V(2) = e0 ap + @P’ V(3) me 4 ge 2 
1 


S31) = 068 ete + TE + ery} 


Be ee ee 
S(82) == 069 aa 


Als zweites Beispiel wird eine Inversion von allgemeinen in kartesische Koordinaten durch- : 
gefiihrt. In parabolischen Koordinaten ist die Lésung y(x!,«*) der Laplaceschen Differential- — 
gleichung zur Beschreibung elliptischer Querschnitte bekannt. Sie lautet a 


p (a a2) = 6 + 6 [at — Gat) (x)? + (x)4], 
und. man erhalt durch Inversion 
ple, 3) = ¢ + e[(e)* — 6)? @)? + @)4]. 


Die Spannungsfunktion kann in der Form 


By 


PENH) = Cle + 1 — OGY + GI] — FIR + OF 


angegeben werden. Setzt man 
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und fihrt ei : . + = 5 
eime neue Konstante 7 und eine neue Funktion K(x, x?) gemaB 


K(z1, x2) =1+- 2 E(a4, x") 
'1—27 Gé@ 


ein, so wird 
(1 — 2m) K(@, ®) — [(@)4— 6&2 (2 + @)4] +. 2 n ((&)? + (@)?] =0. (21) 


In Abb. 5 ist fiir K = 1 [F(«1,x?) = 0] und verschiedene Parameterwerte 7 eine Schar von Kurven 
dargestellt, die als Berandung von Torsionsquerschnitten angesehen werden kénnen. 


92-50 


Abb, 5, Querschnittsrandkurven fir K = 1. Abb. 6. Torsionsquerschnitt fiir 7 = — 2,63. 


/ 


Hin spezieller Stabquerschnitt ist samt seinen Schubspannungslinien fiir 1 — — 2,63 in Abb. 6 
wiedergegeben. Die zu der Funktion (x1, x”) gehérige Verwélbungsfunktion folgt zu 
xX? Fae (xaye 2 
= 1) — +8, 


g(x", 2) = — 
und man erhalt mit ¢ = 0 die Komponenten des Verschiebungs- und Schubspannungsvektors aus 
Vil) =—@028, V2Q)=a2988, 

V3) =<" BB (a — (Oy, 


5(31) = aol (@2— 3&4] —1] B, 


(32) = —aco | 1) — 36 "]—1] 8. 


a 8. Zusammenfassung. Ausgehend von den Grundgleichungen der Elastizitatstheorie in all- 

gemeinen Koordinaten werden durch einen speziellen Ansatz die Grundgleichungen der Torsions- 
theorie abgeleitet. Sowohl fiir eine Verwélbungs- als auch fiir eine Spannungsfunktion werden die 
partiellen Differentialgleichungen aufgestellt und die Randbedingungen angegeben. Fiir spezielle 
krummlinige Koordinatensysteme erhalten die Differentialgleichungen in allgemeinen Koordinaten 
dieselbe Form wie in kartesischen. Daher kénnen bekannte Lésungen von einem Koordinaten- 
system in das andere invertiert werden und liefern neue geschlossene Lésungen des Torsions- 


problems. 
i (Eingegangen am 3. Februar 1960.) 


Anschrift des Verfassers: Dr. H.-J. Franeck, Dresden A 27, Regensburger Str. 9. 
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Klassifikation gedampfter Schwingungssysteme 
und Eingrenzung ihrer Eigenwerte 


Von S. Falk 


1. Einleitung. Die Klassifikation gedampfter Schwingungssysteme mit n Freiheitsgraden ge- 
lingt nach Quade? mit Hilfe der Elementarteiler der Polynommatrix (18). Beschrankt man sich auf 
positiv definite Energieformen, so gibt es fir n = 1 nur drei, fiir n = 2 bereits vierzehn verschie- 
dene Bewegungsvorgange; ihre Anzahl wachst mit n rasch ins Uniibersehbare. In der vorliegenden 
Arbeit wird daher eine Einteilung in nur fiinf verschiedene Typen angegeben, die die eventuell auf- 
tretende Vielfachheit von Eigenwerten unberiicksichtigt laBt. Nebenbei ergeben sich praktisch 
wertvolle Eingrenzungen fiir die Kigenwerte. 


2. Gedimpfte Schwingung mit einem Freiheitsgrad.. Die Bewegungsgleichung des Schwingers 
der Abb. 1 oder eines anderen, ihm mechanisch oder elektrisch Aaquivalenten Systems lautet 
bekanntlich 

axtbxtc=0. (1) 
Dabei ist 


Masse a positiv reell, 


Dampfungszahl 6 , (2) 


Abb. 1. Gedaimpftes Schwingungs- 
system, mit einem Freiheitsgrad. Federzahl c 


beliebig reell. 


Nur a wird als positiv vorausgesetzt. Negatives b bedeutet Anfachung, negatives c wird zum Bei- 
spiel mit c = — a @? bei Vorhandensein von Zentrifugalkraften mechanisch verwirklicht. 
Der Lésungsansatz 


a(t) =e" . (3) 
fihrt auf die charakteristische Gleichung 


a#ti+bA4+c=0 (4) 
oder wegen a + 0 auf 
H+ BA+C=0, (5) 
—#=BA+C (6) 
mit den Quotienten 
b c 
see arg agate (7) 


Die Wurzeln der quadratischen Gleichung (5) sind 


ha=—z + Ve—c. (8) 
Damit lautet die Lésung von (1) 
a) fir komplexe Wurzeln J,,. =a + iw 
~ x(t) = e*' (K, cos wt + Ky sin w t) (9) 
b) fiir reelle Wurzeln A,, A, 
x(t) = K, eh! + K, e%! (10) 
mit beliebigen Integrationskonstanten K, und K,. 


Nach (5) bis (8) gibt es einen ,,Kennpunkt“ P in einer (B, C)-Ebene, der den Bewegungsab- 
lauf eindeutig charakterisiert, und zwar wollen wir nach Abb. 2 (in der das Wort ,,alle“ vorlaufig- 


1 Ing.-Arch. 6 (1935), S. 15. 
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See eae zu ersetzen ist) die folgenden fiinf Bewegungstypen unterscheiden (siehe auch 
a 


| Schwingen Kriechen 


gedampft | Typ I, stabil Typ I, stabil 


Typ V, instabil 


angefacht | Typ III, instabil Typ IV, instabil 
C 


Alle Realyeile Alle Rea/teile 
positiv neganv 


Ypm Ypt 


IN, 


Alle Wurzeln Komplexe Wurzeln Alle Wurzeln 
reel! und postiv reel! und negativ 
"yp Wp 


Alle Wurzeln reell, nicht Null und nicht alle von gleichem Vorzeichen 
lypv 
Abb. 2, Aufteilung der (B, C)-Ebene in fiinf Bereiche. 


G 
S Schwi ich B 
=. chwingungsbere/c C 
\ 
| | | | | | | | | | ‘Stabler Bereich 
B 8 
Abb. 3. Schwingungsbereich und stabiler Bereich in der (BC)-Ebene. Abb. 4. Zeichnerische Ermittlung eines komplexen 


Wurzelpaares mit Hilfe des Kennpunktes P. 


Abb. 5. Zeichnerische Ermittlung eines reellen Wurzelpaares mit Hilfe des Kennpunktes P. 


Eine Bewegung heift stabil, wenn ‘die Koordinate x(t) fiir beliebige Anfangsbedingungen und 
beliebig groBe Werte von t beschrankt bleibt, d. h. also fiir durchweg negative Realteile der Wur- 
 geln A. Wir sagen kurz: 
, ,,Nur der erste Quadrant der (B, C)-Ebene ist stabil.“ ~ 
Die Wurzeln selbst findet man mit Hilfe des Kennpunktes P nach (7) und den Vietaschen Wurzel- 


satzen 
alZ Athi=—D, (11) 
Ai he = Gx (12) 


_ auf folgende Weise: 
30 
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a) Komplexe Wurzeln: der Kreis mit dem Halbmesser /C um den Nullpunkt der komplexen 
Zahlenebene schneidet die Vertikale durch den Punkt « = — B/2 in den beiden gesuchten Werten 
A, und A, (Abb. 4.) 

b) Reelle Wurzeln. Nach (6) schneidet die Gerade y, = BA + C die Parabel y, = —/? in zwei 
Punkten S und T, deren Abszissen die gesuchten Wurzeln A, und A, sind (Abb. 5). 


3. Gedimpfte Schwingungen mit n Freiheitsgraden. Nach diesen Vorbereitungen betrachten 
wir ein beliebiges gedimpftes Schwingungssystem mit n Freiheitsgraden (Abb. 6). Bei hinreichend 


te — 
ml AAA 
UZ, 
es 
-#——> -———_> Pr 
DL, Zp Z3 Zn, 


Abb. 6. Gedaimpftes Schwingungssystem mit n Freiheitsgraden, schematisch, 
In der entspannten Lage des Systems sind alle Koordinaten gleich Null. 


kleinen Auslenkungen lassen sich kinetische und potentielle Energie sowie die ,,Zerstreuungs-— 
funktion“ F der Dampfung darstellen als quadratische Formen der gewahlten Koordinaten x; bzw. 


1: er zeitlichen Ableitungen x;; in Formeln: 
2B == tie 2 a en (13) 


“Vendet man hierauf die Lagrangesche Methode an, so bekommt man das homogene lineare Diffe- 
eatialgleichungssystem zweiter Ordnung mit konstanten Koeffizientent 


Ay +Be+Cr=0 (14) 


mit den drei reellsymmetrischen Matrizen 


Tragheitsmatrix positiv definit , 
Dampfungsmatrix § ! Pen (15) 
i ae © Rang und Signatur beliebig 
und dem Koordinatenvektor 
$= (4, Ree (16) 
Der Ansatz 
u(t) =he', Yy =konst. (17) 
fiihrt auf das lineare homogene Gleichungssystem n-ter Ordnung 
(WH7+ BA4+C)y=0. (18) 


Notwendige Bedingung fiir das Vorhandensein nichttrivialer Lésungen ist das Verschwinden der 
Determinante 

2? + BA + C| = p(s) =0, (19) 
und‘diese Gleichung ist vom Grade 2 n. Das System (18) ist lésbar fiir jede ihrer 2 n Wurzeln, und 
fiir jede einfache von ihnen gibt es genau eine Eigenbewegung von der Form 


u()=hier . (20) 


Die allgemeine Lésung der Bewegungsgleichungen (14) mit 2 n willkiirlichen Integrationskonstanten — 
1aBt sich stets aus 2n Eigenbewegungen linear kombinieren; doch treten bei Mehrfachwurzeln 
von (19) auSer Eigenbewegungen der Form (20) auch solche der Form (30) auf. 


4. Der Wertebereich der Kennpunkte. Um méglichst allgemein zu bleiben, wollen wir die drei 
ate (15) nicht nur als reellsymmetrisch, sondern sogar komplex hermitesch zulassen, d. h. es 
gelte 


YW = positiv definit , 8 =%8, © =€ 3 (21) 


1 R. Zurmiihl, Matrizen 1958, S. 431, Berlin/Gottingen/Heidelberg, 1958. 
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wo Querstriche Ubergang zu konjugiert komplexen Werten bedeuten. Multiplikation der Glei- 
chung (18) von links mit {)’ ergibt 


(2% + BA7A46) pHa 4 bAtc=0. (22) 


Da  definit ist, die zugehérige hermitesche Form a also nicht verschwinden kann, diirfen wir 
dividieren und bekommen 


Be BA C0, ; (23) 
und das ist genau die Gleichung (5) mit den beiden Rayleigh-Quotienten 
y By y Gy 
B= fee Gn = 
y/ Uy y Ay ” ee 


die fiir n = 1 in die KenngréBen (7) tibergehen. Zwar sind uns die Kigenvektoren {);, mit denen 
die Quotienten (24) und damit die Kennpunkte P, der zugehérigen Eigenbewegungen om berechnen 
sind, unbekannt; doch liegen diese Rayleigh-Quotienten 
jedenfalls innerhalb der Wertebereiche der Matrizen- 
paare %, 2 bzw. ©, Y. Diese Wertebereiche sind nach 
den Voraussetzungen (21) reelle Strecken, und es gilt, 
wenn wir die Eigenwerte / und y der Paare %, YU und€, 
Q{ der GréRe nach ordnen, 


py = Brin < B = Bax Se o) (25) 

MSS Cnin = Cc = Ginase Sane (26) 

~ Damit haben wir in der (B, C)-Ebene einen rechteckigen 
Kennbereich, innerhalb dessen die Kennpunkte P; aller 


méglichen Eigenbewegungen liegen miissen (Abb. 7). 
Liegt dieser Kennbereich ganz innerhalb eines der fiinf 


rin 


Gebiete der Abb. 2, so 14Bt sich ohne Kenntnis der Abb. 7. Kennbereich in der (B, C)-Ebene 
; A 3 g fiir positiv definite Matrizen § und ©. 
-Wurzeln / der Typ des Schwingungssystems eindeutig 
_ bestimmen. 


Ehe wir die Abb. 2 aus Abschnitt 2 wbernehmen, miissen wir allerdings noch zeigen, daB beim 
Typ V nicht alle reellen Wurzeln vom gleichen Vorzeichen sein kénnen. Fiir n = 1 war das trivial, 
da ja das Produkt der beiden Wurzeln gleich Cist [Gleichung (12)], und Cist im Bereich V negativ. 
Der Beweis fiir beliebiges n ist umstandlich und sei hier unterdriickt. Keinesfalls kann man ihn 
dadurch erbringen, da man allein die Gleichung (23) heranzieht, die ja auf jeden Fall eine negative 
und eine positive reelle Wurzel besitzt (woraus dann sogar folgen wiirde, und tibrigens bei entkoppel- 
baren Systemen auch tatsachlich folgt, da8 genau n Wurzeln positiv und n negativ sind) ; denn wohl 
hat die mit einem Kigenvektor 4); aufgestellte quadratische Gleichung (23) als Wurzel den zugehéri- 

gen Eigenwert /;, ob aber die andere Wurzel /,; dieser Gleichung ebenfalls Eigenwert ist, steht 
dahin und kann ohne Einblick in die Struktur der Polynommatrix (18) nicht entschieden werden. 


5. Der Wertebereich der Eigenwerte. Die Abb. 4 und 5 zeigten, wie man durch eine einfache 
Konstruktion aus dem Kennpunkt P das zugehérige Wurzelpaar /,,/, bekommt. Man iiberlegt sich 
nun leicht, daB bei gegebenem rechteckigem Kennbereich die gleiche Konstruktion nur fiir die vier 

- Eckpunkte durchzufiihren ist. Damit gewinnt man die folgenden beiden Aussagen: 

< a} Komplexe Wurzeln (C positiv): Der Kreisring mit den Halbmessern V/'Cmin und Cmax 
um den Nullpunkt der komplexen Zahlenebene hat mit dem durch die beiden Abszissen — Bpin/2 
und — Bra,/2 begrenzten vertikalen Streifen ein (eventuell in zwei zur reellen Achse spiegelbildlich 
gelegene Teile zerfallendes) Gebiet gemeinsam, in dem simtliche komplexen Wurzeln der Gleichung 


(19) zu finden sind (Abb. 8). : 
b) Reelle Wurzeln: Die vier Halbgeraden 
yA Bie One | (27) 
schneiden die Parabel y = —/? in vier (oder zwei) reellen Punkten, deren Abszissen E, F; G, H 


(baw. E, F) die beiden Wertebereiche (baw. den Wertebereich) derreellen Wurzeln der Gleichung (19) 


_ begrenzen (Abb. 9). 
30* 
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Nun wird man die Wertebereiche der Matrizenpaare %, 2{ und ©, Q{ nicht immer exakt angeben 
kénnen. Da sich aber bei VergroBerung des Kennbereiches in der (B, C)-Ebene auch der Wurzel- 
bereich in der komplexen Zahlenebene vergroBert, geniigen auch dufere Schranken fiir B und C. 
Solche Schranken findet man zum Beispiel am einfachsten mit Hilfe der Umformung 


y8y 
peewee ee, (28) 
et ee 
yy 
ves 
Go ee Oe (29) 
y AY g 
yy 


denn hier lassen sich die drei Rayleigh-Quotienten a, b und c stets mit minimalem Rechenaufwand 
eingrenzen!. 


Abb. 8. Zeichnerische Ermittlung des komplexen Wurzelbereiches Abb. 9. Zeichnerische Ermittlung des reellen Wurzelbereiches 
mit Hilfe der vier Eckwerte des Kennbereiches, mit Hilfe der vier Eckwerte des Kennbereiches. 


6. Stabilitatsfragen. Zu einem q-fachen Eigenwert A gehért als Verallgemeinerung von (20) 
die Kigenbewegung 


tlt) =e + bitte yee (sSq—l). (30) 


Schon bei einem Freiheitsgrad und reeller Doppelwurzel tritt uns diese Erscheinung in ihrer aller- 
einfachsten Form als sogenannter aperiodischer Grenzfall entgegen. Der Koordinatenvektor (30) 
wachst unbegrenzt an fiir positive Realteile von /, aber auch fiir) = 0 und / rein imaginar. Um also 
ganz sicher zu gehen, schlieBen wir die beiden den ersten Quadranten begrenzenden Halbachsen vom 
stabilen Bereich aus und haben im Hinblick auf Abb. 2 und als Verallgemeinerung der Aussage des 
Abschnittes 2: ; 
»»Der erste Quadrant der (B, C)-Ebene unter AusschluB der beiden begrenzenden Halbachsen ist der 
einzige Bereich, in dem samtliche Eigenbewegungen stabil sind.“ 

Und weiter, da nur fiir negativ definite Federmatrix€ der Rayleigh-Quotient C durchweg negativ 
ist (Typ V): 
»,Bei negativ definiter Federmatrix © sind alle Eigenbewegungen von kriechendem Verhalten, und 
zwar ist mindestens eine von ihnen stabil und mindestens eine instabil.“ 

Dieser Satz ist auf den ersten Blick tiberraschend, leuchtet aber ein fiir n — 1. Hier iibersieht 
man sofort, da der Massenpunkt im Zentrifugalkraftfeld in den Nullpunkten hineinkriechen und 
dort fiir t > co liegen bleiben kann, einerlei, ob die Bewegung angefacht oder gedimpft ist. 


7. Feder- und Dimpfungsmatrix sind positiv definit. Der praktisch weitaus wichtigste Fall ist 
der positiv definiter Matrizen § und ©. (Q{ ist stets positiv definit.) Der Kennbereich liegt jetzt 


1 R. Zurmiihl, Matrizen, S. 203—210. 
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“ae im ersten Quadranten (Abb. 7), es kommen allein die Typen I und II in Betracht, alle Eigen- 
ewegungen sind stabil, eine bekannte Tatsache, die hiermit nochmals bewiesen ist. Weiterhin gilt; 


Ist 4C,, i : a 
s Can DE Be ioci Typ I vor (schwache Dampfung) (31) 
tet. AC 


= min 


Typ IL vor (starke Dampfung) 


Trifft keine der beiden Bedingungen zu, so sind Eigenbewegungen vom Typ I (gedampfte Schwin- 
gungen) oder vom Typ II (Kriechbewegungen) oder von beiden Typen méglich. 


Abb. 10. Ermittlung des reellen Wurzelbereiches Abb. 11. Eingrenzung der Eigenwerte 
fiir positiv definite Matrizen § und ©. mit Hilfe der Abschatzung (35). 


Der Wertebereich fiir etwa vorhandene reelle Wurzeln besteht nach Abb. 10 nur aus dem einen 
Stiick E F; er wird begrenzt durch die beiden leicht zu berechnenden Werte 


AEF =— Brax/2 a / B2ax/4 a. Gan . é (32) 


-Wegen der vorausgesetzten positiven Definitheit aller drei Matrizen ist Null eine untere Schranke fiir 
_ Bund C. Obere Schranken sind ebenfalls leicht angebbar, wenn 2{ von Diagonalform ist (was in fast 


allen praktischen Fallen zutrifft). Dann namlich bilden wir zunachst die Zeilen — Betragssummen 


der Matrizen § und © 
b; =) |b; ’ é;= Dd} le;x| (33) 
k= k=l 
und daraus die stets positiven reellen Quotienten 


pee ge. SU. (34) 


Der EinschlieBungssatz von Gerschgorin liefert dann fiir die beiden Rayleigh-Quotienten (24) 


" 


die Abschatzung 
. Cie bien 0 Oe, bbc (35) 


Der zugehérige Wurzelbereich besteht aus einem Halbkreis der komplexen Zahlenebene und der 
Strecke EF auf der negativen reellen Achse (Abb. 11). 


8. Verallgemeinerung auf beliebige Matrizen. Jede quadratische Matrix R laft sich darstellen 
in der Form 


R=K +R, (36) 


wo die Anteile §{ und i beide hermitesch sind. Die quadratische Gleichung (22) lautet daher bei 
beliebigen (reellen oder komplexen) Matrizen 2, 6 und © 


(a+i@)et(6+ib)a+(e+id) =o. (37) 


In dieser Gleichung treten sechs reelle hermitesche Formen auf; mit Hilfe ihrer Wertebereiche (oder 
wenigstens duGerer Schranken fiir sie) lassen sich auch jetzt wieder alle 2 n Wurzeln in Gebiete ein- 
schlieBen, doch gehen wir nicht darauf ein. 
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Besonders einfach ist noch der Fall, daB 2 und © hermitesch, % dagegen schiefhermitesch ist 
(Kreiselwirkungen). Dann wird zunachst aus (23) wegen 1 = —# : 
—?/2—#PAB4C=0. (38) 


Jetztisti B= B reell; setzen wir daher iA =, so geht (38) tber in 
: go ge pecs (39) 


und diese Gleichung wird wieder wie friher behandelt. 


9. Beispiele. a) Erstes Beispiel. Gegeben ist das gediampfte Schwingungssystem von Abb. 12 
mit n = 3 Freiheitsgraden. Die Matrizen (15) lassen sich hier direkt aus der Zeichnung ablesen; es 


ist 


fa 0-0 {Fees | Se | 2—l 0 
Tmto 2 0), BHdlo- 3 —11,. © =e) —1] sae 
0 20% 0—l 6 0-07 10 


0c 


Abb. 12, Gedampftes Schwingungssystem mit drei Freiheitsgraden, (Beispiel 1). 


Alle drei Matrizen sind positiv definit, 2{ hat Diagonalform; also kénnen wir die rohe Alschateuas 
(35) benutzen. Zunachst bilden wir die Zeilen-Betragssummen nach (33) 


ae Ie Cy = an; 
| pals” ee 
daraus die Quotienten (34) 

= 0 = 

=> =0, S23, 
ef 4 ez 
==, => =2, 
= ai = 

oe 1h, Ga 2 


und haben damit nach (35) die Aussage 
(eo Oe ees 


Die zugehdrigen Bereiche sind in Abb. 13 und 14 gestrichelt eingezeichnet. 
Wesentlich engere Grenzen bekommen wir mit den exakten Eigenwerten f; bzw. y; der Matrizen- 
paare %, Whzw. ©, Y, die sich hier leicht berechnen lassen. Man findet (in dimensionslosen GréBen) 


8 —p U| = (0 —f) (og*—278 +17) =0, 
pb, =9, fg=1, Bgl 
Boe ==0% ape ree ss Liles 
|E —y U| = (2y2?— Ty +5) (10—5y) =0, 
Yy=1, Yy=2, Naas a 


das hei®t nach (25) 


somit nach (26) 
Cas = 9 Gras = 259) . 
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Den Kennbereich in der (B, C)-Ebene zeigt Abb. 13; es liegt der Typ I vor, alle Eigenbewegungen 
sind stabil (starke Dampfung). Der zugehérige Wurzelbereich besteht aus einem zerfallendem 


Kreisringgebiet (Abb. 14). 


47 2,0 


Abb. 13. Kennbereich und Kennpunkte fiir das 
Schwingungssystem von Abb, 12. 


Abb. 14 (rechts). Eigenwertbereich und Eigen- 
werte fiir das Schwingungssystem von Abb. 12. 


Zum Vergleich sind auch die exakten Werte aus der Gleichung (19) berechnet und in die Werte- 
bereiche eingetragen worden, namlich aus 


p(A) = 10 AS + 27 45 + 72 44 + 102 23 + 129 22 + 904 + 50 = 0 


+ 


B, [a/m] 


| 
| 
| —0,753 + 0,858 i 


1 1,506 1,303 
2 —0,514-£ 1,247i | 1,028 1,819 
3 —0,083 £1,4507 | 0,166 2,109 


Abb. 15. Kennbereich fiir das Schwingungssystem Beispiel 2. 


Abb. 16 (rechts). Eigenwertbereich fiir das Beispiel 2. 


b) Zweites Beispiel: Y{ und © positiv definiert, 8 indefinit, n beliebig. Der Kennbereich sei 


-bekannt: 
Dis ae 1,00 ? Crtin a 0,81 ? \. Gea = 0,90 ° 


Dinax = + 4,60 , Cmax = 2,25, i Coe = 1,50. 


2S Hand von Abb. 15 stellen wir fest, da die Typen I, II und HI in Betracht kommen. Der 
komplexe Wurzelbereich ist nun leicht zu zeichnen (Abb. 16), den reellen Bereich berechnen wir nach 


(32) Bee ea aes 
: a + 2,3? — 0,81 , 


Ap = — 4,92, Arp = —9,18. 


eas 
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c) Drittes Beispiel: Schwingungssystem mit n Freiheitsgraden nach Abb. 17. Die 2 n Eigen- 
werte A; = + iw; des ungedampften Systems seien bekannt. Die Diampfungen werden nachtraglich 


angebracht. 
Frage 1: Wie stark diirfen die d; héchstens sein, damit noch alle Eigenbewegungen zum Typ I 
gehéren (gedampfte Schwingungen) ? 
Frage 2: Wie stark miissen die d, mindestens sein, damit alle Eigenbewegungen zum Typ I 


gehéren (Kriechbewegungen) ? 


Abb. 17. Gedimpftes Schwingungssystem mit n Freiheitsgraden, 2{ und & von Diagonalform, Beispiel 3. 


In den gewahlten Koordinaten sind die Matrizen 2f und % von Diagonalform, also sind ihre 
Kigenwerte 


Bb; = O:/7,. Smit sb. =(5 , ies =(3] A 


Nach (31) muB sein 
4 Crate = Br as r) 
das hei®Bt hier 


2 rs d ? 
4 Mmin = ae > 
also als Antwort auf Frage 1: 
d 
2 Omin = i 


max 


Die Antwort auf Frage 2 lautet dementsprechend 


d 
: sae = ele : 


(Eingegangen am 16. Februar 1960.) 


Anschrift des Verfassers: Dr.-Ing. S. Falk, Braunschweig, Techn. Hochschule, Institut fiir Technische Mechanik 
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